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INTRODUCTION 



Dans les annees 50, lorsque Bardeen, Cooper et Schrieffer ont developpe 
une theorie pour expliquer la supraconductivite (theorie BCS), ils ont montre que 
celle-ci permettait d'expliquer l'effet Meissner observe lorsque Ton soumet un 
supraconducteur a un champ electromagnetique [8]. Leur derivation fut effectuee 
en utilisant la jauge de Landau. Cependant, lorsque leur calcul etait reproduit 
dans une autre jauge, on obtenait une distribution de courant differente ne 
respectant pas la conservation du nombre de charge [94]. En d'autres mots, leur 
theorie brisait l'invariance de jauge. Ce probleme fut resolu simultanement par 
Anderson [4] et par Bogoliubov et al. [15], puis de fagon plus complete par 
Rickayzen [93]. L'explication repose sur le fait qu'il faut inclure dans le calcul des 
courants la contribution provenant des modes collectifs longitudinaux afin de 
restaurer l'invariance de jauge. Sans cette contribution, les parties longitudinales 
des courants paramagnetiques et diamagnetiques ne se compensent pas lors d'un 
changement de jauge. 

Suite a la solution de ce probleme, Baym et Kadanoff [10,11] ont propose 
un formalisme permettant d'assurer qu'une solution meme approximative d'un 
modele a N-corps puisse satisfaire les lois de conservation globales du nombre de 
particule (invariance de jauge), de Penergie et de la quantite de mouvement. Ce 
formalisme consiste a assurer une certaine coherence entre la description des 
excitations elementaires du systeme et ses proprietes de transport. 

Encore de nos jours, beaucoup de recherches ont pour objet le 
developpement de nouvelles methodes analytiques visant a predire les proprietes 
de modeles relativement simples comme le modele de Hubbard. L'interet pour ce 
domaine s'est accentue depuis la decouverte de materiaux a proprietes exotiques. 
La forte anisotropie de la plupart de ces systemes a motive beaucoup d'etudes 
sur les systemes unidimensionnels et bidimensionnels. En particulier, pour les 
systemes quasi-unidimensionnels, plusieurs techniques telles la bosonisation [36], le 
groupe de renormalisation [16], la theorie des champs conformes [39], l'ansatz de 
Bethe [120] ont permis de donner une description detaillee de l'etat fondamental, 
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des etats excites et des fonctions de correlation de differents modeles et ce pour 
des cas a forte comme a faible interaction. Plus recemment, la theorie du champ 
moyen dynamique [44] a permis de resoudre exactement des modeles a dimension 
infinie. 

Pour l'etude du cas bidimensionnel, a faible couplage la phase ordonnee est 
assez bien decrite par une theorie de type BCS. De meme, la limite a fort 
couplage (limite atomique) du modele de Hubbard est assez bien comprise et son 
comportement critique correspond a la condensation de bosons 1 [84]. Cependant, 
en 1' absence de theories couvrant les deux limites beaucoup de recherches ont 
ete menees dans le domaine [79]. En plus des tentatives pour generaliser les 
methodes fructueuses a une dimension, plusieurs techniques nouvelles ont ete 
developpees: approche par les moments (approximation de type n-poles) [100], 
approximation dynamique par amas [50], approche fonctionnelle [47], approche 
GW (Approximation auto-coherente pour la fonction de Green, G, faisant appel 
au potentiel ecrante, W) [5], approximation de parquet [35], approximation FLEX 
("Fluctuation exchange" ) [45] , approximation TPSC ("Two-particle self-consistent") 
[116,117]. On a vu aussi le developpement d'approches diagrammatiques a partir 
de la limite fort couplage [86]. 

Dans cette these, nous developpons une generalisation de 1' approximation 
TPSC [116,117] qui permet de traiter le modele de Hubbard attractif. Parmi 
les methodes mentionnees precedemment, celle-ci semble etre Tune des plus 
prometteuses. En effet des comparaisons avec les donnees Monte Carlo pour de 
petits systemes ont donne des resultats tres satisfaisants [80]. En particulier, 
1' application a des regions du diagramme de phase peu sujettes aux effets de 
taille a mis en evidence la presence d'un pseudogap dans le modele de Hubbard 
repulsif [116,80]. 

1 Dans la limite fort couplage, les electrons forment des paires locales a haute 
temperature. A la temperature critique, on observe une condensation des bosons dans 
l'etat fondamental appelee condensation de Bose-Einstein. 
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L'approche est basee entre autres sur des regies de somme qui permettent 
d'assurer une coherence entre un vertex irreductible approximatif et les fonctions 
de correlation a deux particules, d'ou le nom de l'approche: approximation 
auto-coherente a deux particules ("Two-particle self-consistent"). Par la suite, on 
peut evaluer la contribution de ces fluctuations a la self-energie. Ainsi, par 
l'application de cette methode nous pourrons obtenir les proprietes a une 
particule, dont le poids spectral. Cette methode sera appliquee entre autres a une 
etude detaillee du probleme du pseudogap dans les supraconducteurs a haute 
temperature critique (SC hTc). II y a deja un grand nombre de methodes qui ont 
ete appliquees a ce probleme, la plupart (matrice-T auto-coherente, FLEX) ne 
predisant pas de pseudogap pour le modele de Hubbard. 

Cette these est divisee de la fagon suivante. Dans le premier chapitre, 
nous expliquons comment notre etude du modele de Hubbard attractif en deux 
dimensions pourrait aider a comprendre la presence du pseudogap observe dans 
les SC hTc. Nous abordons l'effet des fluctuations supraconductrices sur le poids 
spectral des excitations fermioniques. Afin d'expliquer la presence d'un large 
domaine de temperature oil un pseudogap pourrait etre present, nous mettons 
l'accent sur l'effet de l'augmentation de la symetrie du parametre d'ordre. 
L'existence d'un tel mecanisme dans les SC hTc a ete predite par Zhang avec sa 
theorie SO (5) [122]. Cette theorie suppose que Pelargissement de la symetrie du 
parametre d'ordre serait present entre la region avec ordre antiferromagnetique et 
la region avec ordre supraconducteur la oil le pseudogap est observe sur un vaste 
domaine de temperature. 

Le chapitre suivant constitue le coeur de la these. Nous y decrivons 
l'approche fonctionnelle developpee. Notre point de depart est un ensemble de 
relations exactes: l'equation du mouvement de la fonction de Green, l'equation 
de Dyson et l'equation de Bethe-Salpeter. Etant donne le formalisme matriciel 
propose ici, nous rederivons ces formules qui prennent des aspects quelque peu 
differents de la forme habituelle. Par la suite, nous suivons le formalisme de 
Baym et Kadanoff qui nous permet d'obtenir le vertex irreductible a partir d'une 
forme approximative pour la fonctionnelle de self-energie. Le systeme d'equations, 
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ainsi etabli, peut etre ferme par l'imposition d'une regie de somme deduite 
du theoreme de fluctuation-dissipation. Ainsi, notre theorie donne une solution 
approximative du modele de Hubbard sans l'ajout d'aucun parametre arbitraire. 
De plus, nous demontrons que notre approximation reproduit exactement le 
premier moment de la susceptibilite de paires et, de fagon auto-coherente, le 
second moment. 

Finalement, dans le troisieme chapitre, nous verifions la validite de notre 
approximation et comparons nos resultats aux donnees Monte Carlo. Nous 
montrons que notre theorie prevoit la presence d'un pseudogap dans le poids 
spectral pres du demi-remplissage et ce pour des temperatures bien superieures a 
la temperature critique. Nous montrons que 1' apparition de ce pseudogap est 
associee a la presence d'un regime classique renormalise. Ces resultats sont 
obtenus dans un regime ou les effets de la taille finie du systeme sont controles 
et ce meme pour les tailles accessibles par simulation Monte Carlo. 

Afin d'aider a la comprehension de certains passages de la these, 
nous avons ajoute quelques annexes a la fin de celle-ci. Par exemple, dans 
la premiere annexe, le lecteur trouvera une table de formules utiles pour 
l'application du formalisme presente au deuxieme chapitre. Egalement en annexe, 
nous avons donne une description de la technique de prolongement analytique 
par maximisation d'entropie. Cette technique est la methode la plus efficace 
permettant d'extraire le poids spectral de la fonction de Green en temps 
imaginaire mesuree par les simulations Monte Carlo. 
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MOTIVATION ET MISE EN CONTEXTE 

1.1 La supraconductivite a haute temperature 

Depuis la decouverte par Bednorz et Miiller [12] de la presence d'une phase 
supraconductrice a haute temperature dans le melange Ba, E La 5 _, E Cu 5 5 (3_, / - ) , on 
a decouvert toute une famille d'oxydes de cuivre presentant cette meme 
phase. Ces composes sont caracterises par une forte conductivite dans les 
plans cuivre-oxygene (Cu0 2 ) et une faible conductivite entre ces plans. Cette 
anisotropic explique pourquoi on considere ces composes comme des systemes 
quasi-bidimensionnels 1 . 

Pour l'etude du diagramme de phase, la variation de la densite de 
porteurs dans ces plans se fait par subtitution d'atomes hors des plans ou 
par modification du nombre d'atomes d'oxygene apical [59]. Typiquement, le 
diagramme de phase est caracterise par une phase antiferromagnetique (AFM) 
dans la region generalement considered comme faiblement dopee (Nd 2 _ ;c Ce :E Cu04 
avec x < 0.13 La2_2.Sr2.CuO4 avec x < 0.02). A plus fort dopage, assez pres de la 
disparition de la phase AFM, se trouve une phase supraconductrice (voir figure 
1). La region entre ces deux phases est differente d'un compose a un autre. 
Certains presentent une phase de verres de spin, d'autres des rayures ("stripes") 
ou se presente une alternance de regions avec ordre local AFM et de regions 
riches en porteurs de charges [119]. 

Une caracteristique commune au differents composes est la presence d'un 
dopage ou la temperature critique entre la phase "normale" et la phase 
supraconductrice est maximale. Ce dopage est appele dopage optimale. La section 

1 La presence d'une interaction entre les plans cuivre-oxygene est toutefois impor- 
tante car elle determine la valeur de la temperature critique 
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comprenant la phase AFM jusqu'au dopage optimale est appelee le regime 
sous-dopee et l'autre section le regime surdopee. 

Une autre caracteristique du diagramme de phase est la presence de 
proprietes anormales pour un vaste domaine de temperatures superieures a la 
temperature critique. Ces proprietes sont dites anormales car elles ne peuvent etre 
decrites par la theorie des liquides de Fermi developpee par Landau. Une des 
caracteristiques ayant attire beaucoup d'attention est la presence d'un abaissement 
de la densite d'etats au niveau de Fermi. Sous la temperature T* identified sur la 
figure 1, l'abaissement est tel qu'on observe la presence d'un minimum local a ce 
point. C'est ce que Ton appelle le pseudogap (PG). Ce phenomene a ete observe 
dans un grand nombre d' experiences: mesures de chaleur specifique [69], de 
conductivity optique [25], resonance magnetique nucleaire [106], diffusion Raman 
[97], diffusion inelastique de neutrons [28], spectroscopic resolue en temps [33], 
spectroscopie par effet tunnel [92], photoemission [34] (voir Timusk et Statt pour 
une revue des differents resultats experimentaux [107]). 

Le diagramme de phase presente a la figure 1 est schematique et change 
legerement d'un compose a un autre. La presence du pseudogap a ete observe 
dans plusieurs de ces composes que ce soit ceux de la famille des YBaCuO, de la 
famille des lanthanes ou des bismuths. Le pseudogap a ete observe au-dessus de la 
phase AFM aussi bien qu'au-dessus de la phase supraconductrice. Typiquement, la 
temperature ou celui-ci apparait diminue au fur et a mesure que Ton s'eloigne de 
la phase AFM, pour tendre vers la temperature critique dans la phase surdopee. 

Notre etude du modele attractif visera a expliquer l'existence de ce regime 
avec pseudogap. Mais avant de presenter notre approche, nous exposons quelques 
theories proposees pour expliquer ce pseudogap. La plupart de ces theories sont 
presentees dans Particle de revue de Loktev et al. [68]. 
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1.2 Pseudogap 

Pour la suite de ce travail, nous defmirons le regime pseudogap par la 
region du diagramme de phase ou on observe dans le poids spectral a une 
particule un minimum local au niveau de Fermi (sans toutefois qu'il n'y ait un 
vrai gap, c'est-a-dire au-dessus de la temperature critique). D'un point de vue 
quantitatif le pseudogap sera donne par la mesure de la separation en energie des 
maxima locaux situes de part et d'autre du niveau de Fermi. 

Partant de l'observation que le pseudogap evolue de fagon continue, 
devenant le gap dans l'etat ordonne (pour un dopage au-dessus de la 
phase supraconductrice) , plusieurs theories sont basees sur l'existence d'un 
mecanisme permettant d'abaisser la temperature critique sans toutefois modifier la 
temperature a laquelle apparaissent les fluctuations supraconductrices. Le gap dans 
les excitations fermioniques ne se forme que dans la phase ordonnee, done a 
basse temperature. L'ecart important entre ces deux temperatures laisse done 
place a un regime ou des proprietes anormales peuvent se presenter. Pour mieux 
identifier ces proprietes, il faut comprendre quel serait le mecanisme responsable 
de l'abaissement de la temperature critique. Une deuxieme question importante a 
laquelle on doit repondre concerne le developpement du pseudogap dans ce 
regime. Jusqu'a maintenant, plusieurs theories ont propose differentes reponses a 
ces questions. 

Kivelson et Emery [37] furent parmi les pionniers de l'etude theorique du 
pseudogap avec leur scenario base sur l'effet des fluctuations supraconductrices. 
lis ont propose que la faible densite de porteurs dans les supraconducteurs 
sous-dopes donne un role important aux fluctuations de phase et provoque 
l'apparition du pseudogap. Supposant que le comportement critique appartient a 
la classe du modele XY, leur theorie a obtenu un accord quantitatif interessant 
avec un certain nombre d'observations. 

Tout recemment, Nozieres et Pistolesi ont propose une theorie faisant 
aussi appel a une faible densite de porteurs [85]. lis ont etudie la possibility 

7 



Motivation et mise en contexte 



qu'un semiconducteur puisse avoir une transition vers un etat supraconducteur. 
Le pseudogap resulte alors d'une competition entre le gap de la phase 
supraconductrice et celui associe a la phase semiconductrice. 

Par ailleurs, Randeria a, tout comme Emery et Kivelson, fait appel a 
l'observation d'une faible rigidite du superfluide dans l'etat ordonne pour proposer 
une nouvelle theorie [91]. Cependant, il a note le fait que pour le modele XY, la 
temperature critique champ moyen et la transition Berezinskii-Kosterlitz-Thouless 
(BKT) sont tres rapprochees et qu'un tel modele ne peut done pas presenter une 
vaste region avec un pseudogap. Dans sa theorie, le pseudogap est un vestige de 
la phase supraconductrice et provient de l'existence de paires a courte portee 
au-dessus de la temperature critique. Cette theorie requiert done l'existence d'une 
interaction assez forte entre les particules. Comme les autres theories faisant appel 
a la force de l'interaction, la dimension du systeme ne joue pas, ici, de role 
primordial pour l'existence du pseudogap. 

La grande difference entre la temperature critique et la temperature ou 
apparait le comportement anormal, T cr ~ 6T C [27], a aussi ete utilisee par Pines 
pour critiquer la theorie de Emery et Kivelson [88]. II a propose avec ses 
collaborateurs, un scenario ou e'est la formation de precurseurs d'ondes de 
densite de spin qui induit un pseudogap dans les excitations fermioniques pres 
des singularites de Van Hove ("hot spots"). Cette theorie est une extension de la 
theorie des liquides de Fermi quasi-antiferromagnetiques developpee par le groupe 
de Pines. 

On ne peut passer sous silence la theorie du groupe de Levin [65]. 
Leur theorie consiste a faire une distinction entre le gap dans les excitations 
fermioniques et le parametre d'ordre. lis associent la difference entre ces 
parametres a la presence d'excitations de paires de faible energie. Au-dessus de la 
temperature critique, le parametre d'ordre est nul, mais l'existence de paires 
preformees pourrait faire appaitre le pseudogap. Un prealable important de leur 
theorie est l'existence d'une resonance dans la fonction de correlation de paires. 
Cependant, la creation de cette resonance par un processus microscopique n'a 
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pas ete demontree. Le groupe de Levin l'obtient par une approximation de 
type matrice-T ou les particules couplees par l'interaction ont des niveaux de 
Fermi differents. Cette difference etant liee a la force du couplage, 5fi oc E, ils 
obtiennent qu'au-dessus d'une valeur critique un pic de resonance apparait dans 
la susceptibilite de paires. C'est cette resonance qui permet l'apparition d'un 
pseudogap au-dessus de la temperature critique. La taille de ce dernier est 
proportionnelle a l'energie de la resonance et done a la difference des niveaux de 
Fermi: S oc A 2 g . II ne faut cependant pas trop pousser la critique de cette theorie 
et garder un esprit ouvert. Certains effets, tel la presence du niveau de Fermi hors 
de la bande d'energie permise, pourraient amener l'apparition d'une resonance. 

En comparaison, la theorie que nous proposons est basee sur la presence de 
fluctuations de phase supraconductrices. Cependant, nous proposons que, dans la 
region sous-dopee, ces fluctuations n'appartiennent pas a la classe BKT tel que 
propose par Loktev et al. [68] mais a une classe a symetrie plus elevee. 

La theorie SO (5) de Zhang [122] fournit un exemple de theorie ou la 
symetrie joue un role important dans les SC hTc. Partant de la proximite entre 
la phase supraconductrice et la phase antiferromagnetique, Zhang a suggere 
que dans la partie du diagramme entre ces deux phases, le parametre d'ordre 
possederait une symetrie approximative de type SO (5). Par ailleurs, il est bien 
connu qu'une augmentation de la symetrie du parametre d'ordre a pour effet de 
faciliter la presence de fluctuations de phase 2 et done d'abaisser la temperature 
critique (particulierement pour un systeme bidimensionnel) . Ainsi, il devient 
possible d'observer une vaste region de fluctuations de phase au-dessus de la 

2 A basse symetrie, deux configurations du parametre d'ordre minimisant l'energie 
peuvent etre topologiquement distinctes. II faut alors surmonter un cout energetique 
important pour passer de l'une a l'autre. Une telle transformation est done improba- 
ble. L'augmentation de la symetrie du parametre d'ordre amene de nouvelles avenues 
par lesquelles nous pouvons transformer les configurations. II s'ensuit que deux con- 
figurations qui etaient topologiquement distinctes a basse symetrie peuvent devenir 
topologiquement equivalentes lorsque la symetrie est augmentee. 
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temperature critique dans la region a symetrie elevee. Ce sont ces fluctuations de 
phase qui seraient responsables de l'apparition du pseudogap. 

II est a noter que meme si la presence ou l'absence d'une symetrie 
supplementaire devrait etre detectable sans ambiguite, l'effet que Ton veut mettre 
en lumiere pourrait se presenter dans un cas ou la symetrie n'est pas exacte. 
Ainsi, l'apparition de cet effet ne se ferait pas de fagon brusque mais serait 
plutot quelque chose de continu. Illustrons ce que Ton veut dire en considerant 
l'exemple du modele XXZ avec un terme d'anisotropie, A, faible [26]. A une 
certaine temperature ce modele presente le meme comportement critique que le 
modele de Heisenberg, comme si le terme d'anisotropie n'etait pas present. Ce 
n'est qu'a plus basse temperature que le systeme suit un crossover vers le 
comportement critique du modele XY. Plus le terme d'anisotropie est faible plus 
le crossover se fera pres de la temperature critique. 

II pourrait en etre de meme pour les SC hTc. A temperature elevee des 
fluctuations a symetrie elevee apparaissent. Cette forte symetrie implique que la 
phase du parametre d'ordre ne peut s'ordonner qu'a basse temperature. II peut 
alors subsister un large domaine de temperature ou, sous l'effet des fluctuations, 
le systeme presenterait des proprietes anormales. Cet effet serait accentue par la 
quasi-bidimensionalite du systeme. Ainsi, le pseudogap qui serait associe a des 
fluctuations a symetrie elevee, pourrait se developper jusqu'a une temperature ou 
Ton observerait un crossover vers un regime critique a plus basse symetrie. Ce 
crossover permettrait d'expliquer revolution continue du pseudogap vers le gap 
qui est associe a un parametre d'ordre a symetrie reduite, celui de la phase 
supraconductrice . 

Ann de verifier cette theorie, nous avons choisi d'etudier le modele 
de Hubbard attractif a deux dimensions. II est clair que ce modele n'est 
pas un modele realiste pour decrire la physique des supraconducteurs a 
haute temperature critique. Cependant, la presence de regions avec differentes 
symetries pour le parametre d'ordre de ce modele nous permettra de verifier si 
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l'elargissement de la symetrie du parametre d'ordre pourrait etre un mecanisme 
responsable de l'agrandissement du domaine ou un pseudogap apparait. 

Apres avoir discute du modele de Hubbard dans la section suivante, nous 
reviendrons sur la question de l'existence du pseudogap dans ce modele. 



1.3 Modele de Hubbard attractif 

Le modele de Hubbard sur reseau decrit par la relation (1.1), presente 
deux parametres: le terme de saut, t, qui est l'amplitude de probability qu'un 
electrons passe d'un site du reseau a un autre site et l'interation entre des 
electrons se trouvant sur un meme site, U. A ces parametres s'ajoutent ceux de 
1 'ensemble grand canonique que nous utiliserons: la temperature, T, et le potentiel 
chimique, fi, qui sont les parametres de Lagrange associes respectivement a 
l'energie moyenne et au nombre moyen de particules du systeme. En plus de ces 
parametres, s'ajoute une liberte dans le choix du reseau. Ann de demontrer 
l'argumentation avancee, nous choississons un reseau qui permet d'assurer la 
variation de la symetrie du parametre d'ordre. Nous considererons done le cas du 
reseau carre bidimensionnel. Nous etudierons le cas ou l'integrale de saut est non 
nulle que pour les sites voisins, £ = t(5 lJ ± s . + 5 l 3 ±$)- Ce choix nous permettra 
d'utiliser la propriete bipartite du reseau. II est a noter que les indices t et j 
indiquent un site du reseau, quant aux vecteurs x et y, ils indiquent les vecteurs 
de la base du reseau cristallin. De plus, la condition que i = t Jt est une 
condition essentielle pour l'hermiticite du Hamiltonien. Comme unite de mesure 
pour l'energie nous utiliserons le terme de saut, t. 




(1.1) 



Hamiltonien du modele de Hubbard 



Le diagramme de phase du modele etudie dependra done de trois 
parametres. De fagon schematique, il sera donne par les figures 2 et 3. Pour 
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mettre en evidence la ressemblance avec celui des supraconducteurs de type oxyde 
de cuivre, nous l'avons represents en fonction du dopage x = 1 — n qui represente 
l'ecart entre un remplissage donne, n, et le cas du demi-remplissage, n = 1. Nous 
avons represente sur la figure 2, la temperature critique, T BKT , la temperature 
critique predite par un calcul champ moyen, T°, la temperature a laquelle 
apparaitrait le pseudogap, T* ainsi que la temperature oil Ton observerait un 
crossover entre le comportement critique a haute symetrie et le comportement 
critique 0(2). Alors que la temperature T° peut etre evalue numeriquement, la 
temperature critique T BKT n'est pas connue. Ici, nous avons represente ces deux 
courbes de fagon schematique. Le crossover quant a lui n'est pas une ligne 
bien definie. Sa representation est egalement approximative. Etant donne que 
le crossover est associe au remplissage, la temperature oil celui-ci se produit 
augmente en s'eloignant du cas demi-rempli (x = 0). Pour un remplissage plus 
grand que le remplissage optimal, cette temperature decroit car alors le 
comportement critique ne se manifeste qu'a basse temperature. Enfin, l'existence 
d'une temperature oil apparait un pseudogap n'est que speculative et nous allons 
tenter de montrer son existence. 

Nous avons donne egalement le diagramme de phase en fonction de 
l'interaction, figure 3. On y a represente les deux temperatures dont on connait 
revolution: la temperature critique, T BKT qui augmente de fagon exponentielle a 
faible interaction puis passe par un maximum pour decroitre en suivant une loi 
de puissance a forte interaction; et la temperature critique champ moyen, T°, qui 
augmente de fagon exponentielle a faible interaction puis tend vers une constante 
a forte interaction. Les deux courbes sont representees de fagon schematique. 
Physiquement, la temperature champ moyen correspond a la formation de paires 
de Cooper. Lorsque le couplage est faible celle-ci est a peu pres identique a la 
temperature critique T BKT . Cependant a fort couplage, cette temperature est bien 
differente de la temperature critique. Cette derniere correspond a la temperature 
oil se fait la condensation des paires locales (condensation de Bose-Einstein). 
Etant donne que l'energie de liaison des paires est proportionnelle a la force du 
couplage, il en sera de meme pour la masse des paires. Du fait que des bosons 
de grande masse ne peuvent condenser qu'a basse temperature ceci explique 
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pourquoi la temperature critique decroit (T BKT oc 1/|C/|) a fort couplage [84]. 
Ainsi dans cette region on observe une grande difference entre T BKT et T°. 
Cette difference explique pourquoi differentes theories du pseudogap font appel 
a la presence d'une forte interaction. Nous avons identifie sur le graphique 
la valeur de couplage oil se trouve la valeur maximale de la temperature 
critique telle qu'estimee par Singer et al. [104: figure 3] (Estimation par simulation 
Monte Carlo). Dans notre travail, nous utiliserons une valeur de U = —4. Celle-ci 
est situee sur dans la region oil la temperature critique croit en fonction du 
couplage. Ainsi, le systeme sera dans la limite couplage faible a inter mediaire. 

II est possible de montrer qu'a demi-remplissage notre modele presente 
une symetrie plus elevee qu'en dehors du demi-remplissage. La particularite du 
cas demi-rempli vient du potentiel chimique qui est alors connu de fagon 
exacte, /i = U/2. Cette valeur particuliere de \i permet d'avoir une symetrie 
supplemental. Pour le demontrer considerons la transformation canonique 
suivante: 

Ai - (-i)Vi < - (-i)Va (i-2) 

Considerant la transformation de chacun des termes du hamiltonien, on peut 
montrer comme resultat intermediaire que: 

De ces resultats on obtient que le terme de saut aux premiers voisins demeure 
inchange sous la transformation tandis que le terme d'interaction change de signe 
et fait apparaitre un terme supplementaire: 

+ ^ E n iT n H -^N^H Q -U^2 Wj. + - (n - U/2) ^(n, T - n ti ) (1.3) 

Ainsi, le modele de Hubbard attractif (U < 0) a remplissage quelconque se 
transforme en modele de Hubbard repulsif a demi-remplissage, n = 1, sous champ 
magnetique, h = [i — U/2. Etant donne que pour le cas demi-rempli le potentiel 
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chimique est donne par fi = U/2, on voit que dans ce cas le champ magnetique 
obtenu par la transformation est mil. 

Le changement du hamiltonien sous la transformation (1.2) pour le 
cas demi-rempli est associe a une symetrie SU(2) du modele. A cette 
symetrie s'ajoute la symetrie de spin (voir [121]). Ce qui donne une symetrie 
SU(2) x SU(2)/Z 2 = SO(4) 3 . Hors demi-remplissage, la symetrie sous rotation de 
spin SU(2) est toujours presente, cependant, l'ajout d'un champ magnetique sous 
la transformation (1.2) brise la symetrie supplement aire qui est reduite au groupe 
U(l). Ainsi la symetrie totale est reduite a SU(2) x U(l). 

La symetrie qui nous interesse est celle associee au parametre d'ordre et 
non pas toute celle presente dans le hamiltonien. Pour le cas attractif, on sait 
qu'a basse temperature il y aura apparition de fluctuations supraconductrices. 
Ainsi le parametre d'ordre sera associe aux operateurs de creation et de 
destruction de paires. Sous la transformation (1.2), ces operateurs se transforment 
de la fagon suivante: 

^^,V^(-i) ,+ 4(sr+^) 

Par ailleurs, on sait que pour le modele de Hubbard repulsif en l'absence de 
champ magnetique, les trois composantes de spin ont des fluctuations critiques. 
Done, en juxtaposant le resultat precedent avec la transformation du hamiltonien 
(1.3), on observe que pour le cas attractif a demi-remplissage, des fluctuations 
de densite de charge s'ajoutent aux fluctuations supraconductrices. Sous ces 
conditions, le parametre d'ordre presente une symetrie SU(2) qui est suffisamment 
importante qu'il soit impossible d'avoir un point critique a temperature finie pour 
ce systeme bidimensionnel 4 . En comparaison, pour le cas hors demi-remplissage, la 

3 La premiere symetrie SU(2) couvre en partie l'invariance sous rotation du spin. 

Ainsi, pour le calcul de la symetrie totale du modele, il faut enlever le facteur Z 2 . 

4 Avec une symetrie SU(2), il n'y a plus de configurations topologiquement dis- 

tinctes en deux dimensions (plus de vortex). 
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transformation du hamiltonien fait apparaitre un champ magnetique qui brise la 
symetrie de spin pour U > et done la symetrie du parametre d'ordre du 
cas U < est reduite a une symetrie U(l). Dans ce dernier cas il existe un 
point critique de type (BKT) a temperature finie pour un systeme en deux 
dimensions. Ceci nous emmene a la conclusion que la temperature critique BKT 
qui augmente avec le remplissage de la bande doit atteindre un maximum puis 
redescendre vers la temperature nulle au fur et a mesure que Ton s'approche du 
demi-remplissage tel que presente a la figure 2. L'abaissement de la temperature 
critique pres du demi-remplissage alors que la temperature ou apparaissent les 
fluctuations est elevee, permet d'avoir un vaste domaine de temperature ou 
existerait un pseudogap. 

Void comment nous prevoyons que pourrait se developper le gap dans le 
poids spectral. Nous nous interesserons a quatre points pres du demi-remplissage 
tels qu'identifies sur la figure 2. A la figure 4 nous avons represents de fagon 
schematique quelle serait, selon notre conjecture, la forme du poids spectral pour 
chacun de ces points. A temperature suffisamment elevee (mais petite par rapport 
a Penergie de Fermi), point a, la theorie des liquides de Fermi peut s'appliquer et 
le poids spectral au niveau de Fermi a la forme d'une lorentzienne dont la 
largeur diminue avec la temperature. En abaissant la temperature sous T* , les 
fluctuations supraconductrices deviennent importantes, mais la symetrie elevee 
du parametre et la bidimensionalite du systeme empechent l'avenement d'une 
transition de phase. C'est alors, qu'apparaitrait un pseudogap dans le poids 
spectral, point b. II est possible qu'avant que ne se developpe le pseudogap, le pic 
present a haute temperature s'elargisse avec l'abaissement de la temperature 
contrairement aux caracteristiques des liquides de Fermi. Un tel phenomene peut 
amener la presence d'un pseudogap a faible frequence dans la densite total 
d'etats sans qu'il n'y en ait dans le poids spectral. II peut aussi y avoir des 
singularites associees a la forme emboitee de la surface de Fermi [63]. 

Une fois que Ton passe sous la temperature critique BKT, le poids spectral 
s'annule au niveau de Fermi, point c. On sait qu'a partir de ce point et jusqu'a 
temperature nulle le parametre d'ordre prend generalement une valeur non-nulle 



15 



Motivation et mise en contexte 



mais avec une phase differente d'une region de l'espace a une autre. II peut 
aussi exister des regions ou il demeure nul (coeur des paires vortex-antivortex). 
Globalement, il ne prendra une valeur non nulle qu'a temperature nulle. Qu'en 
est-il pour le gap (c'est-a-dire, le domaine de frequences ou le poids spectral 
s'annule)? II pourrait prendre une valeur finie pour toutes temperatures sous 
T BKT , mais ce n'est pas stir. Ce que Ton peut affirmer, c'est qu'a temperature 
nulle, point d, le gap est fini. 

Dans la section suivante, nous verrons que les deux premiers spectres 
presentes a la figure 4 ont ete obtenus par certaines etudes du modele de 
Hubbard a demi-remplissage. 

1.4 Quelques etudes anterieures du modele de Hubbard 

En terminant ce chapitre, nous presentons les resultats obtenus par 
differentes methodes qui ont ete appliquees au probleme du pseudogap dans le 
modele de Hubbard attractif. 

Etant donne la difficulte a construire une approximation analytique 
satisfaisante pour decrire le modele de Hubbard, certaines etudes ont utilise 
des techniques Monte Carlo afin d'etudier l'importance relative des differents 
diagrammes obtenus d'un developpement perturbatif [18,22]. Ces etudes etant 
appliquees au modele repulsif, on peut a demi-remplissage les extrapoler au 
modele attractif. Parmi les resultats obtenus, Chen et al. [22] ont montre qu'une 
approximation de type RPA avec un vertex renormalise permet de bien expliquer 
le facteur de stucture magnetique observe par simulations Monte Carlo pour 
une interaction intermediaire U ~ 4. On verra au chapitre 3 que cet accord 
est egalement assez bon pour la susceptibilite de paires du modele attractif. 
L'approximation TPSC a l'avantage de bien predire la valeur de ce vertex pour 
le modele repulsif. 

L 'accord entre la susceptibilite magnetique d'une approximation RPA et les 
resultats Monte Carlo a aussi ete observe par Bulut et al. [18]. lis ont aussi 
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montre que ce genre d'approximation permet de bien determiner la self-energie en 
frequences imaginaires mettant ainsi en importance les corrections de vertex, 
lis ont compare la valeur effective de l'interaction particule-particule avec un 
calcul du second ordre oil la correction de ce vertex provient de l'echange de 
fluctuations de spin. Leur resultat montre que l'accord n'est pas tres bon et n'est 
pas ameliore de fagon significative par l'ajout du troisieme ordre. 

Malgre ces resultats peu encourageants, plusieurs auteurs ont utilise 
1' approximation FLEX pour l'etude du modele de Hubbard. II faut mentionner 
que cette approximation va au-dela de l'unique contribution de l'echange de 
fluctuations de spin, elle tient compte egalement des fluctuations de charge et des 
fluctuations de paires de Cooper. Recemment, Deisz, Hess et Serene [32] ont 
applique cette approximation au probleme du pseudogap dans le modele de 
Hubbard attractif. lis ont observe que cette approximation predit des proprietes 
anormales dans le phase normale, tel la suppression du poids spectral a faibles 
energies avec l'abaissement de la temperature. Cependant, aucune evidence pour 
l'apparition d'un pseudogap n'a ete observee et ce meme pour des systemes de 
taille 128 x 128. 

La technique Monte Carlo permet d'aller plus loin qu'une simple 
comparaison avec les approches analytiques en frequences imaginaires. Utilisant 
des methodes pour faire un prolongement analytique des donnees Monte Carlo, 
on peut ainsi verifier si ces donnees predisent ou non la presence d'un pseudogap 
dans le poids spectral des excitations fermioniques. Un resultat reproduit par 
toutes les simulations Monte Carlo est la presence d'un pseudogap pour les 
systemes de petites tailles. Le comportement en taille observe par Vekic et White 
[118] tend a montrer que ce pseudogap disparait lorsque le systeme est de taille 
suffisamment grande. Plus recemment, nous avons publie une etude du pseudogap 
dans le modele de Hubbard repulsif a demi-remplissage [80]. La symetrie entre ces 
modeles montre que ce resultat serait retrouve pour le modele attractif. Nous 
avons observe la presence d'un pseudogap pour des systemes de differentes tailles 
et notre etude en taille a montre que celui-ci n'est du a des effets de taille finie. 
Nous avons montre que l'approche TPSC developpee par Vilk et Tremblay 



17 



Motivation et mise en contexte 



[116,117] montre un accord parfait avec les resultats Monte Carlo pour les 
systemes de petites tailles. Cette observation suggere d'utiliser cette methode afin 
de determiner s'il existe ou non un pseudogap dans la limite thermodynamique. 
Notre conclusion est qu'il en existe un. Dans cet article, nous expliquons pourquoi 
le resultat de Vekic et White differe du notre. En effet, nous avons montre que 
plus le bruit est important dans les donnees Monte Carlo, moins on observe 
de details dans le poids spectral obtenu par l'application de la technique de 
maximisation d'entropie aux donnees Monte Carlo (cette technique est expliquee 
en annexe D). Ainsi, l'observation que les resultats de Vekic et White sont plus 
precis pour les systemes de petites tailles que pour les systemes de grande taille 
renforce notre argumentation selon lequel il y aurait bien un regime pseudogap 
dans la limite thermodynamique. 

D'autres approches analytiques prevoient egalement l'existence d'un 
pseudogap. Parmi celles-ci, il y a une approche par transformation de 
Hubbard-Stratonovich [19,20] ou Ton observe l'apparition d'un pseudogap a une 
temperature legerement superieure a la temperature critique. II y a aussi 
l'approximation de type matrice-T [56,64]. On remarque pour cette derniere 
approximation que lorsque le calcul est entierement auto-coherent le pseudogap 
n'est plus observe. Tout recemment, on a vu le developpement de l'approximation 
d'amas dynamique (DCA). Celle-ci a egalement ete appliquee au probleme du 
pseudogap [50]. II semble qu'en ameliorant la precision de cette approximation, la 
presence d'un pseudogap dans le spectre devient de plus en plus marquee. Ainsi, 
cette methode renforce notre resultat. 
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FORMALISME 

Dans ce chapitre, nous explicitons les differentes etapes menant a 
1' approximation auto-coherente a deux particules. Pour notre derivation nous 
adoptons un formalisme developpe dans les annees 50 et 60. Nous referons le 
lecteur aux articles de Martin et Schwinger [76], de Kadanoff et Martin [55], 
de Baym et Kadanoff [10,11] et de Luttinger [70,71,72,60]. Pour des ouvrages 
pedagogiques nous suggerons le livre de Mahan [73] ainsi que les notes de cours 
de Tremblay [110]. 

Pour resumer les principales particularites de notre formalisme, mentionnons 
que nous utilisons Papproche en temps imaginaire avec integration sur le contour 
de Matsubara. De plus, nous employons la notation compacte [76] ou les chiffres 
arabes sont utilises pour representer un point dans l'espace a trois dimensions 
forme par le reseau cristallin et le temps imaginaire. Etant donne que Ton 
s'interesse au modele de Hubbard attractif, nous faisons usage de la notation 
matricielle de Nambu. 

Le chapitre est separe en quatre sections. Dans la premiere section nous 
obtenons des resultats exacts et definissons les differentes quantites utilisees par 
la suite. Nous considerons la fonction de Green qui permet de decrire l'etat 
thermodynamique d'un systeme. En outre, en presence d'un champ externe, cette 
fonction permet de deduire les proprietes de transport du systeme. On peut alors 
affirmer qu'elle contient toute l'information necessaire pour decrire le systeme. La 
premiere relation exacte que nous etudierons est l'equation du mouvement de la 
fonction de Green. Pour le modele de Hubbard, cette equation nous emmene a 
considerer le hamitonien en deux parties. La premiere partie est diagonalisable de 
fagon analytique et assez facilement. La deuxieme introduit la self-energie qui est 
liee a la fonction de Green par la relation de Dyson. En fait, cette derniere 
equation peut etre consideree comme une definition de la self-energie et est done 
applicable a tout hamiltonien. 
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L'equation du mouvement permet de faire un lien entre les proprietes dites 
a une particule (fonction de Green, self-energie) et la fonction de correlation a 
deux particules. C'est pourquoi nous presentons dans la deuxieme section une 
relation exacte permettant d'evaluer cette fonction de correlation en l'absence de 
champ externe et dans la phase desordonnee. Cette relation est l'equation de 
Bethe-Salpeter. Celle-ci fait appel a un nouvel objet qui est le vertex irreductible. 
II sera done essentiel d'approximer ce vertex pour evaluer la fonction de 
correlation. Le lecteur pourra noter que pour la derivation de cette formule nous 
ne faisons appel qu'a la relation de Dyson et pas a l'equation du mouvement. 
Ainsi l'equation de Bethe-Salpeter pourrait s'appliquer a d'autres hamiltoniens. 

La troisieme section contient le developpement de l'approximation. Pour ce 
faire nous suivons un formalisme qui s'apparente a celui de Baym et Kadanoff 
[10,11] en ce sens que nous developpons une approximation pour la self-energie en 
presence du champ externe. De cette meme approximation nous estimons la 
fonction de Green et le vertex irreductible. Etant donne que Ton s'interesse aux 
proprietes du modele en l'absence de champ externe, notre approximation n'est 
developpee qu'au premier ordre afin de determiner le vertex. Notre approximation 
pour la self-energie fait appel a une quantite inconnue: la valeur moyenne de 
double occupation (deux electrons sur un meme site du reseau). Suivant la 
proposition de Vilk et Tremblay [116,117], nous determinons cette valeur de fagon 
auto-coherente en imposant une regie de somme. Par la suite, utilisant la relation 
de Bethe-Salpeter et l'equation du mouvement, nous pouvons raffiner l'estimation 
de la self-energie. 

Enfin, a la derniere section, nous montrons que notre approximation permet 
de satisfaire une nouvelle regie de somme, qui est Panalogue de la regie de somme 
f, et le theoreme de Mermin- Wagner, selon lequel il est impossible d'avoir une 
brisure spontanee de symetrie continue a temperature finie en deux dimensions. 
Le fait que notre approximation respecte les moments de la susceptibilite de 
paires nous permet de faire un parallele entre l'approche des moments et la regie 
de somme choisie pour determiner la valeur moyenne de double occupation. 
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2.1 Relations exactes 
2.1.1 Definitions 

Nous cherchons une approximation pour decrire le modele de Hubbard 
attractif. Nous savons qu'a basse temperature ce modele presente des fluctuations 
supraconductrices importantes. II est alors approprie d'adopter la notation de 
Nambu oil on definit des operateurs vectoriels: 

*(1)=(^))) ^(1) = (^{(1), ^(1)) (2.1) 

Les fleches f et [ servent a indiquer les deux orientations possibles du spin 
electronique. La dependance temporelle de ces operateurs est definie de fagon 
conventionnelle par l'operateur devolution en temps imaginaire: 

= ^(r^Ti) = e ^( w -^)^( ri )e- ri ^-^ 
^t(l) = !pt(r 1>Tl ) = e n(W-^)^t( ri ) e -n(W-AiiV) ^ 

oil N est l'operateur du nombre de particules, \i, le potentiel chimique et 
les variables de temps imaginaire sont definies dans l'intervalle (0,(3), (3 etant 
l'inverse de la temperature. Suivant les definitions (2.2), l'operateur de creation 
W\t) ne sera pas le conjugue hermitique de l'operateur d' annihilation &(t). 
Pris au temps r = 0, il le sera et ces deux operateurs obeiront a la regie 
d'anticommutation: 

{^(r 1 ),^(r 2 )} = 5 i)i 5(r 1 -r 2 ) 



Dans la notation de Nambu la fonction de Green prend une forme 
matricielle. 

G(l,2;0 = -<tf(l)tft( 2 )), 

_/ (^(1)^(2)), (^(l)Vi(2)>^ (2-3) 

WliM®)* (^I(i)^(2)) J 

On reconnait le long de la diagonale la fonction usuelle pour les deux 

composantes de spin. Hors diagonale, on retrouve les fonctions de Green 
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anormales presentees par Abrikosov et al. [1]. La moyenne thermodynamique a ete 
definie en presence d'un champ source externe: 

<o>, = Z- 1 [0]Tr{e-^-^T T exp(-^t(T)0(T,2)^(2))o} (2.4) 

ou 

Z[<j>] = Tr (e'^-^T^p (-!pt (1)0(1, 2)lP(2))) (2.5) 

Pour indiquer que Ton integre sur un point, on place une ligne au-dessus du 
chiffre associe a ce point (ex.: 1). L'operateur d'ordonnance temporelle T T est 
utilise pour indiquer que les operateurs doivent etre ordonnes selon leur indice de 
temps imaginaire. II faut d'abord faire le produit direct puis ordonner, dans 
chaque element de matrice, les operateurs en suivant le contour de Matsubara, 
c'est-a-dire en les plagant de la droite vers la gauche en ordre croissant de temps 
imaginaire. 

Le champ externe est choisi afin de pouvoir evaluer les correlations 
supraconductrices, c'est-a-dire: 



Les indices inferieurs 12 et 21 servent pour la notation des deux elements (1,2) 
et (2, 1) de la matrice 0. Ainsi, nous avons introduit deux champs distincts qui 
seront couples aux deux termes hors diagonaux de la fonction de Green. On peut 
verifier que la derivee de la fonction de partition, (2.5), permet de generer ces 
termes hors diagonaux. 

flnZ[0] / G 12 (l,2;4 



ou on a defini: 



50(2,1) V G 2i(l,2;0) 
S ( 



<ty(2,l) 




(2.7) 



(2.8) 



\ 5012(2,1) 

Les indices 21 et 12 servent a indiquer clairement par rapport a quelle 
composante du champ la differentiation est faite. Cette definition de la matrice 
de derivation nous permet d'ecrire: 



0(1, 2) = 5(1 -4)*(3-2)(j J) 



50(4,3) 
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ou pour appliquer la differentiation nous devons faire le produit matriciel et 
deriver chaque composante a tour de role. 

Avant de terminer cette section, il faut remarquer qu'il reste un point a 
determiner pour completer la definition de la fonction de Green. Le lecteur 
aura note, par la presence du terme fiN, que nous utilisons l'ensemble grand 
canonique ou le potentiel chimique est determine afin d'obtenir le remplissage 
desire. Cependant, que devient cette condition en presence d'un champ externe? 
Nous choisissons de traiter le potentiel chimique de la meme fagon que la 
temperature (les deux etant des parametres de Lagrange permettant d'ajuster le 
nombre moyen de particules et l'energie moyenne respectivement). Ainsi, il n'est 
pas modifie par la presence du champ externe et il est ajuste pour que le 
nombre d'occupation moyen en l'absence de ce champ ait la valeur desiree. 

2.1.2 Equation du mouvement 

Afin d'evaluer la fonction de Green d'un modele donne, il est courant 
de considerer l'equation du mouvement de cette fonction. On precede done en 
differential G((f>) par rapport au temps. Trois termes sont generes par cette 
derivation. Les deux premiers proviennent de l'ordonnance temporelle: il y a 
l'ordonnance d'un vecteur de Nambu par rapport a l'autre (voir equation (2.3)) 
et il y a l'ordonnance temporelle par rapport au champ externe. Le troisieme 
terme vient de la dependance temporelle de Foperateur (2.2). On peut ecrire la 
derivation ainsi: 



La notation utilisee dans l'equation (2.9) est un peu ambigue. II faut bien 
comprendre ce que Ton veut representer par l'anticommutation entre bra et ket. II 
faut dans un premier temps faire le produit direct des deux elements dans l'ordre 
ou ils apparaissent dans l'anticommutateur, puis appliquer l'anticommutation 




+ <[tft (3)0(5, 4)^(4), ^(l)]^t (2) ) 5 ( Tl _ r _) 



(2.9) 
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element par element. Cette ambiguite n'apparait que dans le premier terme du 
cote droit de 1' equation (2.9). 



Le troisieme terme peut etre evalue en utilisant l'equation de Heisenberg: 



= [H-fiN,¥(l)] 



En appliquant (2.9) au modele de Hubbard (1.1), on obtient: 
0G(1,2;, 



=5(1-2) 



Mr { - t-) [ ' (r '- ( ^ " 



0(1,3) 



G(3,2) 



U 



(V T (l)^(l)^(l+)^(2)) (^ T (l)^(l)^(l + )^(2))^ 

(v;(i)v}(i^ T (i-)^(2))^ (^(1)^(1)^(1-)^ (2)>^ 



(2.10) 

Etant donne que differents operateurs apparaissant dans une moyenne sont 
evalues a une meme valeur de temps imaginaire, il devient necessaire d'indiquer 
explicitement l'ordre selon lequel il faut ordonner ces operateurs. Pour ce faire 
nous plagons un indice superieur + ou _ qui indique que Poperateur doit etre 
place plus a gauche ou plus a droite lors de l'arrangement en ordre. 



Le dernier terme de la relation (2.10) provient du terme d'interaction dans 
le modele de Hubbard. En l'absence de ce terme, il serait facile de determiner la 
fonction de Green. C'est pourquoi nous allons traiter separement ce dernier terme 
en introduisant la self-energie: 

^(1,2;0) = f {^ T (1)^(1)W( 1+ )^T(3)>0 (V'tW|(1)^(1 + )^(3))^ g _ 1( 3 52;0) 



<^(1)V4(1)V t (1-)^(3))^ (^(1)^(1)^(1-)^ (3)> . 

(2.11) 

On voit par cette relation que le calcul de la self-energie fera appel aux fonctions 
de correlation. Cette definition, avec (2.10), nous donne la relation de Dyson: 



(Gq-^I, 3) - 0(1, 3) - 17(1, 3; 0))G(3, 2; 0) = 6(1 - 2) 



(2.12) 
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ou on a defini G par les autres termes presents dans l'equation du mouvement: 
G-(l, 2 ) = -,(1 - 2)|_ - ^ - . 2 ) ("^.^ t{r J i)+ ^ 

Ce resultat justifie notre choix selon lequel le potentiel chimique est fixe 
par le remplissage en l'absence du champ externe. Ainsi, G est independant de 
ce champ externe. 

Jusqu'a maintenant, nous nous sommes employe a regarder ce que Ton 
obtenait en etudiant la differentiation de G par rapport a un de ses arguments. 
Si Ton avait considere la differentiation par rapport a l'autre argument, la 
relation de Dyson aurait pris la forme suivante: 

G(l, 3; 0)(Gq 1 (3, 2) - 0(3, 2) - 17(3, 2; 0)) = 6(1 - 2) (2.13) 

Etant donne que la relation de Dyson nous servira a estimer la fonction de 
Green a partir d'une approximation de la self-energie, on peut se demander si le 
resultat serait le meme si au lieu de considerer la relation (2.12), nous utilisions 
(2.13). La reponse a cette question nous vient simplement de l'algebre lineaire: 
tant que la fonction de Green ne sera pas singuliere son inverse sera unique. 

2.2 Equation de Bethe-Salpeter 
2.2.1 Definition 

On a vu de l'equation du mouvement que revaluation de la self-energie 
passe par le calcul d'une fonction de correlation a quatre points. Or, il est 
possible de montrer que cette fonction de correlation est en fait la reponse de la 
fonction de Green au champ externe. Ainsi, on definit 

x (l,2,3,4;0) = -^- T yG(3,4;0) 

= / (^(1)^(2)^(3)^(4))^ (^(1)^(2)^(3)^(4))^ (2.14) 
1^1(1)^(2)^(3)^(4))^ (V|(1)V}(2)^(3)^(4))J 
-G HD (1,2;0)G(3,4;0) 
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ou l'indice HD indique que seule la partie hors diagonale de la matrice est 
conservee, les elements de la partie diagonale etant remplaces par zero. Nous 
utiliserons aussi l'indice D qui ne conserve que la partie diagonale de la matrice. 
Le dernier terme de l'equation (2.14) vient de la presence de Z~ x \<$\ dans 
la definition de la moyenne thermodynamique (2.4). La matrice de derivation 
utilisee est la meme que celle introduite precedemment (2.8). La fagon de 
l'appliquer consiste a faire la multiplication matricielle dans un premier temps 
puis d'appliquer la derivation appropriee a chaque element. Nous illustrons par 
l'exemple de la derivation d'une matrice A(<p) quelconque: 



50(1,2) 



A{4>) 



/ 6A 21 (<t>) SA 22 (<t>) 
.._ [ <5</> 2 i(l,2) <S02i(l,2) 
- Miiffl 8A 12 (<t>) 



i(l,2) 5</, 12 (l,2) 



La comparaison avec (2.11), nous montre que la self-energie sera alors 
donnee par: 



17(1, 2; cj>) = U X (1, 1,1^3; ^G"^, 2; <j>) + 176(1* - 2)G HD (1, 1; , 



(2.15) 



Nous avons utilise comme indice superieur ± pour indiquer que dans la ligne 
superieure de % le premier operateur sera considere a un temps infinitesimalement 
superieur a r x tandis que dans la deuxieme ligne de la matrice, il sera considere 
comme etant a un temps inferieur. 

La definition de la fonction de correlation (2.14), comme etant la reponse 

de la fonction de Green au champ externe, nous indique comment nous pourrions 
l'estimer. De la definition: 

G(l,3;</>)G- 1 (3,2; ( />) = ( 5(l-2) (2.16) 

nous voyons que la dependance de G sur <j> devrait compenser exactement celle 
de G _1 . Au premier ordre en </>, la fonction de Green est donnee par: 

G(l, 2; <P) = G(l, 2) - 0(4, 3) X (3, 4,1,2) + O(0 2 ) (2.17) 
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II est a noter que le terme correspondant a la transposee de n'apparait pas 
dans cette expression car nous sommes interesses par la phase normale oil 



50(4,3) 



G(l,2) 



= 



Ainsi, la relation de Bethe-Salpeter que nous derivons ne sera valide que dans 

l'etat desordonne ou les termes hors diagonaux des matrices G et % sont nuls. 

Pour G -1 , on peut utiliser Pequation de Dyson (2.12), pour montrer qu'au 
premier ordre on a: 



G-^l, 2; 0) = G-^l, 2) - 0(1, 2) - 0(4, 3) 



-= _ 11(1,2;. 
L 50(4, 3) V 



+ O(0 2 ) (2.18) 



Encore une fois, il devrait apparaitre un terme supplemental. Cependant, on 
peut montrer que ce terme s'annule au premier ordre en 0: 



(50(4,3) 



E(l,2) 



= 



!>=0 



ou T indique la transposee de la matrice. 



Introduisant (2.17) et (2.18) dans (2.16), on trouve une relation pour \ en 
annulant le coefficient du terme proportionnel a (f>: 



X (l,2,3,4) =^6(3,2)^6(1,4) 

~d 



a x G(3, 5)<T a 



50(2, 1) 



17(5,6) 



G(6,4) 



(2.19) 



b=0 



ou a x est la matrice de Pauli correspondant a la composante x du spin. Cette 
matrice provient de la commutation de la matrice diagonale G et de la matrice 
hors diagonale 0: 

G(1,2)0(4,3) = 0(4,3KG(1,2K 

La matrice G est diagonale par le fait qu'on l'evalue en champ externe nul et a 
temperature finie (pas d'ordre a longue portee). 
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2.2.2 Transformation de Legendre 

On voit de la relation (2.19) que la determination de la fonction de 
correlation passe par l'estimation de la variation de la self-energie par rapport au 
champ externe. C'est ce terme qui permet d'introduire le vertex irreductible. Mais 
pour y arriver, nous allons considerer une transformation de Legendre. L'utilite de 
cette transformation s'explique par le fait que la dependance de la self-energie sur 
(j) se fait au travers de la fonction de Green £{4>) = lg=G H D(</>)" Ainsi, de 
fagon generale, on peut definir la transformation de Legendre de la fonction de 
partition: 

W[G]= (-lnZ[0]-Tr(G HD (T,2;0)0(2,T))) (2.20) 

V / Gun \<p\=y 



avec la condition 

SlnZ[(/)} 



50(2,1) 



Ghd [4>]=Q 

= -0(1,2) (2.21) 



GhdO)=G 

L'indice apres la derniere parenthese du cote droit de l'egalite signifie que la 
fonction a l'interieur des parentheses est evaluee a une valeur de <fi telle que les 
composantes hors diagonales de G(0) soient egales au nouveau champ introduit 
Q. Cette condition ne s' applique que sur les composantes hors diagonales car 
comme le champ <fi ne possedait que deux composantes non nulles, il doit en 
etre de meme du champ Q, celui-ci etant introduit par une transformation de 
Legendre (Les details de cette transformation sont donnes en annexe B). Pour 
montrer que toute la dependance de U sur <fi peut s'exprimer a travers Q (voir 
annexe B), il suffit d'examiner la derivation de W par rapport a Q: 



dU =-*(l,2;fi?)- ( 5 —^m4-,g))( 5 —^\nZ[ ( f>} + g(4,3) 



5g(2,l) v ' ,z " \8Q{2,\) v ' '"7V 4) 

= -*(l,2;a) 

oil la matrice de derivation a ete definie par: 

8 ( 



=*(£?) 

(2.22) 



5021(2,1) 



W(2,l) VfelTO) J 

et &{G) est la valeur du champ <j> obtenue en posant la condition (2.21). Cette 
derniere presuppose que pour une valeur de <fi donnee une seule valeur possible 
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de G HD (</>) soit donnee par la derivation de la fonction de partition (2.7) et que 
pour une valeur de G HD (</>) = Q donnee une seule et unique valeur de <fi 
permettre d'obtenir cette fonction de Green. 



Utilisant la transformation de Legendre introduite, on est a meme d'evaluer 
la variation de la self-energie. Nous appliquons la regie de derivation en chaine 
qui est demontree en annexe A 



50(4,3) 



27(l,2;</>) 



■G HD (5,6;0) 



(p=0 



==-27(1, 2; Q) 

50(5,6) 1 ; 



^=o L 50(4, 3) 

II semble alors opportun d'introduire le vertex irreductible qui est defini par: 

S 



r(l,2,3,4;£) 



^(2,1) 

Ceci nous emmene a la relation de Bethe-Salpeter 



27(3, 4; Q) 



(2.23) 



X (l, 2, 3, 4) = -(7,6(3, 2KG(1, 4) 

+ a x G(3, 5)a xX (h 2, 7, 8)r(8, 7, 5, 6)G(6, 4) 



(2.24) 



2.3 Estimation des vertex irreductibles 

2.3.1 Approximation de la self-energie 

Ann d'estimer les vertex irreductibles nous allons utiliser une approche qui 
s'apparente au formalisme de Baym [11] et faire dans un premier temps une 
approximation de la self-energie en presence d'un champ externe. Puis le vertex 
sera estime en utilisant la definition (2.23). II faudra done exprimer 27 en 
fonction de Q. Etant donne que Ton s'interesse uniquement a la phase normale ou 
les termes hors diagonaux des matrices sont nuls, on peut se contenter d'evaluer 
27 au premier ordre en Q. 

Comme approximation, Vilk et Tremblay [117] ont propose dans l'esprit de 
l'approche de Singwi et al. [105], de considerer la self-energie au premier ordre en 
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developpement perturbatif et de corriger cette approximation par un facteur 
multiplicatif qui soit tel que l'approximation puisse, en principe, etre exacte 
lorsque Ton considere les correlations locales. Pour bien indiquer que ce sera la 
notre premier niveau d'approximation nous placerons un indice superieur ^ a 
toutes les quantites approximees. L'approximation est la suivante: 

2^(1, 3; G)G(3, 2; Q) « UA(1, £)G(1, 2; Q) (2.25) 

La matrice A sera determinee par la valeur des fonctions de correlations locales. 
Nous pouvons par exemple poser: 

2<-V(l, 3; £)G«(3, 1+; Q) = U X {1) (1, 1, 1 ± , 1 ++ ; G) + UG(1, 1)G«(1, 1+) (2.26) 

L'indice superieur ++ sert a indiquer l'ordonnance de l'operateur correspondant 
par rapport aux operateurs evalues au meme temps t x et egalement a t±. Nous 
avons clairement indique qu'il faut exprimer % et G comme des fonctions de Q 
afin d'obtenir A comme une fonction de Q. Encore une fois on peut se contenter 
du premier ordre. Pour la fonction de Green, c'est assez simple: 

G«(1,2;£) = G«(1,2) + 0(1,2) 

Pour la fonction de correlation, on utilise le principe de Pauli considere 
localement qui, par le fait que Ton ne s'interesse qu'aux termes locaux, permet 
de simplifier l'expression: 

xc'd, i, i* s) = ( { "< n ^ I^a-^X 5 ) - e(1 - 1)G(1)(1 - 1+) 

= ( < T' ) Kni ° ) _ 3 ) +6l2 (i,iK-e(i,i)G«(i,i +) 

Ainsi, la relation (2.26) nous donne une relation permettant d'approximer 
A apparaissant en (2.25) et done la self-energie S. Cependant, on peut se 
demander si l'approximation serait la meme si Ton imposait la relation suivante: 

17^(1, 3; £)G«(3, 1-;G) = U X {1) (1, 1, 1 ± , 1~ ; G) + UQ(1, 1)G«(1, 1") (2.27) 
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On peut montrer que ces deux approximations donnent le meme resultat pour 
ce qui est de la partie diagonale du vertex irreductible obtenu par derivee 
fonctionnelle. Cependant, meme apres avoir eteint le champ externe, la partie 
diagonale de la self-energie demeure differente selon l'approximation choisie. Ann 
de determiner l'approximation qui serait la plus appropriee, considerons ces 
relations exprimees comme des sommes des composantes de Fourier. On obtient 
pour la partie gauche de l'egalite de (2.26): 

lim —S"s( 1 \k;g)G {1) (k;G)e~ iknri 

k 

Le facteur de convergence est essentiel car aux hautes frequences la self-energie 
est constante, lim fc ^ S(k, ik n ) ~ tandis que la fonction de Green decroit en 
loi de puissance liim.^^ G(k, ik n ) ~ Cette lente convergence de la somme 
nous indique que les composantes a hautes frequences contribuent de fagon 
importante a la somme. De meme, lorsqu'on considere la transformer de Fourier 
de (2.27), on obtient un resultat identique a l'exception du signe dans 
l'exponentiel du facteur de convergence. Cette difference de signe nous indique 
que si Ton considerait la somme de (2.26) et de (2.27), le facteur de convergence 
ne serait plus necessaire. On le verifie par le fait que: 

lim (U(k, ik n , G)G(k, ik n , Q) + 27(k, -ik n , G)G(k, -ik n , G)) oc — 

Ainsi, pour la somme sur les frequences de Matsubara, une contribution plus 
importante vient des basses frequences lorsqu'on considere la somme de (2.26) et 
de (2.27). C'est pourquoi nous choisissons A tel que: 

A(i, G) =0(1, i) + \ (x (1) (i, i, i*, £)(G (1) (i, i + ; G)Y' 

x (2-28) 

+ x (1) (i,i,i ± ,i--;a)(G( 1 )(i > i- ; a))- 1 ) 

De ce resultat, on est a meme de deduire la self-energie et sa dependance 
au premier ordre en Q. Ainsi la definition du vertex irreductible (2.23) nous 
permet de creer un lien entre ce vertex et les correlations locales. Introduisant 
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l'equation (2.28) dans 1' approximation pour la self-energie (2.25), on obtient, 
apres avoir pose = 



^ (1) (1,2) =£/<5(l - 2) 



(n T )(l - (n T )) 
= i(c/-(l-n)C/ pp )<5(l-2) ( 7 z 



et en utilisant la definition du vertex irreductible (2.23): 



r (1) (l,2,3,4) =U5(1 - 


2)5(1- 


3)6(1- 


4) T w 1 1 
(n i )(l -n T ) 


=U pp 6(l 


-2)5(1 


-3)5(1 


-4)1 



(2.29) 



(2.30) 



ou <j z est la matrice de Pauli correspondant a la composante z du spin. On a 
introduit la variable U pp qui correspond a la valeur renormalisee du vertex. 



Malgre que nous ayons demontre une fagon de deriver une approximation 
pour la self-energie et le vertex, ces relations devraient etre considerees comme 
des ansatz qui seront utilises dans le calcul des fonctions de Green a une et a 
deux particules. 



2.3.2 Determination de la fonction de correlation 

Comme on s'interesse a la phase normale sans champ externe, toutes 
les matrices sont diagonales et nous pouvons evaluer les composantes (1,1) 
independamment des composantes (2,2). C'est pourquoi a partir de cette section 
et pour le reste du chapitre nous nous concentrerons sur les composantes (1, 1) de 
chaque matrice. Nous ne l'indiquons pas de fagon explicite afin de ne pas alourdir 
la notation. Toutefois, le lecteur attentif, remarquera que nous n'utiliserons plus 
les caracteres gras. 

Le calcul de E et de F nous permet d'estimer la fonction de correlation. 
Mais avant de pouvoir le faire nous devons determiner la fonction de Green en 
utilisant la relation de Dyson, (2.12). Pour la determiner, nous procederons de la 
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meme fagon que dans le formalisme de Baym et Kadanoff [10,11] en utilisant 
la meme self-energie que celle que nous avons utilisee pour deriver le vertex 
irreductible. On peut ainsi dire que Ton considere la fonction de Green et le 
vertex irreductible a un meme niveau d'approximation. 



On peut remarquer que suivant notre premiere approximation, la 
self-energie n'est en fait qu'une correction du potentiel chimique. En appliquant 
une transformation de Fourier on peut diagonaliser les termes sans interactions. 
Ce faisant, on obtient la fonction de Green: 



(2.31) 



oil k = (k, ik n ) et k n = (2n + l)nT est une frequence de Matsubara fermionique 
avec n£ (—00, +00). Nous avons place un indice superieur W au potentiel 
chimique pour indiquer que celui-ci sera ajuste afin que la fonction de Green 
ait le bon remplissage, c'est-a-dire tel que 1 : 




lim G (1) (k)e- ik ^ = 

— 0- N ^ 



n 



(2.32) 



II est interessant de verifier le cas demi-rempli oil la valeur du potentiel 
chimique est connue exactement, fx = U/2. On sait que pour que le nombre moyen 
de particules correspondant a soit n = 1, il faut avoir /i^ 1 -* — S^ 1 -* = 0. Or, a 

demi-remplissage la self-energie obtenue, eq.(2.29), est donnee par T,^ = U/2. 
Done, le potentiel chimique prend sa valeur exacte /i^ = U/2. 

Maintenant nous pouvons determiner la fonction de correlation. Etant 
donne que dans notre approximation le vertex irreductible est local et statique, 
eq.(2.30), nous pouvons nous restreindre a revaluation des fonctions de correlation 

1 Cette condition accompagnee du fait que la self-energie, S^ 1 ) n'est qu'un 
deplacement du potentiel chimique nous montre que la valeur exacte de cette derniere 
quantite n'aura pas d'importance dans l'estimation de 
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a deux points: 1, 2, 2) = x(r x — r 2 , r x — r 2 ). Cette fonction correspond a la 

fonction de correlation pour les paires locales. En utilisant la relation de 
Bethe-Salpeter, (2.24), on peut montrer qu'apres avoir applique la transformation 
de Fourier, on obtient: 



(2.33) 




ou 



X {1) (q) = J d 2 r j\rx(T,T)e^ 



r-tq-r 



avec q n = 2nirT une frequence de Matsubara bosonique. La fonction de correlation 
obtenue est de la meme forme que celle obtenue d'une approximation de type 
matrice-T a l'exception du fait que le vertex est renormalise. De plus, le 
numerateur est donne par 

T 
N 



xo 1) (9) = 47E G(1) (-^+^) G(1) ^) 



(2.34) 



Comme on peut le remarquer, x^ a 1 & nieme forme que la fonction de 
correlation en l'absence d'interaction. Cependant, il est a noter que la fonction de 
Green qui y apparait est la fonction de Green corrigee par la premiere estimation 
de la self-energie. 



2.3.3 Determination de la double occupation 

Le developpement presente precedemment nous a permis d'ecrire le vertex 
irreductible comme une fonction de la valeur moyenne d'occupation double (2.30). 
Cependant, cette quantite est inconnue. Differentes approches permettent de regler 
ce probleme: on peut, par exemple, considerer la valeur obtenue par simulation 
Monte Carlo, ou, comme nous le verrons dans le chapitre suivant, celle obtenue 
d'un calcul champ moyen a temperature nulle. La solution proposee ici consiste 
a considerer une regie de somme arm de fermer le systeme d'equations sans 
appel a d'autre information exterieure. Ce choix nous assure ainsi que notre 
approximation respectera cette regie de somme. 
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La regie de somme que nous allons considerer donne une autre relation 
liant le vertex irreductible a la valeur moyenne de double occupation. Cette regie 
de somme decoule du theoreme de fluctuation-dissipation. 

A ce niveau, on substitue des expressions approximatives qui permettent d'obtenir 
(n^n^) de fagon auto-coherente. 

Etant donne le comportement de la susceptibilite a haute frequence, 
lim _ +oo x (zq' n ) oc l/(iq n ), le facteur e~ iqnV est essentiel pour assurer la 
convergence de la somme. Cette regie de somme est en fait la serie de 
Fourier-Matsubara et ce facteur de convergence sert a indiquer l'ordonnance 
temporelle des operateurs. 

Ainsi, les relations (2.30), (2.31), (2.32), (2.34) et (2.35) forment un 
systeme d'equations qui nous permet de determiner, de fagon approximative, les 
proprietes a deux particules du modele de Hubbard. Parfois, on a caracterise 
cette approche par "imposition d'une regie de somme exacte". Cependant, il est a 
noter qu'il serait plus precis de parler de regie de somme auto-coherente car 
comme nous l'avons mentionne precedemment la double occupation n'est pas 
connue de fagon exacte. Les valeurs obtenues ainsi pour la double occupation sont 
presentees au chapitre suivant ou Ton compare avec les resultats obtenus par 
simulation Monte Carlo. 

2.3.4 RafRnement de l'approximation 

Connaissant les proprietes a deux particules il est possible de rafliner 
notre approximation pour les proprietes a une particule. Utilisant l'equation du 
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mouvement, (2.11), et l'equation de Bethe-Salpeter, (2.24), on peut demontrer la 
relation exacte: 



On voit qu'etant donne qu'a notre premier niveau d'approximation le vertex 
irreductible est local, la seule fonction de correlation contribuant a la self-energie 
sera la fonction de correlation de paires a deux points que nous avons evaluee a 
la section (2.3.2). Utilisant la transformation de Fourier, on obtient ce que Ton 
appellera le deuxieme niveau d'approximation pour la self-energie: 



La relation (2.37) est utile car elle permet de separer la self-energie en deux 
termes. Le premier terme est constant et est done present aux grandes frequences 
comme aux basses frequences. Le deuxieme terme contient les contributions basses 
frequences que nous evaluons au premier niveau d'approximation. 

On peut remarquer que ce resultat differe du resultat obtenu par Vilk et 
Tremblay [116] par le fait qu'il n'apparait dans l'expression de la self-energie 
qu'une seule fonction de correlation et un seul vertex renormalise. On devrait 
s'attendre qu'a demi-remplissage on retrouve leur resultat (a cause de la symetrie 
particule-trou pour ce remplissage) . Cette difference vient du fait que notre 
calcul est fait dans le canal transverse et leur calcul a ete fait dans le canal 
longitudinal. II serait possible de faire le calcul du cas attractif dans le canal 
longitudinal, on retrouverait alors la formule de Vilk et Tremblay pour la 
self-energie. Pour le cas attractif, e'est la susceptibilite de charge qui serait 
critique a demi-remplissage (voir section 1.3). Ainsi la symetrie entre les 
correlations de charge et cedes de paires qui semble etre brisee dans notre 
approximation pourrait etre retrouvee en traitant le canal longitudinal et le canal 
transverse de la meme facon [80]. 



17(1, 2) = -U5(l - 2KG(1^, IK + Ua x G(l, 5)a xX (h 1, 7, 8)r(8, 7, 5, 2) 



(2.36) 




(2.37) 
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La fonction de Green se determine simplement par l'utilisation de l'equation 
de Dyson 



G {2) {k) = - 



(2.38) 



- e(k) + /i( 2 ) -E( 2 )(fe) 

En accord avec l'argument de Luttinger [71], le potentiel chimique est evalue au 
meme niveau d'approximation que la self-energie. Nous devons done egalement 
raffiner l'estimation de ce parametre. 



lim S^G{k)e~ lknn = n 



ik n r) 



N ^ 



lim — > — 

„-»o- N Y ik n - e(k) + - S( 2 )(fc) 



= n 



Avec le formalisme presente jusqu'a maintenant, il semble difficile 
de poursuivre les iterations et d'obtenir une meilleure estimation du vertex 
irreductible. Pour y arriver, il faudrait ameliorer notre estimation des termes 
hors-diagonaux de la self-energie. De plus, l'equation de Bethe-Salpeter derivee 
n'est valide que pour les termes diagonaux (Nous expliquons en annexe C 
comment il serait possible de le faire). Ainsi, 1' approximation de la self-energie 
par la relation (2.37) sera notre approximation finale. Nous presentons au chapitre 
suivant la fonction de Green et la self-energie obtenues par les relations (2.38) et 
(2.37). 



2.4 Regies de somme 

Maintenant que nous avons termine le developpement de notre 
approximation nous allons verifier si celle-ci permet de satisfaire certaines regies 
de somme ainsi que le theoreme de Mermin- Wagner. 
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2.4.1 Moments de la susceptibilite 

On peut se demander si notre result at serait le meme si au lieu de prendre 
la regie de somme (2.35) on avait determine la valeur moyenne de double 
occupation avec la regie de somme exacte suivante 2 : 

E X(<i)e- iqnv = ((1 - n T )(l - n A )> = 1 - n + (n^) (2.39) 

On peut montrer en fait que la reponse est affirmative si notre approximation 
satisfait le moment d'ordre zero de la susceptibilite de paires. En fait, nous irons 
plus loin, nous allons demontrer que notre approximation satisfait les deux 
premiers moments de la susceptibilite de paires. 



Pour la demonstration, nous allons, dans un premier temps, determiner les 
resultats exacts. Definissons la partie imaginaire de la susceptibilite de paires: 

X "(q,t) = i<[A(q,t),At(q,0)]) (2.40) 

Utilisant la representation de Lehmann [110] pour la fonction de correlation et 
pour la susceptibilite, on peut demontrer la representation spectrale suivante 

M =r±M (2 . 41) 

J-oo TT uj-iq n 



Le developpement haute frequence de x(q) permet de faire le lien entre 
cette fonction et les moments de la fonction spectrale x"{q^)- 

-1 /»00 7 -1 /»00 T 

lim X (q, iq n ) = — / — x"(q, u>) - — ^ / — wx"(q, w) + 0{\/&) (2.42) 
iq n J-oo TT (iq n y J^oo 7T 

On voit que les deux premiers termes du developpement seront donnes par les 
deux premiers moments de la partie imaginaire de la susceptibilite. On peut 

2 La difference par rapport a la regie de somme (2.35) vient du changement de la 
limite. Ici, on pose r\ — > 0+. 
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verifier que le moment d'ordre zero correspond exactement a la difference entre 
(2.35) et (2.39): 



/. 



00 rlt 1 

— X "(q,u;)=([A(q),At(q)]> = l-n 



w . . x (2-43) 

— >0+ — v / 



In 



Ainsi, etant donne que nous avons utilise (2.35) comme condition pour determiner 
la double occupation, on voit que si notre approximation satisfait le moment 
d'ordre zero, la double occupation que nous aurions obtenue aurait ete la meme 
si nous avions impose la regie de somme (2.39). 

Pour ce qui est du premier moment, qui est une generalisation de la regie 
de somme f au cas du modele de Hubbard attractif, il sera donne par: 
00 j 

^X ,/ (q,w)« = -<[[^A(q)],At(q)]> 

-00 n 

1 ^ ( 2 - 44 ) 
= N S( £ ( k ) + £ (" k + 0) - 2m + C0(1 - 2n(k)) 

k 

oil n(k) = X^ C r(' i / ; i(k)^ ; cr (k)) est la densite d'etats renormalisee. 

On voit de par les expressions (2.43) et (2.44) que la principale 
contribution a ces sommes vient du comportement haute frequence de la fonction 
de correlation. C'est pourquoi, afin de verifier si notre approximation satisfait les 
moments, nous procedons au developpement en l/iq n de l'expression (2.33). 
Regardons dans un premier temps le comportement de Qui es t donnee par 

l'equation (2.34). 

Um tfW)- lim iE ' ,W y'; t t ?i'' 

g„->oo u g„^oo N iq n - e(k) - e(— k + q) 

= ^ + 7-^4 E( 2 ^( £ ( k )) " !)( £ (k) + £ (" k + «!)) + °(V<£) 

ou e(k) = e(k) — fi^ + et n F est la fonction de distribution de Fermi-Dirac. 
Ce resultat vient du fait que le potentiel chimique fj,^ est ajuste de fagon a ce 
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que ait le bon remplissage n. On peut en effet verifier que si le potentiel 



chimique est ainsi ajuste, satisfait la condition: 

U ™ Jf S ^o 1} (9) (e~ iqnn - e iqn ^ =l-n (2.45) 



On peut maintenant montrer en developpant le denominateur apparaissant 
dans x* 1 ) que la limite pour la fonction de correlation renormalisee sera: 

9 lim o X (1) (q,^)= ! ^ 

+ (^j2 E( 2 ^( £ ( k )) - x )( £ ( k ) + £ (- k + 1)) - (! - ™) 2r(1) ) + °(v^) 

(2.46) 

Comparant la relation (2.42) au resultat (2.46), on peut evaluer les deux premiers 
moments de la susceptibilite. Pout le moment d'orde zero on a 



" ^x" (1) (k,u,) = l-n 



On satisfait done de fagon exacte le moment d'ordre zero de la susceptibilite de 
paires. On voit qu'une condition importante vient du comportement a haute 
frequence du numerateur dans la forme de type RPA que prend la susceptibilite 
dans cette approximation. 



Pour ce qui est du premier moment le resultat est le suivant: 

/OO 7 -1 
— X //(1) (q, u)u> ~ J> ~ 2n F ( e (k)))(e(k) + e(k) - - E«)) 

-oo " k 

+ (i-n) 2 r« 

En comparant avec le resultat exact (2.44), on voit qu'on satisfait de fagon 
autocoherente ce moment car dans notre approximation la fonction de distribution 
n(k) est reduite a la distribution de Fermi-Dirac. De plus, on utilise le fait que 
dans cette approximation le potentiel chimique est approxime par n ~ /i*- 1 - 1 et que 
la correction du potentiel chimique provenant de notre approximation de la 
self-energie (2.29) permet de satisfaire la condition: 

U = 2E (1) + (1 - n)r (1) (2.47) 

Done, la demonstration est faite que le comportement hautes frequences de la 
fonction de correlation sera tel qu'il permet de satisfaire notre generalisation de 
la regie de somme f. 
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2.4.2 Comparaison avec l'approche des moments 

II est interessant de noter la ressemblance entre la methode utilisee pour 
determiner le vertex et les approches de type approximation par les moments 
[57,74]. Dans ces approches, on considere une approximation de type n-poles pour 
la fonction de Green. Ce qui correspond a une approximation ou apparaissent 2n 
parametres associes a la position et au poids de chaque pole. Ces parametres sont 
determines par la connaissance de 2n moments spectraux. Une telle realisation 
est possible par le fait que les moments du poids spectral contiennent toute 
l'information necessaire pour construire la fonction de Green. Lonke [57] a montre 
que la connaissance de n de ces moments permet de definir un domaine ferme du 
plan complexe auquel appartiendra la fonction de Green. 

Le probleme auquel nous nous sommes confronted pour la determination de 
la fonction de correlation est different pour trois raisons. Primo, dans notre cas, 
on ne tronque pas le nombre de poles. Nous utilisons les poles que nous pouvons 
determiner de fagon exacte pour la fonction x ( u; ) (exacte pour une valeur de /j, 
donnee). Les poles de x seront differents de ceux de Xo ^ cause du parametre T 
que nous avons introduit. Mais nous gardons un nombre infini de poles. Ce qui 
nous donne un certain avantage sur la methode habituelle de l'approche des 
moments. 

Secundo, les moments de la susceptibilite de paires ne sont pas connus de 
fagon exacte, contrairement a ceux du poids spectral a une particule. II faut done 
ajouter d'autres relations pour fermer le systeme d'equations. Ici, on a utilise la 
relation (2.30) que Ton a derivee d'une approximation du type Baym-Kadanoff. 

Tertio, la susceptibilite de paires n'est pas un poids spectral car elle n'est 
pas positive dans tout son domaine de definition. En effet, elle est positive pour 
les frequences positives et negative pour les frequences negatives, e'est-a-dire 
qu'elle mene a une dissipation positive. Si on voulait utiliser la methode des 
moments, il serait preferable de definir une fonction de distribution qui soit 
positive. 
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Ann de s'en definir une, considerons la representation spectrale de la 
fonction de correlation en temps imaginaire. 

x(q^)=y_ oo -x (q-^ir^ 

da; X "(q,") e ( 2-)- 
J_ 00 2vr tanh(/3w/2) cosh(/?w/2) 

= 7_ 00 2^^ (q ' W) cosh(/5u;/2) 

De par la positivite de la dissipation, x"(q, w) = — x"^ ~ w )? u es t facile de 
montrer que le poids spectral, A x est positif pour toutes les frequences. De plus, 
on sait que la susceptibilite s'annule a frequence nulle. On en deduit que A x sera 
integrable a ce point. On peut montrer que le premier moment de ce poids 
correspond a la condition que nous avions posee pour determiner le vertex 
irreductible. Plus precisement, nous obtenons la somme de (2.35) et de (2.39), 
mais comme mentionne precedemment si Ton satisfait Tune de ces conditions 
l'autre est satisfaite aussi. 

/ ^-A x (d,u) =T^ X (q,,gj(e-^ 0+ +e-^°-) 

= (n T n i ) + ((l-n T )(l-n i )) 

Ainsi, on peut presenter la regie de somme choisie pour determiner la valeur 
moyenne de double occupation (2.35), comme une fagon de determiner le 
deplacement des poles de x^ P ar rapport a ceux de Xo^ en imposant le respect 
du premier moment du poids spectral A % . 



2.4.3 Theoreme de Mermin- Wagner 

II est interessant de verifier si notre systeme d'equations permet de 
satisfaire le theoreme de Hohenberg-Mermin-Wagner-Coleman [49,77,24] selon 
lequel aucun parametre d'ordre invariant sous symetrie continue peut acquerir une 
valeur moyenne non- nulle a temperature finie dans un systeme bidimensionnel. II 
est a noter que ceci n'est pas en contradiction avec la presence d'une transition 
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de type Berezinskii-Kosterlitz-Thouless [13,61] a temperature finie car comme 
nous l'avons mentionne precedemment cette transition se fait sans apparition 
d'ordre a longue portee. 

Considerons l'expression pour le vertex irreductible (2.30), pour la fonction 
de correlation (2.33) et la regie de somme (2.35). Ces trois relations sont 
suffisantes pour demontrer le respect du theoreme de Mermin- Wagner. 

Considerons d'abord la limite haute temperature ou Ton connait 
approximativement la valeur moyenne de double occupation. On peut montrer que 
dans cette limite (n^n^) est pres de la valeur Hartree, n/4. Alors, le vertex est 
faiblement renormalise T ps U (U etant negatif le vertex le sera egalement). Done, 
si la temperature est suffisamment elevee par rapport a U, le denominateur dans 
(2.33) est toujours positif. En abaissant, la temperature Xq{q = 0) augmente, il en 
sera done de meme de x^Hq = 0) e * de (n^n^). L'expression du vertex (2.30) 
nous indique que celui-ci devrait alors diminuer (en valeur absolue). 

Pour la demonstration de l'absence de point critique, on procedera par 
l'absurde. Supposons done, qu'a une temperature donnee l'augmentation de 
X (q = 0) est suffisante pour qu'il y ait une singularity Tx (q = 0) = — 1. Alors, 
etudions le comportement autour de ce point singulier: 1 + r%(q, iq n = 0) oc q 2 . Un 
tel comportement n'est pas integrable pour un systeme bidimensionnel (etant 
donne qu'il y a au numerateur un facteur q provenant de la mesure). On obtient 
que pour toute temperature finie, la somme (2.35) ne convergera pas. Ce qui 
indiquerait un taux de double occupation divergeant. De par la relation (2.30), on 
aurait un changement de signe du vertex qui serait infiniment repulsif ce qui 
est en contradiction avec la supposition de depart. On en conclut que notre 
formalisme empeche l'existence de point critique et qu'au fur et a mesure que la 
temperature est abaissee, la renormalisation du vertex devient de plus en plus 
importante de telle fagon que sa valeur absolue soit inferieure a la valeur critique. 
La seule supposition que nous avons faite est qu'il existe bien une valeur de 
{n^n^) qui permette de satisfaire simultanement la regie de somme (2.35) et 
l'expression pour le vertex irreductible (2.30). 
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II faut porter attention au cas d'un reseau fini car alors la fonction 
Xoi^ i-Qn = 0) P ren d un nombre fini de valeurs distinctes. II est alors possible de 
trouver une valeur de F telle que pour certaines valeurs de q, le denominateur de 
(2.33) puisse etre negatif alors que la somme (2.35) demeure finie et positive (le 
denominateur ne s'annulant pour aucune des valeurs possibles de q). Le systeme 
d'equations pourrait alors etre entierement satisfait. 3 

Notre preuve du respect du theoreme de Mermin- Wagner est base sur le 
fait que dans notre approximation la susceptibilite de paires ne peut pas etre 
singuliere. Ainsi, on en deduit que notre approximation ne peut pas predire la 
presence d'un point critique BKT et ce meme hors du demi-remplissage. Ceci sera 
done une limite de notre approximation qui ne pourra pas s'appliquer dans la 
limite tres basse temperature hors demi-remplissage. 

2.4.4 Est-ce que notre approche est conservative? 

La ressemblance entre notre approche et une approche de type 
Baym-Kadanoff va plus loin que ce que nous avons mentionne a la section 2.3. 
En effet, tout comme la fonction de partition en presence d'un champ externe 
permettait de generer les termes hors diagonaux de la fonction de Green, la 
fonctionnelle W[Q] permet de generer les termes hors diagonaux de la self-energie 
(voir annexe B). Ainsi, le resultat demontre par Baym et Kadanoff [10,11], se 
generalise dans notre cas a la definition d'un vertex irreductible dont les termes 
diagonaux peuvent s'exprimer comme une double derivation d'une fonctionnelle 6 
par rapport a la fonction de Green. Pour determiner cette fonctionnelle, on doit 

3 II est a noter qu'une etude en temperature devrait alors soit montrer que le 
vertex admet une discontinuity ou qu'a basse temperature il existerait deux valeurs 
distinctes de T qui permettent de satisfaire simultanement la regie de somme (2.35) 
et l'expression du vertex irreductible (2.30) dont une ou Tx (q = 0) < 1 et l'autre ou 
rXo(? = 0) > 1- -Des resultats preliminaries obtenus par Kyung semble montrer que 
cette deuxieme solution serait observee pour des interactions intermediaires [62] . 
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considerer la partie hors diagonale de la self-energie approximee par (2.25) et 
(2.28). Au premier ordre en C/, on obtient: 

s 12 (i, 2; g) = J & jf u = a 12 (i, i)rS>*(i - 2) 

S 21 (l, 2; a) = A / 9[ f n = G 21 (l, l)r}?5(l - 2) 
<5£ 12 (2, 1) 

Or, nous avons obtenu r n = r 22 . Ainsi, il est facile de developper une 
fonctionnelle qui permet de generer les parties hors diagonales de la self-energie 
au premier ordre en Q et de generer la partie diagonale du vertex irreductible. 
Celle-ci est donnee par: 

e[G] = |Tv|^(T,2)r«(2,T, 4, 3)<?(3,4)j 

II est a noter que la fonctionnelle ne permet de generer que les termes 
non-diagonaux de la self-energie. 

II est pertinent de se demander si notre approche conserve de fagon 
globale l'energie, la quantite de mouvement et le nombre de particules. Par 
conservation globale, on veut dire que les fonctions de correlation seront telles 
qu'elles pourront bien evaluer la variation de la fonction de Green par rapport a 
l'introduction d'un champ infinitesimal approprie. Par exemple, si Ton fait une 
petite variation du champ de jauge, on peut verifier si 1' approximation satisfait la 
loi de conservation globale du nombre de particules [11]. Si Ton se refere a la 
derivation de Baym [11] pour le developpement d'une approximation conservative, 
il faut montrer que la self-energie (partie diagonale dans notre cas) peut s'ecrire 
comme la variation d'une fonctionnelle par rapport a une variation de la fonction 
de Green. Cependant, dans notre cas, le champ que nous avons introduit est 
couple a des fonction de Green anormales. Ainsi la fonction de reponse que nous 
evaluons, x(q) ne P eu t P as decrire la reponse du systeme a une variation globale 
du nombre de particules ou de l'energie. 

Pour pouvoir affirmer si notre approximation est conservative, il faut 
d'abord determiner comment nous decidons d'evaluer les fonctions de reponse 
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pertinentes. Nous pourrions considerer la relation de Bethe-Salpeter, (2.24), que 
nous avons utilisee pour evaluer la fonction de correlation de paires, et, 
presumant que les symetries de croisement sont satisfaites, evaluer les fonctions de 
reponses d'interet. Malgre la ressemblance entre notre approche et le formalisme 
prescrit par Baym, il n'y a aucune raison de penser qu'une telle approximation 
pourrait etre conservative. 

La meilleure fagon de proceder serait de considerer notre formalisme en 
presence d'un autre champ externe qui permettrait de decrire la reponse du 
systeme a une variation globale du nombre de particule (changement de jauge), 
de l'energie ou de la quantite de mouvement. Ainsi, on pourrait s'assurer d'avoir 
une estimation des fonctions de reponse qui soit conservative. Cette etude a ete 
developpee par Kyung [62] pour le calcul de la susceptibilite magnetique et la 
susceptibilite de charge du modele de Hubbard attractif. 
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COMPARAISON AUX 
DONNEES MONTE CARLO 

Maintenant que nous avons developpe une approximation pour l'etude du 
modele de Hubbard, nous pouvons verifier la validite de cette approximation et 
etudier les proprietes predites par celle-ci. Nous pouvons entre autre verifier ce 
que nous avions presume au chapitre 1, c'est-a-dire la presence d'un large regime 
avec pseudogap pres du demi-remplissage. 

Pour notre etude nous utiliserons une valeur d'interaction situe du cote des 
faibles couplages par rapport au couplage maximisant l'interaction (voir figure 3), 
soit U = —4. Ainsi, notre etude ne devrait pas presenter des effets de couplage 
fort. Malgre que le cas a demi-remplissage a deja ete etudie pour le modele 
repulsif dans le cadre de l'approximation TPSC [116] nous l'etudierons a nouveau 
car cela nous permettra de comparer les resultats du canal longitudinal avec ceux 
du canal transverse. La raison du choix de ce cas est qu'il est alors beaucoup 
plus facile de faire une etude en taille par simulation Monte Carlo quantique sans 
avoir des effets parasites tel que l'absence de points sur la surface de Fermi. De 
plus, dans notre etude nous presenterons de nouvelles quantites telle que la 
densite de suprafluide. Malgre le choix du cas demi-rempli, le lecteur pourra noter 
que les equations derivees dans le cadre de la methode TPSC ne seront pas 
restreintes a ce cas. 

Ce dernier chapitre est divise en cinq parties. Dans la premiere partie, nous 
presentons la technique Monte Carlo quantique (MCQ) utilisee pour le calcul des 
proprietes du modele de Hubbard. Nous expliquons les sources d'erreurs qui sont 
reliees a cette technique. 

La technique MCQ est utile car elle permet de verifier la validite de notre 
approximation. C'est ce qui est presente dans la deuxieme partie du chapitre, 
sections 2 et 3. Nous y demontrons que la forme de type RPA utilisee pour la 
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susceptibilite de paires est bien observee dans les resultats MCQ et que le vertex 
renormalise est bien approxime par notre approche analytique. Dans la section 3, 
nous montrons que le spectre a une particule est bien estime par notre approche 
analytique. Nous montrons que ce spectre presente un pseudogap au niveau de 
Fermi et ce a des temperatures bien superieures a la temperature critique. 

Dans la troisieme partie, section 4 et 5, nous etudions d'autres proprietes 
mesurees par MCQ arm d'etudier le mecanisme qui serait responsable de 
l'apparition du pseudogap. Comme nous le verrons, les donnees MCQ demontre 
clairement qu'aux temperatures ou le pseudogap est observe, les fluctuations de 
paires dominantes sont des fluctuations classiques. De plus, nous montrons que 
la densite de suprafluide y demeure tres faible et qu'ainsi le mecanisme a 
la base du regime pseudogap n'est pas lie a la presence d'ordre a longue 
portee. Nous argumentons que le pseudogap correspond a un presurseur de l'etat 
supraconducteur present a plus basse temperature. Ce resultat montre bien que 
les fluctuations supraconductrices a longue portee jouent un role important dans 
la disparition du liquide de Fermi. 

La partie suivante, section 6, presente l'etude du mecanisme pour la 
formation du pseudogap tel que predit par notre approximation. Nous voyons qu'il 
correspond aux proprietes que Ton a observe dans l'etude MCQ. La longueur 
de correlation de paires devient tres importante a basse temperature amenant 
l'apparition d'un regime classique renormalise. Nous montrons que lorsque cette 
longueur devient plus importante que la longueur d'onde de de Broglie, il peut 
apparaitre un pseudogap dans le spectre a une particule. 

Dans la derniere section, nous montrons que ce que nous observe par 
l'application des methodes de MCQ au modele de Hubbard demi-rempli peut 
s'extrapoler en dehors du demi-remplissage. On y observe egalement un pseudogap 
ainsi que l'apparition d'un regime classique renormalise. 
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3.1 Monte Carlo quantique 

Les methodes Monte Carlo visent a resoudre des problemes complexes 
par Putilisation des nombres aleatoires et des statistiques. De fagon generale, 
elles consistent a resoudre des problemes d'integration ou de somme de 
fagon partielle en selectionnant par un processus aleatoire un certains nombre 
d'elements de la somme. C'est pourquoi ces methodes font parties d'une branche 
des mathematiques dite experiment ale. Etant donne que Ton ne fait qu'un 
echantillonage plus ou moins representatif de la totalite des observations possibles, 
la reponse obtenue par Papplication des methodes Monte Carlo est statistique 
dans sa nature et est done entachee d'une certaine incertitude. Cette derniere 
donne a l'experimentateur une mesure de la vraisemblance de la reponse observee. 
Comme ouvrage general sur les techniques Monte Carlo nous recommandons la 
lecture des ouvrages de Hammersley et Handscomb [46], de Kalos et Whitlock 
[58] et de Robert et Casella [95]. 



En mecanique statistique, il est courant d'utiliser les methodes Monte Carlo 
pour revaluation des moyennes thermodynamiques. II faut dans un premier temps 
choisir des variables aleatoires permettant d'explorer l'espace des phases. II faut 
s'assurer que le le domaines de valeurs admises par les variables permettent 
d'explorer l'ensemble de l'espace des phases de fagon a permettre le calcul des 
moyennes thermodynamiques d'interet. Pour un systeme quantique, il n'est pas 
simple d'y parvenir. Nous utilisons le developpement de Trotter [111] qui permet 
de transformer notre systeme quantique bidimensionnel en un systeme classique 
tridimensionnel: 

e -/?(«„+Wi) =exp - IH 'dT(H +Hj) 



= lim If 



-ArH„ e -ArH; 

aC^ooAA- (3.1) 



1=1 

N, 



:f[(e- Arn °e- ATH ' +0(At 2 )) 
i=i ^ ' 



ou N t At = (3. La troisieme dimension correspond au temps imaginaire. La 
discretisation de ce temps introduit une erreur de l'ordre de At 2 dans le calcul 
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des moyennes thermodynamiques [41]. Typiquement, nous choisirons At ~ 1/10, ce 
qui implique que cette approximation induira une erreur systematique de l'ordre 
de 1% qui sera environ, pour notre echantillonage, du meme ordre de grandeur 
ou inferieure a l'erreur statistique. On voit de notre resultat (3.1), que la 
dimension supplemental sera discrete avec N T points. 

Ayant ajoute cette dimension, on peut introduire la variable aleatoire 
qui nous permettra d'explorer l'espace des phases. On peut y parvenir par 
l'utilisation d'une transformation de Hubbard-Stratonovich ou la variable auxiliaire 
[82] correspond a la variable aleatoire. Appliquant cette transformation de fagon 
a decoupler le terme d'interaction, H l7 on peut alors faire la trace sur les 
fermions de fagon analytique. La somme sur les configurations du champ aleatoire 
sera quant a elle faite par la methode Monte Carlo. Pour l'application de la 
transformation de Hubbard-Stratonovich au modele de Hubbard, nous referons le 
lecteur aux papiers originaux [14,48,118] ainsi qu'a la presentation de Tremblay 
[109]. Le lecteur y trouvera aussi l'expression permettant d'evaluer les observables 
en fonction des configurations du champ aleatoire x i ou i correspond aux points 
du reseau tridimensionnel espace et temps imaginaire. II pourra noter que la 
variable aleatoire choisie permet d'evaluer directement un nombre important 
d'observables. Les details de l'application des methodes Monte Carlo au modele 
de Hubbard peuvent aussi etre retrouves dans le rapport technique de Touchette 
et Poulin [108]. 

A l'erreur systematique s'ajoute l'erreur statistique mentionnee 
precedemment. Celle-ci est liee a l'efficience de la simulation a explorer l'espace 
des phases. Ann de connaitre la vraisemblance de nos simulations, il sera 
important d'evaluer l'erreur statistique. L'importance de cette erreur depend, entre 
autres, de la technique d' echantillonage choisie. Par exemple, si Ton fait N 
observations d'une observable O, la moyenne de notre echantillonage sera donnee 
par 

„ _ Eili QjPi 
Ei=i Pi 
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ou p i est la densite de probability associee a l'evenement Si les observations 
sont independantes les unes des autres, il est alors justifie d'utiliser comme 
estimateur de l'incertitude sur la moyenne la variance de cette observable reduit 
par un facteur lie au nombre d'observations: 

j2_ 1 / E,1iO, 2 p. 



Dans notre cas, on utilise un processus du genre marche aleatoire (chaine 
de Markov) ou le point de l'espace des phases considere a une etape donnee est 
entierement determine par le point precedent. Dans un tel processus, il y a une 
forte correlation entre des observations rapprochees. Cette correlation entre les 
differentes mesures complique l'estimation de la variance d'une observable. Ann de 
se ramener a des observables independantes, on regroupe les N observations en 
sous-groupes de M observations (N etant un multiple de M). On peut alors 
definir de nouvelles variables: 

2^j = l+(l-l)MPj 

Si M est suffisamment grand par rapport a la longueur des correlations entre les 
mesures, on peut alors supposer que les N/M observations de la variable O sont 
independantes les unes des autres. Dans ces conditions, nous pouvons utiliser les 
estimateurs suivants pour la moyenne de l'observable O et pour son incertitude: 



2^i=i Pi l^i=\Pi 



V=T77^-r[ E '-^°! P ' -<' a 



N/M-l\ J2»Jm Pi 
ou, on a introduit une nouvelle densite de probabilite 

IM 

Pi= E Pi> le[\,N/M\ 

j=l+(l-l)M 

Typiquement, les correlations entre les mesures couvrent 50 mesures pour une 
temperature de T = 1/5 (moins etendues a plus haute temperature). C'est 
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pourquoi nous avons regroupe les mesures en blocs de 125 a 350 mesures. II est a 
noter que revaluation de la probabilite p i (ou de Pi j) n'est pas faite de fagon 
explicite. Celle-ci est partiellement compensee par l'utilisation d'un echantillonage 
biaisee (voir [108]). 

En terminant, mentionnons qu'il existe un certain nombre de conditions ou 
il est difficile de faire des echantillonages efficaces et done ou la convergence vers 
la moyenne statistique est lente. Une premiere de ces conditions est liee a ce que 
Ton appelle le probleme de signe [109]. Ce probleme est observe parfois lorsque la 
convergence de revaluation de la fonction de partition est lente et que celle-ci 
prend une valeur presque nulle. Cependant, avec la variable aleatoire que nous 
utilisons, ce probleme ne sera pas present dans le modele de Hubbard attractif. 
Un autre probleme important est le probleme du collement. Dans le modele de 
Hubbard on l'observe a basse temperature et pour de forts couplages. II resulte 
de l'incapacite a bien explorer l'espace des phases, la variable aleatoire prenant 
des valeurs toutes situees dans une meme region. Pour surmonter ce probleme, 
nous faisons plusieurs simulations avec des conditions initiales differentes ce qui 
devrait permettre, pour un nombre de simulations suffisamment eleve, d'explorer 
une bonne partie de l'espace des phases. 

3.2 Validation des approximations de la methode analytique 

Maintenant que nous avons decrit brievement la technique Monte Carlo 
et ses sources d'erreurs ou d'incertitude, nous pouvons comparer les resultats 
obtenus lorsque appliquee au modele de Hubbard avec les resultats de l'approche 
a N-corps developpee au chapitre 2. Dans le reste de ce chapitre nous refererons 
aux resultats Monte Carlo quantique par 1' abbreviation MCQ et a l'approche 
analytique du chapitre 2 par l'abbreviation TPSC ("Two-particle Self-consistent"). 

Dans un premier temps, nous allons verifier la validite des approximations 
faites. La premiere approximation a verifier est l'approximation du vertex 
irreductible. Est-il approprie de le supposer constant? Pour le verifier regardons la 
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fonction de correlation obtenue par simulation Monte Carlo et comparons-la 
a une approximation de type RPA avec une interaction renormalisee. Nous 
presentons a la figure 5 la susceptibilite a frequence nulle en fonction du vecteur 
d'onde. Les donnees proviennent d'un cas avec interaction intermediaire U = —4, a 
demi-remplissage n = 1 et l'inverse de la temperature j3 = 4. Le vertex renormalise, 
T, a ete ajuste de fagon a ce que la valeur a q = (0,0) obtenue de la forme 
RPA, representee par une croix sur la figure, concorde avec celle mesuree par 
Monte Carlo, representee par un carre. On voit que cette valeur du vertex permet 
de bien retrouver la susceptibilite pres de q = (n,n) = (4,4). Cependant, pour les 
points au voisinage de (0,0), les donnees Monte Carlo seraient plus approchees 
avec un vertex plus eleve en valeur absolue. On en deduit qu'un vertex avec une 
dependance vectorielle telle que sa valeur est plus petite pres de q = (0, 0) 
qu'ailleurs pourrait donner une meilleure estimation de la susceptibilite de paires. 

A la figure 6, on presente la dependance en frequences de Matsubara de 
la susceptibilite de paires uniforme, q = (0, 0) dans les memes conditions. Les 
donnees Monte Carlo ont ete obtenues par une transformation de Fourier des 
resultats obtenus en temps imaginaire. II faut noter que les barres d'incertitude 
sur ces donnees sont approximatives car nous n'avons pas mesure la fonction de 
covariance de la susceptibilite en temps imaginaire. On remarque que la valeur du 
vertex ajustee par le resultat MCQ a frequence nulle est trop petite pour les 
frequences voisines mais a peu d'effet a haute frequence. II est a noter que l'ecart 
observe est relativement insignifiant par rapport a la valeur a frequence nulle qui 
n'a pas ete placee sur la figure car elle est environ 20 fois plus importante que la 
frequence voisine. Le fait que la susceptibilite a frequence nulle soit beaucoup plus 
importante qu'aux autres frequences de Matsubara est du au fait que le systeme 
est dans un regime classique renormalise comme nous le verrons a la section 3.4. 

On deduit de nos resultats qu'une approximation de type RPA avec un 
vertex constant n'est pas si mauvaise. II est alors interessant de comparer la 
valeur du vertex ajustee par valeur de la susceptibilite uniforme et statique 
mesuree par MCQ et celle obtenue par notre approche. Dans le tableau 1, on 
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Tableau 1: Comparaison de la valeur du vertex 
choisie arm que la susceptibilite de paires d'une 
approximation du type RPA concorde avec les 
resultats Monte Carlo (MCQ) avec celle obtenue par 
1' approximation proposee (TPSC) pour un reseau de 
taille 8x8 a demi-remplissage, n = 1 et U = —4. 



p 


MCQ 


TPSC 


3 


-2.084 


-2.1183 


4 


-1.940 


-1.9977 


5 


-1.784 


-1.8189 



presente les resultats pour un reseau 8x8 a demi-remplissage avec interaction 
intermediaire pour trois temperatures differentes. 

L'accord est assez bon. Cependant, on doit faire remarquer que dans la 
forme RPA, lorsque la valeur du vertex devient telle que le denominateur est 
petit a faible frequence, une petite variation sur la valeur du vertex peut avoir 
un effet important sur la susceptibilite de paires. II est done interessant de 
comparer d'autres donnees telle la valeur moyenne de double occupation. Etant 
donne que revaluation de celle-ci est beaucoup plus rapide que celle de la 
susceptibilite de paires, nous presentons au tableau 2 la valeur mesuree pour une 
plus grande diversite de conditions. 

A premiere vue, il semble y avoir un assez bon accord entre ces valeurs. 
Pour les cas a plus faible remplissage, on peut denoter que l'accord n'est pas 
aussi bon a basse temperature. En terminant, on peut remarquer que les donnees 
MCQ presentees au tableau 2 sont tres precises, en fait, l'erreur statistique est 
inferieure a l'erreur systematique qui doit etre de l'ordre de 1% (car At = 0.1). 
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Tableau 2: Comparaison des valeurs moyennes 
de double occupation predites par 1' approximation 
proposee (TPSC) et celles mesurees par Monte Carlo 
(MCQ) pour des reseaux 8x8. Les chiffres entre 
parentheses servent a indiquer l'incertitude sur le 
dernier chiffre significatif. 



u 


n 


P 


MCQ 


TPSC 


-4 


1. 


2 


0.3675(2) 


0.364 


-4 


1. 


3 


0.3671(2) 


0.370 


-4 


1. 


4 


0.3689(2) 


0.377 


-4 


1. 


5 


0.380(5) 


0.386 


-4 


0.8 


4 


0.2671(1) 


0.273 


-4 


0.8 


6 


0.2703(1) 


0.288 


-4 


0.55 


4 


0.1561(1) 


0.157 


-4 


0.55 


6 


0.1590(1) 


0.169 


-6 


1. 


3 


0.4227(5) 


0.4038 


-6 


1. 


4 


0.4258(8) 


0.4140 


-6 


1. 


5 


0.4267(4) 


0.4233 



3.3 Proprietes a une particule 

3.3.1 Poids de quasiparticules 

Pour l'etude des proprietes a une particule, une observation directe des 
donnees MCQ nous montre que notre systeme n'a pas le comportement habituel 
des liquides de Fermi. Cette observation decoule de la mesure du poids de 
quasiparticule tel qu'introduit par Vilk et Tremblay [116]. Generalement, lorsqu'on 
parle du poids de quasiparticule on parle du saut observe a temperature nulle 
dans la fonction de distribution au niveau de Fermi. Luttinger a propose une 
premiere generalisation a temperature finie qui consiste en le residu du pole de la 
fonction de Green au niveau de Fermi. Cette quantite etant difficile a extraire des 
simulations MCQ, Vilk et Tremblay ont propose de prendre une mesure du poids 
spectral a basse frequence donnee par: 

m - r P = "2G(k,r = 0/2) (3-2) 

V ; 2tt cosh(/?u;/2) V f ' J y J 
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La derniere egalite se demontre par la representation spectrale de la 
fonction de Green en temps imaginaire (3.3). On peut montrer que pour un 
liquide de Fermi, on retrouve a temperature nulle la definition usuelle du poids 
de quasiparticules. 

On a represents a la figure (7) le poids spectral a basse energie, z, 
mesure par MCQ en fonction de la temperature. On a utilise les conditions 
U = —4, n = 0.95 et At = 1/10. II est possible de verifier que le resultat est 
sensiblement identique a demi-remplissage. Ann de bien montrer l'effet de taille 
nous avons presente a basse temperature les resultats obtenus pour des systemes 
de differentes tailles (voir les indications dans la legende). Dans la partie du haut 
de la meme figure nous avons presente le facteur de structure de paires au 
vecteur d'onde q=(0, 0). Celui-ci a egalement ete obtenue par simulation MCQ. 
Pour chaque point, nous avons fait entre 5000 et 20000 mesures (plus de mesures 
a plus basse temperature et pour les plus petits systemes). 

Pour un liquide de Fermi, on devrait observer une augmentation du poids 
spectral au niveau de Fermi, z(k F ), avec l'abaissement de la temperature 
indiquant que le poids spectral est pique pres des frequences nulles a basse 
temperature. Cependant, les donnees presentees a la figure 7 montrent que le 
poids spectral diminue lorsque Ton abaisse la temperature demontrant que notre 
systeme n'est pas du type liquide de Fermi. On observe que l'abaissement du 
poids spectral s'accompagne d'une augmentation des correlations de paires. II 
est done normal de penser que ce sont les fluctuations de paires qui sont 
responsables de la disparition du liquide de Fermi et que s'il y a formation 
d'un pseudogap dans le systeme celui-ci corresponde a un precurseur de l'etat 
supraconducteur present a temperature nulle. 

3.3.2 Spectre a une particule 

L'abaissement de poids spectral a basse frequence accompagnant une 
diminution de temperature etant revelateur d'un comportement different des 
liquides de Fermi, n'est cependant pas suffisant pour demontrer la presence d'un 
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pseudogap. Pour le demontrer il faut obtenir le poids spectral. Celui-ci s'obtient 
par prolongement analytique de la fonction de Green en temps imaginaire. Une 
des techniques les plus efficaces pour le prolongement analytique des donnees 
MCQ est la methode de maximisation d'entropie (MME) decrite en annexe D. 
Celle-ci se fait par Futilisation de la representation spectrale de la fonction de 
Green en temps imaginaire: 

Nous presentons aux figures 8 et 9 le poids spectral obtenu par 
le prolongement des donnees MCQ (partie gauche de chaque figure) et par 
prolongement des donnees obtenues par notre approche analytique (partie droite 
de chaque figure). Les resultats analytiques ont ete obtenu par Futilisation des 
equation de et de E^ 2 \ equation (2.38) et (2.37). Pour ces resultats nous 

avons applique la meme technique de prolongement analytique que pour les 
donnees MCQ. De plus, etant donne que cette technique est sensible a la 
precision des donnees nous avons ajoute un bruit gaussien dans nos donnees 
analytiques permettant ainsi d'ajuster leur precision a celle des simulations MCQ. 
II est a noter que l'ajout de ce bruit diminue le nombre de details que 
Ton peut observer dans le spectre et qu'un bruit trop important pourrait 
eventuellement camoufler l'existence d'un pseudogap. C'est pourquoi nous avons 
fait des simulations MCQ tres precises avec de 125 000 a 160 000 mesures pour 
chaque vecteur d'onde et chaque temperature. Pour chaque poids spectral presente 
nous avons indique le vecteur d'onde correspondant du cote gauche du graphique. 

On observe qu'a la temperature la plus elevee aucun pseudogap n'est 
observe et que la relation de dispersion n'est que legerement modifiee par rapport 
a celle des electrons libres. Les differences deviennent plus importantes au fur et 
a mesure que Ton s'eloigne du niveau de Fermi. Lorsqu'on abaisse la temperature 
un pseudogap apparait au niveau de Fermi. Ce pseudogap correspond a un 
minimum local a basse energie tel qu'on observe au centre de la figure 9. Le 
pseudogap est plus important pres du vecteur d'onde (0, tt) que dans la region du 
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vecteur (ir/2, 7r/2). Cette caracteristique n'est pas reliee a la symetrie du gap qui 
est bien de type "s-wave" dans le modele de Hubbard attractif mais est due a la 
variation de la vitesse de Fermi. Comme on l'expliquera a la section 3.6.2 il est 
plus facile de creer un pseudogap dans une region de la surface de Fermi 
caracterisee par une plus petite vitesse de Fermi. 

II faut remarquer que la discretisation du temps imaginaire limite notre 
capacite de prolonger le poids spectral a haute frequence. Pour s'en convaincre, il 
suffit d'observer la forme du noyau dans la relation (3.3). De par sa forme, on 
conclut que tres peu de points en temps imaginaire nous permettent de deduire 
le poids spectral a haute energie, d'ou la difficulte. 

On peut remarquer que l'abaissement du poids de quasiparticules precede 
l'ouverture du pseudogap. Ce resultat nous indique qu'avant l'apparition du 
pseudogap, le pic dans le poids spectral au niveau de Fermi (voir par exemple 
(3 = 2 et k = (0, 7r)) s'elargit avec l'abaissement de temperature. Ce phenomene 
semble etre plus prononce pres du niveau de Fermi. Ainsi, on pourrait voir 
l'apparition d'un pseudogap dans la densite totale d'etats avant d'en voir dans 
le poids spectral. A la figure 10, on presente la densite d'etats obtenus par 
maximisation d'entropie. La partie de gauche de la figure nous montre le resultat 
obtenu des donnees MCQ et la partie de droite celui de notre approche 
analytique. La densite d'etats, D(lu) etant simplement la somme sur les vecteurs 
d'onde du poids spectral, on peut l'obtenir par l'application de la methode de 
maximisation d'entropie a la formule: 



ou N est le nombre de vecteur d'onde du reseau. Pour l'estimation de 
l'incertitude sur les points en temps imaginaire, nous avons suppose que la 
matrice de covariance etait diagonale dans l'espace de Brillouin. 

On voit qu'au fur et a mesure que Ton abaisse la temperature (courbe du 
bas a celle du haut) un pseudogap apparait dans la densite d'etats a basse 
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frequence. Celui-ci apparait a des temperatures ou le poids spectral ne montre pas 
de signe de la presence d'un pseudogap. II est a noter que la dependance en taille 
de cet effet n'est pas completement verifiee est qu'une etude plus detaillee devrait 
etre publiee sous peu [62]. Differentes raisons peuvent expliquer le fait que le 
pseudogap est legerement surestime par la methode analytique. Premierement, 
il n'est pas toujours facile d'ajuster le bruit des donnees analytiques a celui 
des donnees MCQ particulierement pour la densite d'etats ou on doit sommer 
la fonction de Green sur les differents vecteurs d'onde avant le prolongement 
analytique. De plus, une comparaison des fonctions de Green obtenues des deux 
approches en temps imaginaire semble montrer que notre approche analytique 
surestime l'importance des correlations (voir figure 12). 

3.3.3 Etude en taille 

Avant de terminer cette section, il nous faut parler de l'effet de differents 
facteurs sur la presence ou l'absence d'un pseudogap dans le poids spectral. II y 
a un facteur que nous avons deja mentionne soit l'effet de la precision de 
nos mesures. On montre a la figure 11 l'effet de cette precision. Pour bien 
demontrer le fait qu'un manque de precision attenue les details dans le poids 
spectral, nous avons considere la fonction de Green au vecteur d'onde (0, 7r) 
tel qu'obtenue par MCQ a temperature T = 1/4. Nous avons ajoute un bruit 
gaussien plus ou moins important a nos donnees et avons fait le prolongement 
analytique par MME. Dans la partie de gauche de la figure nous montrons la 
fonction de Green utilisee pour chacun des prolongement analytique, le bruit 
etant de plus en plus important du bas au haut de la figure. La premiere courbe 
correspond aux donnees obtenus par nos simulations sans l'ajout de bruit. La 
partie droite montre les poids spectraux correspondant. On voit que plus le bruit 
est important, moins le pseudogap est prononce. Ce resultat explique pourquoi 
lors de comparaison des resultats MCQ et ceux de la methode analytique, il est 
important d'ajouter du bruit dans les donnees de la methode analytique. De plus, 
ce resultat renforce notre conclusion selon laquelle il y a bien un pseudogap dans 
le poids spectral du modele de Hubbard. 
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Cependant, une etude en taille montre que le pseudogap est moins 
prononce lorsque Ton considere un reseau de taille plus importante. Ainsi, il est 
important de faire une etude en taille afin de predire s'il existe ou non un 
pseudogap dans la limite thermodynamique. White et al. [118] ont conclu que 
l'effet de taille finie serait responsable de la presence du pseudogap et que pour 
un reseau de taille suffisamment importante le pseudogap disparaitrait. Cependant, 
ils ont neglige l'effet de la precision des mesures et les donnees pour les reseaux 
de grandes tailles etaient peu precises. Nous presentons ici une etude plus precise 
de revolution du pseudogap avec la taille du reseau. 

A la figure 12, nous presentons la fonction de Green en temps imaginaire a 
temperature T = 1/4 au vecteur d'onde k= (0, 7r) pour differentes tailles de 
reseau. La partie de gauche nous montre les resultats MCQ et la partie de droite 
le calcul analytique (TPSC). La convergence etant plus rapide pour les plus 
grands systemes 1 , le nombre de mesures necessaires en MCQ depend de cette 
taille: 300 000 mesures pour le 6 x 6, 150 000 pour le 8 x 8 et 90 000 pour le 
10 x 10. Nous voyons que l'accord entre les resultats MCQ et notre approche n'est 
pas parfait, la fonction de Green du calcul analytique etant generalement plus 
petite en valeur absolue a r = (3/2 = 2 que celle obtenue par MCQ 2 . Cependant, 
revolution en taille semble etre sensiblement identique (par exemple, on peut 
comparer les reseaux de taille 8 x 8 et 10 x 10 de part et d'autre). 

Afin de s'assurer que revolution en taille soit bien evaluee nous avons 
ajuste le bruit dans les donnees MCQ et les resultats de la methode analytique a 
un ordre de grandeur moyen de l'ordre de 0.002 ~ 0.5% — 1% identique pour 

1 La convergence est plus rapide pour les systemes de plus grandes tailles car pour 
ces systemes il y a plus de configurations distinctes dans une plage donne d'energie. De 

plus, pour de plus grands systemes le probleme de collement semble moins important. 

2 Nous avons montre que l'accord pouvait etre ameliore lorsque le calcul est fait 

dans le canal longitudinal. Ainsi, on separe le calcul de la self-energie entre les contribu- 
tions provenant d'une susceptibilite critique et d'une susceptibilite non-critique. Pour 
la separation entre les differentes contributions, la forme symetrisee donne les meilleurs 
resultats a demi-remplissage [80] 
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toutes les tailles de reseau considerees. L'ajout de bruit se fait par l'utilisation 
d'un generateur de nombre aleatoire pour generer un bruit gaussien. L'ordre de 
grandeur du bruit est environ trois fois plus important que l'estimation de 
l'incertitude sur les simulations les moins precises. Nous avons fait ce choix car 
l'autocorrelation entre les mesures peut induire une sous-estimation des erreurs 
sur les donnees MCQ. Pour s'assurer de la precision respective de nos donnees, 
nous avons evalue pour chacune le x 2 entre la fonction de Green correspondant 
au poids spectral obtenu et la fonction de Green de depart. Le resultat etant du 
meme ordre de grandeur pour chaque taille du reseau nous assure ainsi que le 
bruit etait de la meme importance. 

Regardons maintenant ce que Ton obtient suite au prolongement analytique. 
La figure 13 nous montre revolution en taille du poids spectral au vecteur d'onde 
k= (0, 7r) et a une temperature T= 1/4. Comme a l'habitude, la partie de 
gauche montre les resultats MCQ et la partie de droite le traitement analytique. 
Une premiere observation est que l'approche analytique surestime l'importance du 
pseudogap. Cependant, revolution avec la taille du reseau semble se faire de fagon 
sensiblement identique dans les deux cas. Ce resultat suggere que le pseudogap 
n'est pas du a un effet de taille finie mais serait bien present dans le poids 
spectral et ce a des temperatures bien superieures a la temperature critique. De 
plus, il est possible d'evaluer le poids spectral dans la limite thermodynamique 
pour notre approche analytique. Tel que publie [80], ce resultat montre qu'il y a 
bien un pseudogap dans cette limite a une temperature T = 1/5. 

3.4 Regime classique renormalise 

3.4.1 Comportement critique de la susceptibilite de paires 

Comme nous avons vu, l'abaissement du poids de quasiparticule est 
accompagne d'une hausse du facteur de structure de paires. C'est pourquoi nous 
allons detailler l'etude du comportement de la susceptibilite de paires dans le 
regime de temperature ou le pseudogap est observe. Une observation directe des 
resultats Monte Carlo nous indique que la susceptibilite de paires devient tres 
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piquee autour de la frequence de Matsubara q n = 0. Pour mettre cet effet en 
evidence nous pouvons observer le rapport 

T X (q,ig n = 0)/S A (q) (3.4) 

entre la susceptibilite de paires a q n = et le facteur de structure de paires 
donne par la relation suivante: 

5 A (q) = T lim £ X (q, iqje™ (3.5) 

A la figure 14, nous presentons le resultat observe pour le cas uniforme 
q = 0. Pour chaque point, nous avons evalue entre 5000 et 30000 mesures par 
la technique MCQ. Nous presentons le comportement haute temperature en 
medaillon. On observe que les effets quantiques deviennent de plus en plus 
importants au fur et a mesure que Ton abaisse la temperature. La contribution au 
facteur de structure provenant de la composante a frequence nulle diminue ainsi 
de fagon lineaire avec l'abaissement de la temperature jusqu'a une temperature ou 
on observe un changement brusque dans ce comportement. Le minimum est 
observe a T = 1/2 soit a une temperature legerement superieure a la temperature 
ou apparait le pseudogap dans la densite d'etats. A partir de ce minimum, 
l'importance de la contribution provenant de la susceptibilite a frequence nulle 
augmente tres rapidement avec l'abaissement de la temperature. Ce comportement 
represente un retour vers un regime classique a basse temperature. A tres basse 
temperature, on voit l'apparition d'un maximum puis un retour vers un regime 
quantique. On peut remarquer que la temperature a laquelle apparait ce 
maximum depend de la taille du systeme. On argumentera a la section 3.6.3 que 
cet effet est probablement du a un effet de taille finie et qu'a tres basse 
temperature, lorsque la longueur de correlation est du meme ordre de grandeur 
que la taille du systeme, les effets quantiques se manifestent a nouveau. 

Le fait que la principale contribution au facteur de structure de paires 
viennent de la susceptibilite a frequence de Matusbara nulle pour T < 1/2 est 
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une caracteristique qui montre que, pour ce qui est des proprietes a deux 
particules, le systeme se comporte comme s'il etait classique 3 . En fait, nous 
qualifierons de regime classique renormalise le regime basse temperature ou le 
rapport Tx(q = 0)/5 A (q = 0) augmente en diminuant la temperature. 

On observe un comportement critique classique lorsque l'ecart entre les 
niveaux d'energie est petit par rapport a la temperature. Ainsi, le systeme se 
comporte comme s'il y avait un continuum d'etats disponibles, comme s'il etait 
classique [21]. On precise que le regime est renormalise pour indiquer que les 
parametres d'ecrivant la physique du systeme sont differents de ceux apparaissant 
dans le modele et done de ceux que Ton observe a haute temperature. 

3.4.2 Frequences caracteristiques 

II y a une deuxieme quantite importante qui peut nous reveler que le 
regime ou on observe le pseudogap est un regime classique renormalise, e'est la 
frequence caracteristique des fluctuations de paires. Celle-ci etant mesuree en 
frequence reelle, on doit encore une fois faire appel au prolongement analytique 
pour mesurer cette quantite. On precede de la meme fagon que pour le calcul du 
poids spectral a une particule: on utilise une methode de maximisation d'entropie 
pour deduire la partie imaginaire de la susceptibilite en frequence reelle de la 
susceptibilite en temps imaginaire mesuree par MCQ. Ces deux quantites sont 
reliees par une relation spectrale : 



3 On peut le comprendre par le fait que si le systeme est classique on peut supposer 
que l'operateur A commute avec l'hamiltonien. Dans ces conditions, cet operateur est 
independant du temps imaginaire A(r) = A et la susceptibilite n'est non nulle que 

pour la frequence de Matsubara nulle, iq n ) ^ $n,o 

4 En utilisant cette relation, on peut montrer que si la frequence caracteristique 

est petite par rapport a la temperature alors x( r ) es t a peu pres constante et le rapport 
(3.4) vaut pres de 1. 




63 



Comparaison aux donnees Monte Carlo 



Contrairement de la fonction de Green a une particule, le noyau 

apparaissant ici a un pole sur le domaine d'integration. L'integrale n'est pas pour 
autant indefmie, car on peut montrer que la partie imaginaire de la susceptibilite 
s'annule a ce pole qui est du premier ordre. Cependant, la presence de ce pole 
rend le calcul numerique difficile. On surmonte facilement ce probleme en utilisant 
le poids spectral introduit a la section 2.4.2: 



Ce poids spectral peut etre deduit de la susceptibilite en temps imaginaire par 
l'utilisation de la maximisation d'entropie et du noyau approprie. On peut par la 
suite aisement inverser la relation (3.6) pour deduire la partie imaginaire de la 
susceptibilite. Tout comme pour le poids spectral a une particule, la discretisation 
du temps imaginaire limite la capacite de deduire le poids spectral qu'au domaine 
de frequence autour de la frequence nulle. Mais comme on pourra le voir, tout 
le poids spectral associe a la susceptibilite de paires est concentre dans ces 
frequences. 

On a represente a la figure 15 la partie imaginaire de la susceptibilite 
divisee par la frequence tel qu'obtenue par l'application de la MME aux donnees 
MCQ. Nous avons utilise les conditions habituelles: U = —4, n = 1, At = 1/10 
avec un reseau de taille 8x8. Pour chaque temperature, nos simulations 
contiennent entre 125 000 et 160 000 mesures afin d'assurer une precision 
suffisamment elevee pour avoir le plus de details possibles dans le poids spectral. 

On voit qu'au fur et a mesure que Ton abaisse la temperature la 
susceptibilite de paires devient de plus en plus piquee autour de la frequence 
nulle. Voyons quel effet cela a sur la frequence caracteristique des fluctuations de 
paires. La fonction de correlation de paires est liee a la partie imaginaire de la 
susceptibilite par la relation: 



On voit que pour le cas uniforme, q = 0, x"(0, u)/ui a la forme d'une 
lorentzienne centree a frequence nulle, il est done justifie de prendre comme 




(3.6) 
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Tableau 3: Frequences caracteristiques des fluctuations 
de paires du modele de Hubbard a U = —4, n = 1 tel 
que mesure par MCQ pour un reseau 8x8 avec une 
discretisation du temps imaginaire At = 1/10. 



T 


u c 


1/3 


0.37 


1/4 


0.24 


1/5 


0.18 



mesure de la frequence caracteristique la demi-largeur a mi-hauteur de ce pic. Ce 
choix nous donne les frequences presentees au tableau 3. 

On voit qu'a haute temperature la frequence caracteristique est superieure a 
la temperature. Avec l'abaissement de la temperature, la frequence caracteristique 
diminue rapidement et devient egale a la temperature lorsque le pseudogap 
s'ouvre dans le poids spectral a une particule. A temperature plus basse, T = 1/5, 
la frequence devient inferieure a la temperature. Cependant, on doit noter que nos 
resultats sont entaches d'une certaine incertitude. Etant donne que la frequence 
caracteristique s'obtient par un prolongement analytique, il est difficile d'estimer 
cette incertitude. On devine que si nos resultats MCQ etaient plus precis, on 
obtiendrait des pics un peu plus etroit sur la figure 15. On en deduit que nos 
valeurs de frequences caracteristiques sont une limite superieure. Par ailleurs, le 
fait que nos donnees sont un peu plus precise a plus haute temperature, nous 
permet de deduire que le fait que uj c varie plus rapidement que la temperature, 
T, est bien observe dans nos simulations. 

Vilk et Tremblay [116] ont argumente que le fait que la frequence 
caracteristique devienne inferieure a la temperature est directement lie a 
l'apparition du pseudogap. Leur explication s'applique aussi bien pour le modele 
de Hubbard attractif. Pour le montrer, nous devons considerer l'expression de 
la partie imaginaire de la self-energie au niveau de Fermi en fonction de la 
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susceptibilite de paires. On peut deduire de l'equation (2.37) la self-energie 
retardee en faisant le prolongement analytique ik n — > u + irj. On obtient la relation 
suivante: 



Ti 

E(k F ,u + iri) =U- 



/ f jd (1 r j. ./ 
J^aJ — [«(<✓) + /(e(-k F + q))] ; 



X"(qy) 



u — u)' + e(-k F + q) + if] 

(3.7) 

ou e(k) est la relation de dispersion (dans notre approche analytique 
e(k) = e(k) — fi^ + T,^\ voir equation (2.31)), f(e) est la fonction de distribution 
de Fermi-Dirac et n(u) est celle de Bose-Einstein. On peut en extraire la partie 
imaginaire en procedant a l'integration sur une composante du vecteur d'onde. La 
limite rj — > + fait apparaitre une fonction delta a l'interieure des integrates pour 
la partie imaginaire de la self-energie. On peut alors integrer une des composantes 
du vecteur d'onde. Separons done le vecteur d'onde en ses deux composantes: 
et q ± . Ann de simplifier le calcul lorsque x" es t pique pres de q = 0, on 
approxime e(— k F + q) ps v F q^ . On obtient alors pour la partie imaginaire de la 
self-energie retardee: 

E£(fc F> a;)«-^ d ^[n{u:') + f^-u;)\ X ''{q^q p u') (3.8) 

ou ps (u)' — oj) jv F . Comme note par Vilk et Tremblay [116], le comportement 
de type liquide de Fermi ou non-liquide de Fermi est determine par la 
largeur relative en frequence de la fonction g(oo') = x"(q> oj')/oj' par rapport a 
celle de h(u>') = n(uj') + f(cu' — u>). Si g(uj') est a peu pres constante sur le 
domaine de frequence ou h(u>') est non-nulle, alors on obtient un comportement 
caracteristique des liquides de Fermi: 

w) oc u) 2 + (vrT) 2 (3.9) 



Cependant, lorsque la frequence caracteristique des fluctuations de paires 
est petite par rapport a la temperature, on ne peut plus approximer g(u>') par 
une constante dans le domaine ou h(cu') est non-nulle. Si la fonction g(w') 
est tellement piquee pour u/ < T que Ton puisse approximer h(uu') par son 
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comportement basse frequence T/oj', on obtient alors pour le comportement a 
basse frequence de la self-energie: 

UU pp T f d q ± . n n . (3 ' 10) 

Une contribution importante a l'integrale sur les vecteurs d'onde vient du 
voisinage de q ± ~ 0. Par nos resultats MCQ, figure 15, il semble a premiere 
vue que la partie imaginaire de la self a frequence nulle augmente avec un 
abaissement de la temperature. Ce qui demontre que celle-ci ne correspond pas a 
un comportement du type liquide de Fermi. Cependant, on peut noter qu'il y a 
un facteur T devant l'integrale dans l'equation (3.10) et que notre deduction 
dependra de la dimension du systeme (a plus haute dimension les vecteurs d'onde 
au voisinage de q ± = ont moins d'importance dans le resultat de l'integral et le 
comportement de type liquide de Fermi peut etre retrouve). 



3.5 Mecanisme pour la formation du pseudogap 

3.5.1 Densite de suprafluide 

Jusqu'a maintenant, nous avons montre que le spectre a une particule est 
caracterise par un pseudogap a une temperature superieure a la temperature 
critique, et que ce regime ou apparait le pseudogap est de type classique 
renormalise. Cependant, on peut toujours se questionner sur le mecanisme qui 
serait responsable de la formation du pseudogap: Est-ce du aux fluctuations de 
paires? A-t-on affaire a des paires locales preformees? 

La densite de suprafluide est une quantite qui est nulle dans Fetat normal 
et ne prend une valeur finie que pour des temperatures inferieures ou egales a la 
temperature critique (c'est egalement le cas pour une transition de type BKT). 
Pour ce qui est de la limite fort couplage, la densite de suprafluide prend une 
valeur finie que sous la temperature critique meme si les electrons forment des 
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paires a haute temperature (nous les appellerons paires locales preformees). A la 
temperature critique BKT, ces paires condensent et on observe alors une valeur 
finie de la densite de suprafluide. 

Pour un reseau de taille finie, on ne s'attend pas a avoir un saut tres 
prononce dans la densite de suprafluide a la temperature critique. Ce saut devrait 
plutot etre etale sur une certaine plage de temperature. Ainsi, selon la definition 
de la temperature critique, la densite de suprafluide pourrait apparaitre a une 
temperature legerement superieure a cette temperature critique. C'est pourquoi 
nous etudierons la densite de suprafluide afin de bien montrer que la temperature 
oil apparait le pseudogap est bien superieure a la temperature critique renforgant 
ainsi notre argumentation presentee par l'etude des effets de taille. 

Scalapino et al. [98] ont suggere une methode permettant d'evaluer la 
densite de suprafluide par simulation MCQ. Cette methode fait appel a la limite 
grande longueur d'onde et basse frequence de la fonction de correlation du 
courant. Utilisant la theorie de reponse lineaire, ils ont montre que la densite de 
suprafluide sera donnee par: 

Tl 

-T = A xx(Qx ~^0,q y = 0, iq n = 0) - A xx (q x = 0,q y ^ 0, iq n = 0) 

oil m* est la masse effective et A xx est la fonction de correlation courant-courant 
dans la direction x: 

KM,iQn) = f dre^-^- n \j x (l,r)j x (n,0)) 

l,n J ° 

avec le courant defini sans le facteur de charge habituel: 

j x (l, t) = J2 (^(1 + x, r)^(l, r) - ^t (1> r )^(i + ^ r) ) 

A la figure 16, on presente la fonction de correlation courant-courant telle 
que mesuree par MCQ pour differents vecteurs d'onde a frequence nulle. Cette 
quantite etant dynamique, nous avons fait 112 000 mesures pour atteindre une 
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precision suffisante: Etant donne le peu de points disponibles, il est difficile 
d'estimer les valeurs limites d'interet. Une simple extrapolation lineaire semble 
montrer que la densite de suprafluide est tres petite. Etant donne que Ton ne 
connait pas la masse effective des porteurs, il est difficile d'avoir une estimation 
quantitative de cette densite. Ce que Ton peut faire est de comparer l'importance 
du poids de suprafluide D s = n s 7re 2 /m* au poids de Drude D = n p Tre 2 /m* oil n p 
est la densite de porteur libre. Le poids de Drude s'obtient en considerant 
d'autres limites de la fonction de correlation de paires: 

n 

— = A xx(Qx -» 0,q y = 0,iq n = 0) - A xx (q x = 0,q y = 0,iq n -> 0) 

lit 

Dans la figure 16, on a represents la limite basse frequence de la fonction 
de correlation courant-courant uniforme par des cercles. Pour ces donnees nous 
avons fait une extrapolation a frequence nulle en utilisant un lissage quadratique. 
Ainsi, nous pouvons conclure que le poids de Drude est beaucoup plus important 
que le poids de suprafluide et d'apres la theorie developpee par Scalapino et al. 
[98], on aurait done affaire a un metal et non a un supraconducteur. Ce resultat 
suggere que la presence d'un pseudogap dans le poids spectral ne serait pas liee 
a la presence d'un etat ordonne. D'autre part, si Ton etait en presence de paires 
locales preformees (comme dans les theories de couplage fort), on devrait observer 
une faible densite de porteurs fibres et done un poids de Drude assez faible. 
Etant donne que Ton a choisi une interaction intermediaire U = —4, nous ne 
sommes pas tres etonnes d'observer un poids de Drude relativement important et, 
done, qu'il n'y a pas de presence de paires locales dans le regime avec pseudogap. 

3.5.2 Precurseurs de l'etat ordonne 

On peut argumenter que le pseudogap est en fait un precurseur de l'etat 
ordonne. Nous avons montre au chapitre 1 qu'a demi-remplissage le point critique 
est a temperature nulle. Cependant, pour un systeme de taille finie, lorsque la 
longueur de correlation atteint la taille du systeme, celui-ci devrait presenter 
certaines caracteristiques de la phase ordonnee. Par exemple, on peut s'attendre 
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a ce que le poids spectral presente un gap tel que prevu dans la phase 
supraconductrice. C'est pourquoi nous avons considere un systeme de taille 4x4 
a (3 = 8 afin de verifier si le pseudogap est un precurseur de l'etat ordonne. 

A la figure 17, on presente le spectre a une particule obtenue par 
prolongement analytique des donnees MCQ a basse temperature. On observe qu'il 
y a alors un gap separant les deux pics caracteristiques de l'etat ordonne. On 
voit que ces deux pics sont situes aux memes frequences que les maxima 
observes dans le regime caracterise par un pseudogap. Ce resultat suggere que le 
pseudogap soit un precurseur de l'etat supraconducteur et que celui-ci serait du 
a la presence de fluctuations supraconductrices. Sur la figure 17, on voit la 
subsistance d'un petit pic a frequence nulle. L'existence de ce pic est a verifier. 
Celui-ci pourrait etre du a des effets de taille finie. 



3.6 Longueurs caracteristiques 

3.6.1 Longueur de correlation de paires 

Nous avons observe que le pseudogap est peu sujet a 1'effet de la taille du 
systeme aux temperatures plus elevees. Malgre tout, il serait interessant de voir 
comment se comporte cette caracteristique avec un changement de la taille du 
reseau. Nous avons trop peu de donnees MCQ pour pouvoir determiner revolution 
en taille. Cependant, l'accord entre notre approche et le MCQ suggere d'utiliser le 
comportement en taille de notre approche pour predire les resultats dans la limite 
t her mo dynamique . 



Pour bien detailler notre etude en taille, il faut trouver les longueurs 
caracteristiques du systeme. Le comportement critique de la susceptibilite de 
paires nous indique qu'une premiere longueur a prendre en compte est la longueur 
de correlation de paires. Pour l'introduire, considerons notre approximation pour 
la fonction de correlation: 

x^a (3 . n) 
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On aurait un point critique si Ton pouvait obtenir 
Upp = U c = — 1/xo\q = 0). La valeur U c est appelee valeur critique du vertex. 
On sait qu'a un point critique on observe une divergence d'une longueur de 
correlation, ainsi il semble approprie de definir la longueur de correlation de 
paires par: 

e L =^j f ^ n - (3.i2) 

oil 

?0 2 Xo (q = 0) 5q 2 
On a place un indice L a la longueur de correlation pour indiquer de facon 
explicite sa dependance sur la taille du systeme. Comme on a vu, precedemment, 
a basse temperature, /? — 4, la susceptibilite de paires est piquee en frequences 
de Matsubara et en vecteurs d'onde autour du point q = (0,0,0) qui est 
caracteristique de l'approche d'un point critique supraconducteur. On peut alors 
developper la fonction de correlation autour de ce point. Utilisant la longueur de 
correlation introduite precedemment, eq. (3.12), on peut ecrire: 



u pp so 



X ( 1 )(q,,g n = 0)«-|^- (3.13) 



II est facile d'evaluer le comportement en temperature de la longueur de 
correlation dans le regime classique renormalise. Etant donne que celle-ci est 
reliee a la valeur du vertex, on doit considerer la relation utilisee pour evaluer 
ce vertex, soit l'equation (2.35). Incorporant l'approximation (3.13) dans cette 
condition on obtient: 

ubzlw?w^ +c=M (3 ' 14) 

ou C contient toutes les contributions des frequences de Matsubara non-nulles 
de (2.35). Ces contributions sont non-singulieres et dependent faiblement de la 
temperature. On peut maintenant demontrer que le comportement critique 
correspond a celui du modele spherique [30]: 



U = A- 1 exp(c/T) 
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On a introduit un parametre, A -1 , pour assurer la convergence de l'integrale. 
Celui-ci correspond a peu pres a la plus petite dimension du systeme: la longueur 
de coherence £ . La variable c est a peu pres constante en temperature et est 
donnee par: 



La faible dependance en temperature de c a n = 1 a ete analysee dans [30]. II est 
possible de faire un lien entre la frequence caracteristique et la longueur de 
correlation. Dare et al. [31] ont trouve que dans le cas repulsif, l'approximation 
TPSC prevoit que la frequence caracteristique des fluctuations magnetiques suit le 
comportement suivant u> c oc avec z = 2 (z est appele exposant dynamique). La 
symetrie entre notre approximation et celle du cas repulsif montre que dans notre 
cas la frequence caracteristique de la susceptibilite de paires a un exposant 
critique avec la meme valeur, z = 2. Ainsi, la frequence carateristique decroit plus 
rapidement que la temperature. 

Pour un reseau de taille finie le resultat est different. Supposons que le 
reseau est un carre de taille L x L. Ici en traitant separement le vecteur q = et 
en remplagant le reste de la somme discrete sur les vecteurs d'onde par une 
integrale, la relation (3.14) donne: 



Supposant que la valeur moyenne de double occupation est la meme pour un 
reseau de taille finie que pour celui de taille infinie (ce qui est une bonne 
approximation car cette quantite est peu sujette aux effets de taille), on obtient 
la relation: 



Ce resultat nous montre que si la longueur de correlation est petite par rapport 
a la taille du systeme, £ L << L, elle est a peu pres identique a celle du reseau 
de taille infinie £ L ps ^ oq . Sinon, on aura £ L < ^ car x > ln(l +x). Dans ce 
dernier cas, la longueur de correlation perd son comportement exponentiel et suit 
plutot un comportement en loi de puissance: 



c = 7r|^ P pl^o(( n T n i) " C ) 





(3.16) 




L 



(3.17) 



T l/2 
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Ce comportement pourrait persister tant que le systeme demeure dans un regime 
classique et qu'il est alors justifie de negliger la somme sur les frequences de 
Matsubara dans la regie de somme (2.35). Cependant, on sait par la diminution 
de l'importance relative de la susceptibilite statique observee a basse temperature 
(figure 14) que le systeme redevient quantique. On montrera que le comportement 
en temperature obtenu (3.17) demontre aussi un retour vers un regime quantique 
lorsque les effets de taille sont importants (voir section 3.6.2). 



3.6.2 Effet des fluctuations sur le spectre a une particule 

Voyons comment 1' augmentation de la longueur de correlation pourrait 
induire un pseudogap dans le poids spectral. Considerons done l'expression de la 
self-energie retardee en fonction de la susceptibilite de paires, eq. (2.37). Comme 
le regime est classique, il est justifie, a basse frequence, de ne conserver que la 
contribution a la self-energie provenant de la susceptibilite statique: 

s(k„ ik n ) = u n - - uu„l £ x(q, 0)-^^—^ 



Utilisant le developpement longue portee presentee a l'equation (3.13), nous 
pouvons faire le prolongement analytique de la self-energie. 

„ n n UT f d 2 q 1 1 

^ FfU + x V )«U---^-J (2?r)2 q2 + p - + - + e( _ kF + q) _ - 



Supposant que Ton puisse lineariser la relation de dispersion, on ecrit 
e(— k F + q) ~ e(— k F ) — q • v F . Ainsi, on peut calculer de fagon analytique la partie 
reelle et la partie imaginaire de la self-energie retardee. On obtient une forme 
semblable a ce qui fut obtenu par Vilk et Tremblay dans le cas repulsif [117]: 



E' R (k F , W ) 



E'A(k F ,a;) 



UT 1 
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On est alors a meme d'evaluer le poids spectral. Celui-ci est relie a la self-energie 
retardee par la relation habituelle: 

-2E£(k, w ) 



A(k,u) = -2ImG fl (k,w) 



(u;- e (k)-S^(k, W )) 2 + (S-(k,a;)) 2 



(3.19) 



Verifions quelle est la condition necessaire pour l'apparition du pseudogap. 
La presence d'un pseudogap a basse frequence signifie que le poids spectral a 
frequence nulle doit etre convexe, c'est-a-dire que: 



d 2 A(k 



Ft 



UJ 



> 



iu=0 



Pour verifier comment on peut satisfaire une telle condition, etudions le 
comportement de la self-energie autour de la frequence nulle. Par la definition du 
niveau de Fermi [71], on a: 

e(k F )- f i + Z' R (k F ,uj = 0) = 

On peut egalement montrer que la partie imaginaire de la self-energie est 
optimale a frequence nulle (pour un liquide de Fermi on y observe un maximum 
alors que d'apres la relation (3.18), on y observe un minimum) et que la partie 
reelle y possede un point d'inflexion. Ainsi, on obtient: 



d 2 A(k,uj) 



duj 2 



uj=0 



E'«(k,0) 



u)=0 



££(k F ,0) 



dV R (k F , 



duj 



LJ = 



- S //2 



(k F! 0) 



^-8)\U pp \Te 



4tt 2 



+ 



8^ 



8C 2 ^ 



%v* F ' nv F \U pp \T£_ 

(3.20) 

ou la deuxieme ligne est obtenue en utilisant les resultats analytiques pour la 
self-energie (3.18). Regardons dans un premier temps, l'information que Ton peut 
deduire de la premiere ligne. Le premier terme a l'interieur des crochets du cote 
droit est positif tandis que le deuxieme est negatif (la partie imaginaire de la 
self-energie retardee etant negative). On reconnait que le terme entre parenthese 
est l'inverse du poids de quasiparticule (definition usuelle). Ainsi, la disparition 
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du poids de quasiparticule ne represente pas l'apparition d'un pseudogap. Au 
contraire, il indique une augmentation de la concavite du poids spectral a basse 
frequence. Ann qu'il puisse y avoir un pseudogap la partie imaginaire de la 
self-energie doit etre elevee a frequence nulle et y avoir une forte convexite. 

L'etude de la deuxieme ligne de l'equation (3.20), nous revele que lorsque 
le troisieme terme est dominant, c'est-a-dire a basse temperature et en presence 
d'une longueur de correlation petite, le poids spectral est concave. Cette concavite 
augmente avec l'abaissement de la temperature. Cependant, lorsque la longueur 
de correlation est suffisamment importante, le terme dominant est le premier 
terme. On a alors affaire a un poids spectral convexe. II y aura done un 
pseudogap a basse frequence. La condition de convexite repose principalement sur 
la comparaison entre le premier et le troisieme terme (mentionnons que le 
deuxieme terme contribue egalement a la presence d'un pseudogap): 



La condition (3.21), nous revele que plus la vitesse de Fermi sera petite 
plus il sera facile d'ouvrir un pseudogap. Ce resultat concorde avec notre 
observation que le pseudogap semble apparaitre d'abord pres des singularites de 
van Hove, (0, 7r), ou la vitesse de Fermi est plus petite. Lorsque le comportement 
de la fonction de correlation est exponentiel en temperature, tel que presente a 
la formule (3.15), il est clair qu'un pseudogap apparait a basse temperature. 
Cependant, qu'arrive-t'il lorsque les effets de taille sont importants? On a montre 
a la relation (3.17), que dans ces conditions la longueur de correlation croit en 
suivant une loi de puissance en 1/T 1 / 2 . Alors, les trois termes apparaissant dans 
la relation (3.20), ont un comportement identique en temperature, c'est-a-dire en 
1/T 1 / 2 . Ainsi, si un pseudogap est present dans le systeme, il demeurera et 
deviendra meme plus prononce avec l'abaissement de la temperature jusqu'a 
donner un gap a temperature nulle. II faut cependant, noter que notre 
approximation classique n'est surement pas valide jusqu'a cette temperature. 




(3.21) 
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3.6.3 Longueur d'onde thermique de de Broglie 

Comme on a montre precedemment, il est plus facile de creer un pseudogap 
lorsque la vitesse de Fermi est plus faible. II est done interessant d'introduire une 
nouvelle quantite, la longueur d'onde thermique de de Broglie [116]: 



En fonction de cette nouvelle quantite, la condition pour l'existence d'un 
pseudogap (3.21) prend la forme: 



II est possible d'estimer la longueur de correlation de paires et la longueur 
d'onde thermique de de Broglie a partir des mesures MCQ. Pour la longueur de 
correlation, e'est assez simple: considerant la susceptibilite de paires statique, 
x(q, iq n = 0) nous pouvons en extraire la longueur de correlation en utilisant le 
developpement autour du vecteur d'onde q = 0, equation (3.13). II est a noter 
qu'etant donne que nous avons des reseaux de petite taille, le resultat est moins 
precis a basse temperature ou la susceptibilite de paires est fortement piquee 
autour de q = 0. 

L 'evaluation de la longueur d'onde thermique de de Broglie est plus difficile. 
Nous devons d'abord faire le prolongement analytique de la fonction de Green 
pour differents vecteurs d'onde afin d'en extraire la vitesse de Fermi. L'evaluation 
de cette derniere n'est pas tres precise. Elle est extraite de la position du pic 
dans le poids spectral pour chacun des vecteurs d'onde choisis (voir par exemple 
la figure 9). La dispersion de ces pics nous permet d'evaluer les differentes 
composantes de la vitesse de Fermi. Par la suite, nous pouvons evaluer la norme 
de la vitesse de Fermi et la longueur d'onde thermique par la relation (3.22). 

Nous presentons a la figure 18 la longueur de correlation et la longueur 
thermique obtenue pour un reseau de taille 6 x 6 et de taille 8 x 8 en fonction de 




(3.22) 
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la temperature. Nous presentons les resultats pour une interaction U = —4 avec 
un remplissage n = 0.95. Les deux longueurs sont donnees en unite de pas du 
reseau. Ann de voir l'effet de taille finie nous avons extrapole la longueur de 
correlation pour un reseau de taille infinie en utilisant la relation (3.16). A 
certaines temperatures le resultat extrapole du reseau 6 x 6 ne s'accorde pas avec 
celui extrapole du reseau 8x8. Nous avons alors considere la moyenne des deux 
resultats. Le calcul de la longueur de correlation ne doit etre considere que 
comme une estimation du premier ordre. Une des raisons en est que le parametre 
A n'est pas connu. Pour notre calcul, nous avons utilise A = 7r. On peut toutefois 
verifier que le changement de la valeur de A change peu l'extrapolation au reseau 
de taille infinie. 

Aux temperatures les plus elevees, la longueur d'onde thermique est plus 
grande que la longueur de correlation et elle augmente avec l'abaissement de la 
temperature en suivant la loi de puissance 1/T. La longueur de correlation quant 
a elle augmente avec l'abaissement de la temperature en suivant un comportement 
exponentiel tel que prevu £ oc exp(c/T). Ainsi a une certaine temperature la 
longueur de correlation devient plus grande que la longueur d'onde thermique. On 
peut l'observer a une temperature legerement inferieure a la temperature ou 
apparait le pseudogap dans le poids spectral. On peut aussi observer que les 
eflets de taille finie changent peu la longueur de correlation car celle-ci est 
relativement petite par rapport a la taille du reseau et ce pour un vaste domaine 
de temperature (£ ~ 1.32 a T = 1/4.). 

En terminant, remarquons que Ton n'observe pas explicitement dans la 
figure 18 le comportement en £ oc 1/T 1 / 2 prevu a tres basse temperature. La 
raison en est simplement que les effets de taille semblent peu importants pour 
l'ensemble des points presentes. S'il existe un regime de temperature ou un tel 
comportement est present alors on s'attend a ce que la longueur thermique puisse 
rattraper la longueur de correlation (la premiere de ces longueurs ayant un 
comportement en 1/T). Ainsi, cet effet pourrait provoquer le retour vers un 
regime quantique observe a tres basse temperature pour des systemes de petite 
taille (voir section 3.4). 
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3.7 Hors demi-remplissage 

II est interessant de se demander si les resultats presenters pour le 
cas demi-rempli peuvent s'extrapoler hors du demi-remplissage. La principale 
difference vient du fait que le comportement critique est alors de type 0(2). 
Ce comportement critique permet l'existence d'un point critique de type BKT 
a temperature finie. Cependant, il est a noter que les systemes etudies par 
MCQ sont beaucoup trop petits pour nous reveler la presence d'un tel point 
critique. Ce probleme est encore plus prononce pres du demi-remplissage oil le 
comportement critique est de type 0(3) pour une large plage de temperature. 
Dans ces conditions, on peut s'attendre a ce que le crossover entre les 
comportements critiques 0(2) et 0(3) ne puisse avoir lieu car la longueur de 
correlation devrait etre du meme ordre de grandeur que la taille du systeme a 
une temperature trop elevee. Ceci peut etre verifie par la comparaison du facteur 
de structure de paires uniforme au facteur de structure de charge pris au vecteur 
d'onde Q = (jt,7t). Ce dernier est defini par: 

S p (Q) = (p(Q)p(Q)) 

ou la densite de particule est donnee par: 

p(Q)= J>t(k)</v(k + Q) 

k,cr 

La comparaison entre le facteur de structure de charge et celui de paires 
est presentee a la figure (19). Sur la partie de gauche de la figure, on peut 
voir qu'au remplissage n = 0.95, ces deux quantites sont presque identiques. La 
difference ne se fait sentir qu'a tres basse temperature ou l'incertitude sur nos 
resultats est importante. Pour fins de comparaison, nous presentons dans la meme 
figure les facteurs de structure pour le remplissage n = 0.8. On voit que dans ce 
dernier cas le comportement critique est de type 0(2) car le facteur de structure 
de charge n'est pas critique. 

Une autre quantite nous revele la brisure de la symetrie 0(3), c'est la 
valeur du potentiel chimique. Comme on a vu au chapitre 1, c'est cette quantite 
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Tableau 4: Potentiel chimique du modele de Hubbard 
a U = — 4 tel que mesure par MCQ pour un reseau 
8x8 avec une discretisation du temps imaginaire 
At = 1/10. 



n 


T 


(jt-U/2 


0.95 


1/3 


-0.082 


0.95 


1/4 


-0.06 


0.95 


1/5 


-0.06 


0.95 


1/8 


-0.07 


0.8 


1/3 


-0.275 


0.8 


1/4 


-0.26 


0.8 


1/5 


-0.25 



qui brise la symetrie 0(3) dans le modele de Hubbard. Dans le tableau 4, nous 

presentons les valeurs du potentiel chimique utilisees dans les calculs MCQ pour 

le reseau 8x8. Aux valeurs presentees, nous avons soustrait le terme U/2. 
Rappelons que pour le cas demi-rempli, on a fj, — U/2 = 0. 

On remarque qu'au remplissage n = 0.95, aux temperatures d'interet, le 
potentiel chimique est bien inferieur a la temperature. En comparant les valeurs 
presentes au tableau 4 a la figure 19, on peut remarquer que la difference entre 
le facteur de structure de paires et le facteur de structure de charge (voir figure 
19) se fait sentir lorsque la temperature est du meme ordre de grandeur que le 
potentiel chimique (T ~ 2/x): T ~ 1/8 au remplissage n = 0.95 et T ~ 1/3 a 
n = 0.8. 

Nous avons publie une etude detaillee du remplissage n = 0.95 [2]. Nous 
avons montre dans ce cas nous observons egalement la presence d'un pseudogap 
et ce a une temperature bien superieure a la temperature critique. Dans la figure 
20 nous montrons le poids spectral obtenu par l'application de la MME aux 
donnees MCQ. Pour le remplissage n = 0.95, nous montrons le poids spectral au 
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vecteur d'onde k= (0, 7r). Le fait que ce point n'est pas exactement sur la 
surface de Fermi est responsable de Fasymetrie observee entre les frequences 
negatives et les frequences positives. Dans la partie de droite de cette meme 
figure nous presentons le poids spectral au remplissage n = 0.8 pour le vecteur 
d'onde k= (0, Sir/A). On voit que ce point est tres pres de la surface de Fermi. 
Pour les deux remplissages nous presentons les resultats pour une temperature 
T = 1/5 avec une interaction U = —4. 

Nous voyons qu'a la temperature T = 1/5, pour chacun des remplissages 
etudies, le poids spectral a une particule presente un pseudogap. Ceci semble 
montrer que pres du demi-remplissage, la temperature a laquelle apparait le 
pseudogap, T* , depend faiblement du remplissage. Si Ton compare la valeur 
du poids de quasiparticule pour chacun de ces remplissages, il semble que la 
temperature T*(n) diminue legerement au fur et a mesure que Ton s'eloigne 
du demi-remplissage. En fait, ce comportement pourrait tres bien etre similaire 
au comportement de la temperature critique champ moyen, T°, en fonction 
du remplissage. Pour illustrer ce fait, nous presentons au tableau 5, la valeur 
de la temperature critique predite par un calcul champ moyen pour differents 
remplissages. Nous donnons egalement la valeur du poids spectral basse frequence 
mesure par MCQ a T = 1/4 (voir equation (3.2)). II faut dire que le poids 
spectral ne demontre pas tout a fait un comportement monotone en remplissage. 
Cependant, il faut denoter que pour le remplissage n = 0.95, le vecteur d'onde 
pour lequel nous avons evalue ce poids, k= (0, 7r), n'est pas situe sur la surface 
de Fermi et se trouve legerement en dehors de la mer de Fermi. Cela pourrait 
expliquer le fait que le poids spectral basse frequence y est legerement inferieur a 
celui mesure au demi-remplissage. Ann d'obtenir une comparaison plus detaillee il 
serait interessant d' avoir une mesure de la temperature ou apparait le pseudogap 
en fonction du remplissage, T*(n). 

Le fait que le comportement en remplissage de T* serait semblable a celui 
de T° renforcerait notre argumentation selon quoi le pseudogap est du aux 
fluctuations supraconductrices. 
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Tableau 5: Temperature critique predite par une 
theorie champ moyen pour differents remplissages pour 
un reseau de taille infinie avec une interaction 
U = —4. Pour fins de comparaison, nous presentons 
egalement la valeur du poids spectral basse frequence 
mesuree a T = 1/4 pour un reseau de taille 8x8. 



n 


- 1 c 


z 


1. 


0.7386 


0.446 ±0.001 


0.95 


0.7374 


0.441 ± 0.003 


0.8 


0.7198 


0.514 ±0.001 
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A basse temperature, lorsque les effets de correlation sont importants le 
modele de Hubbard presente des caracteristiques pouvant etre assimilees a un 
regime classique renormalise. Pres du demi-remplissage, ce modele possede un 
parametre d'ordre qui presente une symetrie plus elevee. Nous avons montre que, 
dans ces conditions, une approximation de type RPA avec un vertex renormalise 
concorde avec ce qui est observe par les simulations Monte Carlo. Ce resultat 
avait deja ete observe et utilise comme justification heuristique de l'approximation 
TPSC appliquee au modele de Hubbard repulsif. 

Ces etudes prealables nous ont fourni une base pour la generalisation de 
l'approximation TPSC. En plus d'enrichir le formalisme a la base de cette 
approximation, nous l'avons appliquee au modele de Hubbard attractif. Pour le 
cas a interaction intermediaire etudie nous avons montre que notre approximation 
permet de retrouver, sans parametre ajustable, la valeur du vertex observee par 
les simulations MCQ. 

Nous avons egalement montre que dans le regime classique renormalise le 
spectre a une particule est caracterise par un abaissement du poids a basse 
energie au niveau de Fermi: un pseudogap. Au debut de cette these, nous avons 
argumente que la presence du pseudogap dans un large domaine de temperature 
dans le modele de Hubbard attractif permet de mieux comprendre la presence 
d'un pseudogap dans les supraconducteurs a haute temperature critique et ce a 
des temperatures beaucoup plus importantes que la temperature critique. Notre 
explication provient de l'elargissement de la symetrie du parametre d'ordre qui 
serait responsable d'un abaissement de la temperature critique sans toutefois avoir 
d'effet marque sur la temperature ou apparaissent les fluctuations critiques. Notre 
etude a ete presente pour une interaction inferieure au couplage maximisant la 
temperature critique (regime d'interaction faible a intermediaire), U = —4. Ainsi le 
pseudogap observe n'est pas du a des effets d'interactions fortes. 

Certains resultats experimentaux sur les supraconducteurs a haute 
temperature critique permettent d'appliquer notre etude a ces systemes. Corson et 
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al. [25] ont observe par des mesures de conductivity a haute frequence que le 
regime ou apparait le pseudogap dans les SC hTc est un regime classique 
renormalise. De plus, leur etude a montre que la temperature critique est associee 
a une disparition de la coherence de phase. 

Malgre les resultats obtenus, il semble evident que notre approximation 
ne permet pas d'expliquer les resultats observes au dela du cas d'interaction 
intermediate. Nous avons decouvert recemment qu'il est possible, dans l'esprit 
de 1' approximation TPSC, d'ameliorer le formalisme en ajoutant un vertex 
dynamique non uniforme (voir l'annexe C). Une telle approximation pourrait 
permettre d'ameliorer l'estimation de la susceptibilite critique. Cependant, il y 
aurait un effet plus important sur l'estimation de la self-energie. De plus, il serait 
interessant de verifier si elle permettrait d'obtenir des resultats interessants a plus 
forte interaction. 

En terminant, mentionnons egalement les ameliorations possibles de la 
methode Monte Carlo. Dans un premier temps, il serait possible d'accelerer la 
vitesse de convergence par un choix approprie d'estimateurs. Assaraf et Caffarel 
[6] ont presente une methode ou, par l'utilisation de simulations prealables, il est 
possible de renormaliser les estimateurs pour de nouveaux estimateurs avec une 
plus faible variance. Une autre amelioration notable pour le cas a forte interaction 
et a basse temperature, consiste a ameliorer l'ergodicite de l'algorithme. Une telle 
amelioration a ete propose par Iwamatsu et Okabe [52]. Leur methode est basee 
sur les travaux de Tsallis et Stariolo [112] et ils l'ont appliquee a l'etude par 
simulation Monte Carlo d'un double puits quantique. 

La technique de prolongement analytique par maximisation d'entropie 
pourrait egalement etre raffinee. Tel que mentionne dans le rapport de Touchette 
et Poulin [108], nous avons developpe un nouvel algorithme base sur une 
technique de recuit simule grandement utilisee dans les problemes d'optimisation 
[43]. La force de cette methode provient de sa capacite de trouver le mimimum 
absolu d'une fonction lorsque celle-ci possede plusieurs minima locaux. Cependant, 
il est difficile d'evaluer l'incertitude sur les resultats obtenus par cette technique. 
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Une methode qui nous permettrait de faire un prolongement analytique tout en 
permettant d'evaluer l'incertitude sur le resultat obtenu serait une technique de 
Monte Carlo bayesien [42,95,96]. 
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Table de formules 



Pour faciliter l'utilisation de la technique developpee, cette annexe constitue 
une table des principales formules utiles. Aucune des formules presentees ici ne 
depend du choix de l'hamiltonien. Le lecteur interesse a un hamiltonien particulier 
est done invite a evaluer les equation du mouvement pour son hamiltonien et a y 
aj outer les formules qui lui seront utiles. 

A.l Definitions 

On considere une fonction de Green de forme matricielle 2x2 dont les 
elements hors diagonaux sont associes aux parametres d'ordre et les termes 
diagonaux sont les termes normaux. On considere le systeme en presence d'un 
champ source externe, <fi, couple aux parametres d'ordre. 

L 'equation de Dyson peut etre considered comme une definition de la 
self-energie 



(Go \1, 3) - 0(1, 3) - 17(1, 3; 0))G(3, 2; 0) = 5(1 - 2) 



Definition de la fonction de correlation: 



X (l,2,3,4;0) = 



S 



G(3,4;0) 



50(2,1) 



oil 




Definition du vertex irreductible: 



r(l,2,3,4;£) 



5 



^(3, 4; a) 
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ou 



50(1,2) 



5512(1,2) 



5021(1,2) 




On peut passer d'une expression en fonction du champs a une en 
fonction de Q par un changement de variable. Par exemple, ce changement pour 
le champ peut s'ecrire: 



*(1,2;£O = 0(1,2) 



G H d[0]=5 



A. 2 Regies exact es 

II est a noter que du a la forme matricielle de l'operateur derivee, la 

derivation d'un produit de matrices ne se fait pas de fagon usuelle. En utilisant 

le fait que cette matrice n'a que des composantes hors diagonales, on peut 
demontrer la relation suivante: 

5 f 5 \ 5 5 

Cette relation se demontre par la regie de commutation entre une matrice 
diagonale, D, et une matrice hors diagonale, H, (matrices 2x2): 

HD = a x T>a x li 

et par la regie de commutation entre deux matrices hors diagonales, H x et H 2 : 

H!H 2 = a x H 2 H 1 a x 



Pour demontrer la regie de derivation en chaine lorsque Ton passe d'une 
derivation par rapport au champ a une par rapport au champ Q, il suffit de le 
faire de fagon usuelle puis de ramener le tout sous forme matricielle. 



5 = / 5G 12 (3,4) \ 5 / 5G 21 (3,4) \ 5 

5<f>(2,l) \ 50(2,1) 7^12(3,4) V 50(2,1) ; 5a 2 i (3, 4) 



-G HD (3,4;0) ) 5 —= f =r+ ( 6 ——-a x G HB (3,4-,ct>)Y-r^=^. 



50(2,1) nuy ' ,T 7^(3,4) V 50(2,1) x auy ' 7 5^(3,4) 

(Al) 
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La definition de 1' inverse de la fonction de Green permet de faire un lien 
entre la dependance de G~ x sur <fi et celle de G. II suffit de considerer la relation: 

G(l,3;0)G- 1 (3,2;0) = <5(l-2) (A2) 

Cette relation nous permet de deriver plusieurs relations interessantes. Par 
exemple, on peut Putiliser pour demontrer la relation de Bethe-Salpeter. 



A. 3 Calcul au premier ordre 

La fonction de Green au premier ordre en Q: 

G(l, 2; Q) = G(l, 2) + 2) + 0(G 2 ) (A3) 

nous permet d'evaluer G _1 au premier ordre: 

G-\l, 2; Q) = G-\\, 2) - G -1 (l, 3)^(3, 3)G -1 (i, 2) + 0(G 2 ) (A4) 

Combinant (A. 4) avec la relation de Dyson (2.12) on obtient une expression pour 
le champ externe au premier ordre en Q: 

~5 



*(l,2;0) = -0(4,3) 



+ 0(G Z ) 
G=o (A5) 

=g(4,3)( ( 7 I G- 1 (l,4) ( 7 I G- 1 (3,2)-r(3 ) 4,l ) 2))+0(e 2 ) 



Pour ce qui est du developpement au premier ordre en </>, on obtient: 

G(l, 2; 0) = G(l, 2) - 0(4, 3) X (3, 4,1,2) + O(0 2 ) (A6) 

qui decoule de la definition de %. Pour G _1 , le resultat correspondant se 
demontre en considerant l'equation de Dyson (2.12) avec la definition du vertex 
(2.23) et la regie de derivation en chaine (A.l): 

G-^l, 2; <P) = G-^l, 2) - 0(1, 2) + 0(3, 3) X (3, 4, 5, 6)r(6, 5, 1, 2) + O(0 2 ) (A 7) 

C'est la combinaison de (A. 6) avec (A. 7) dans (A. 2) qui permet de demontrer la 
relation de Bethe-Salpeter. 
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A. 4 Calcul au deuxieme ordre 

Pour assurer une autocoherence dans 1'estimation du vertex, il faut aller au 
deuxieme ordre. Ceci correspond au premier terme de la dependance en (f> de la 
partie diagonale des matrices. L'observation que la partie hors diagonale des 
matrices peut s'exprimer comme une somme de puissances impaires de simplifie 
grandement le probleme. On peut ainsi obtenir G(G) au deuxieme ordre en Q 
sans introduire de nouvelles fonctions: 

G(l,2;0) =G(1,2) + 0(1,2) 

- g(E, 5)(a x G-\4, BKG" 1 ^, 3) - r(5, 6,4, 3)) XhD (3 A, 1, 2; Q) + 0(G 3 ) 

(A.8) 

Pour G _1 , on utilise l'equation de Dyson (2.12), les regies de la section 
(A. 2) et le fait que la matrice &(G) etant hors diagonale n'a dans la phase 
normale que des puissances impaires en Q. 

Gr\\, 2; Q) =G" 1 (1, 2) - G^l, 3)0(3, ^G" 1 ^, 2) 

(A9) 

-0(4,3)r HD (3,4,l,2;£) + O(£ 3 ) 



Afin de bien obtenir les relations (A. 8) et (A. 9) a l'ordre quadratique en 
Q, il nous faut obtenir r HD et % HD a l'ordre lineaire en Q. Ceci peut se faire en 
utilisant les definitions (2.23) et (2.14) et que seule les parties diagonales de G et 
de U ont des termes de puissances paires dans leur developpement en Q. On 
obtient pour P HD une expression ou apparait un vertex a trois particules: 

r HD (l,2,3,4;£) = 0(5,6) 



50(5,6) 50(2,1 



-S(3, 4; 0) 



O(0 3 ) 



(A10) 



De la meme facon l'expression pour % HD fait apparaitre une fonction de 
correlation a trois particules: 

"5 5 



X HD (l,2,3,4;0) = -#(5,6;0) 



.(50(5,6) 50(2,1) 
= ^(8,7)(- ( 7 :c G- 1 (5,8) ( 7 :i; G- 1 (7,6) + r(7,8,5,6)) 



G(3,4;< 
5 



50(5,6) 50(2,1) 



G(3,4; 



J 4>=o 
(All) 
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ou on a fait appel au developpement de ^ au premier ordre en Q, (A. 5). 
Introduisant les developpement du deuxieme ordre (A. 8) et (A. 9) dans (A. 2), on 
peut montrer qu'il existe une relation liant le vertex a trois particules et la 
fonction de correlation a trois particules: 



5£(6,5) 50(4,3) 



= - a x G- l (3, Q)a x G~ 1 (l, 4)G'\5, 2) + G _1 (l, 7) 



x a x ^G-^g, 4KG- 1 (3, 10) - T(3, 4, 9, 10) ) a a 
x [rr,.G : {TT.(>)a,.G '(o. " : r(5Ji.TT.12) 
"5 5 



50(11,12) 5^(9,10) 



G(7,8;<£) 



G- x (8,2) 



6=0 



(A12) 



En annexe C, on utilise ces resultats afin de proposer une fagon d'ajouter 
la coherence entre les proprietes a une et a deux particules. 
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Transformation de Legendre 



B.l Definition 



Avant de voir l'application de la transformation de Legendre, voyons la 
definition de cette transformation. La description presentee ici s'inspire de celle 
retrouvee dans Encyclopcedia of Mathematics [38]. 

Supposons que Ton connaisse une fonction / : A — > M reguliere sur 
l'ensemble ouvert A d'un espace norme X. On peut considerer cette fonction 
comme une surface dans l'espace. La surface etant alors decrite par l'ensemble des 
points (x,f(x)) oil x € A 

La possibilite d'evaluer la transformation de Legendre, /*, de la fonction / 
repose sur l'existence d'une seconde description de la surface. Celle-ci correspond 
a la famille des plans tangents a la surface. On doit alors faire appel a une 
nouvelle variable servant a parametriser les plans tangents, x* G B et a une 
fonction tangente affine pour la description de ces plans, f*(x*). L'ensemble 
ouvert B appartient a un espace norme X* . L'equivalence entre les deux 
descriptions provient de l'existence d'une application biunivoque de A sur B. 
Cette application permet de definir la transformation de Legendre. Supposant que 
/ est une fonction sur R n et que le determinant det(<9 2 / jdxfix^) est finie sur A, 
la transformation de Legendre sera donnee par les formules suivantes: 



Ainsi, lorsqu'on s'interesse a la transformation de Legendre d'une fonction, 
il faut connaitre la derivee de Frechet correspondante fix) et la definition de la 



f(x) = x*, f*(x*) = (x,f(x))-f(x) 



oil 
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norme dans l'espace d'interet. La transformation de Legendre introduit done une 
nouvelle variable prenant des valeurs dans un ensemble B qui peut etre relie de 
fagon biunivoque a l'ensemble A des valeurs admises par la variable x. 



B.2 Transformation de Legendre de la fonction de partition 

Dans notre cas, nous nous interessons a la transformation de Legendre de 
la fonction de partition Z[<f>]. Par rapport a la description de la section precedente 
nous pouvons identifier / au logarithme de la fonction de partition lnZ[0] et la 
variable x au champ 0. L'operation < x,y > est quant a elle remplacee par la 
trace Tr{xy} telle qu'introduite en mecanique statistique quantique. La derivee de 
Frechet de la fonction de partition correspond a la partie hors diagonale de la 
fonction de Green G HD (</>). 

Nous devons done trouver une transformation qui permet d'appliquer 
l'ensemble ouvert des valeurs admissibles par la variable de fagon biunivoque 
sur un ensemble ouvert des valeurs d'une variable Q. Utilisant la definition 
presentee a la section precedente nous pouvons introduire la transformation de 
Legendre que nous avons presentee a la section 2.2: 



W[G] = (-lnZ[$] - Tr(G HD (l, 2; 0)0(2, 1))) 



ou 

Ghd(1 ' 2 '^ ) = -^i) 



Tout comme il nous a ete possible de montrer que la partie hors diagonale 
de la fonction de Green pouvait etre obtenue par la variation de la fonction de 
partition par rapport au champ externe cf), eq.(2.7), la partie hors diagonale de la 
self-energie pourra etre generee par la variation d'une fonctionnelle par rapport 
au champ introduit par la transformee de Legendre, Q. Nous avons montre (2.22) 
que: 



8W[Q] =-•(!, 2; 0) 



dg(2, i) 
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Or, de 1' equation de Dyson, (2.12), nous avons: 

2; Q) = -G-^(l, 2; Q) - I7 HD (1, 2; 0) 



Si Ton considere l'approximation du premier ordre en terme anormal, on 
obtient la relation suivante: 

G H ^(l,2;a) = -G- 1 (l,3)C?(3,4)G- 1 (4,2) 

= " ^^ I y^(^(3,4)G- 1 (4,5)^(5,6)G- 1 (6,3)) 



Des lors, on peut montrer qu'au premier ordre en terme anormal les 
termes hors diagonaux de la self-energie seront donnes par: 



Wl, 2; Q) =|^i) ( W[G] + l T < Q ^ 4 ) G_1 ( 4 ' 5 )^( 5 ' 6 ) G_1 ( 6 ' 3))) 

6 -W 



'50(2,1) 



Ainsi les termes diagonaux du vertex irreductible sont donnes par la 
derivation seconde de cette meme fonctionnelle par rapport aux termes anormaux 
de la fonction de Green. 



r(l,2,3,4) = 



50(2,1)50(4,3) 



e[G] 



g=o 
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Autocoherence 

C.l Comparaison a un hamiltonien avec interaction non-locale 

Une critique importante de notre approche vient du fait que celle-ci ne 
presente pas d'autocoherence entre les proprietes a une particule et les proprietes 
a deux particules. Ainsi, la self-energie au deuxieme niveau d'approximation est 
differente de celle au premier niveau. Ann d'ajouter de 1' autocoherence dans notre 
approximation il faudrait pouvoir utiliser la self-energie obtenue au deuxieme 
niveau pour reevaluer le vertex et la fonction de correlation qui pourrait alors 
etre utilises pour passer au niveau suivant et ainsi de suite, iterer jusqu'a 
convergence de l'approximation. 

Dans cette annexe, nous proposons une methode pour ajouter 
l'autocoherence. Le fait qu'au deuxieme niveau la self-energie est dynamique et 
delocalise peut sembler compliquer les choses si Ton veut garder une certaine 
coherence. Cependant, avant de debuter les calculs, on peut se guider en etudiant 
des relations simples mais appliquees a un hamiltonien un peu plus complexe 
que celui de Hubbard. Nous allons considerer un hamiltonien avec interaction 
non-locale. Ce hamiltonien presente deux termes. Le premier etant le terme 
d'energie cinetique pour les electrons. Le deuxieme pourrait etre considere comme 
un terme d'energie cinetique pour les paires de Cooper. 

H = -^(T)*(T,2>(2) + At(l)V(T,2)A(2) 

Pour un hamiltonien de ce type, la self-energie au premier ordre en 
interaction sera donnee par: 

17(1, 2; G) « V(l, 3)0(3, 3)5(1 - 2) + a x V(l, 2)G(1, 2)a x + 0{Q 2 ) 
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ou la matrice d'interaction est definie par: 

V(l,2) 



V(1 ' 2) ~» -^(2,1) 



II semble done possible d'approximer la renormalisation du vertex de fagon 
telle que la partie diagonale de la self-energie soit non-locale tandis que la partie 
hors diagonale demeure locale. Une telle approximation donnerait un vertex dont 
Pelement (1,1) pourrait etre approxime par: 

r(l, 2, 3, 4) « a x V(3, l)aj(l - 2)5(3 - 4) 

Utilisant une telle approximation avec l'equation de Bethe-Salpeter (2.24), on 

obtient une fonction de correlation de paires de forme RPA mais avec un vertex 
dynamique et delocalise au denominateur. 

Xo(o) 



x(q) = 



^(.q)Xo(q) 



C.2 Equation autocoherente pour le vertex irreductible 

Voyons maintenant comment il serait possible d'ajouter Pautocoherence. On 
considere que Ton a obtenu notre premiere approximation pour le vertex et la 
fonciton de correlation: et donne respectivement par les relations (2.30), 

et (2.33). Ainsi, utilisant l'equation du mouvement et l'equation de Bethe-Salpeter 
nous avons obtenu et G^ 2 ) donne par (2.37) et (2.38). On voit de par sa 

definition que pour reevaluer JT, il nous faudrait connaitre au premier ordre 

en Q. Pour y arriver nous pouvons proceder de la meme fagon que pour le terme 
independant de Q: en utilisant l'equation du mouvement (2.10). 

17(1, 2; G) =UQ(1, 1)6(1 - 2) - U X (1, 1, 1 ± , 3)G _1 (3, 4)^(4, 5)G _1 (5, 2) 

(C.l) 

+ U X (l,l,l ± ,3;G)G- 1 (3,2) + 0(g 2 ) 

Pour estimer le terme x(0) a l'ordre lineaire en Q, nous utilisons les 
developpements du deuxieme ordre pour G(Q) et G _1 (Q) que nous avons 
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presentes en annexe A. Utilisant le fait que les termes du deuxieme ordre doivent 
s'annuler du cote gauche de l'identite suivante: 

G- 1 (l,3;^)G(3,2;^) = ( 5(l-2) 

on peut montrer que: 

2, 3, 4; Q) = X (1, 2, 3, 4) - Q(l, 4)G(3, 2) 

+ X (l, 2, 5, 6) (r(6, 5, 7, 8)0(8, 4) + T HD (6, 5, 7, 8; 0)G(8, 4)) G(3, 7) 
+ O(0 2 ) 

(C.2) 

Incorporant (C.2) dans (C.l), nous obtenons: 
2; H D(l,2;^) = f/ ( 5(l-2)^(l,l) + f7x(l,l,5,6)r HD (6,5,7,2;e)G(l,7) + 0(C; 3 ) (C.3) 



La definition du vertex irreductible (2.23), nous permet en combinant (C.3) 
et (A. 10) d'obtenir une relation determinant la partie diagonale du vertex en 
l'absence de champ externe en fonction d'un vertex a trois particules: 



r(l, 2, 3, 4) =175(1 - 2)5(1 - 3)5(1 - 4) 
+ tfG(3,7K*(3,3,5,6)<7 a 



50(2,1)50(5,6) 



^7^(7, 4, Q) 



(C7.4) 



Notre schema d'approximation pourrait done etre developpe de la facon 
suivante: on approxime la self-energie au deuxieme ordre en Q. Partant de cette 
approximation on evahie tel que presente au chapitre 2, e t 
On y ajoute le calcul de la derivee seconde de la self-energie. Etant donne que 
cette derniere correspond a un vertex a trois particules, on pourrait garder 
cette quantite fixee et assurer l'autocoherence entre les proprietes a une et a 
deux particules en utilisant les relations de Dyson, eq.(2.12), notre expression 
exacte pour la self-energie, eq.(2.37), l'equation de Bethe-Salpeter, eq.(2.24), et 
l'expression du vertex irreductible en fonction de vertex plus complexe, eq.(C.4). 
Cette fagon de proceder s'apparente a la methode hierarchique proposee par 
Kadanoff et Martin [55] qui consiste a utiliser les equations du mouvement pour 
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les fonctions de correlation a une et a deux particules qui permettent de faire le 
lien entre elles et avec la fonction de correlation a trois particules. Par la suite, 
utilisant une approximation pour la fonction de correlation a trois particules sous 
forme d'une combinaison de fonctions a une et a deux particules, ils proposent 
une solution autocoherente pour les fonctions de Green a une et deux particules. 

C.3 Suggestion pour l'ajout de l'autocoherence 

Dans cette derniere section, nous proposons une approximation pour le 
vertex a trois particules. II est possible de developper une approximation qui 
respecte l'esprit de la technique presentee au chapitre deux. Nous utilisons la 
relation (A. 12). Encore une fois, on introduit une matrice, A, qui sera determinee 
par les correlations locales: 



Pour justifier cette approximation, on peut considerer l'equation du 
mouvement pour la fonction de correlation a deux particules: 



On voit apparaitre dans cette relation une fonction de Green a trois particules 
via le terme de derivation. On propose done une approximation pour cette 
fonction qui permet de resoudre l'equation du mouvement a partir des quantites 
introduites precedemment. 



C/G(l,7)x(3,4,9,T0) - 




J <p=o 



(C.5) 



a x Go\l, 5K*(3, 4, 5, 2) = - 5(1 - 4)G(3, 2) + UG(1, 2) X (3, 4, 1, 1) 




U 



5 6 

(50(4,3) 50(M) 



G(l^,2;0) 



J <j>=0 



a x U(l, 5)a,x(3, 4, 5, 2) - U X (S, 4, 1, 1)G(1, 2) 



^G(1,4KG(3,1)G(1,2) 



(C.6) 



96 



Annexe C 



On peut verifier qu'au premier ordre en U l'egalite entre les deux cote est 
verififee de fagon exacte. Ainsi, en introduisant 1' approximation (C.6) dans la 
premiere ligne de (C.5) on obtient la deuxieme ligne de cette meme equation ou 
on a ajoute le facteur multiplicatif A. Combinant les relations (C.4) et (C.5), il 
devient alors possible d'approximer le vertex irreductible. 

r(l,2,3,4) « 175(1 - 2)5(1 - 3)5(1-4) - UAS(3 - 4)a x G(3, 5)a x G(6, 3)r(l, 2, 5, 6) 

(C.7) 

La relation peut se simplifier si Ton considere en premiere approximation 
le resultat obtenu au chapitre deux et que Ton utilise la relation (C.7) afin 
d'ameliorer l'approximation. 

r (2) (l,2,3,4) = C/5(l-2)5(l-3)5(l-4)-;7A5(l-2)5(3-4)a :E G (1) (3,l)cT :E G (1) (l,3)r (1) 



Supposant que les iterations permettent de converger vers une 
solution unique, on obtiendrait alors un vertex a deux points, 
r(l, 2, 3, 4) = T(l, 2, 3, 4)5(1 - 2)5(3 - 4) = r*(l, 3)5(1 - 2)5(3 - 4), correspondant a 
la solution de la relation: 

T*(l,2) = US(1 - 2) - t^il(7 a .G*(2,3)<7 a .G*(3,2)r*(l,3) 



On voit que Ton obtient une forme RPA pour la renormalisation du vertex. 
Etant donne que le coefficient devant le vertex droit est negatif dans l'espace de 
Fourier, on remarque que le vertex sera inferieur, en valeur absolue, a l'interaction 
nue. De plus, il est facile de montrer que le comportement aux grandes frequences 
de ce vertex est identique au resultat exact: lim _ >oo r(q, iq n ) — > U. II ne reste 
ainsi qu'a determiner le facteur multiplicatif A. 



Considerant la relation ou nous avons introduit le facteur A. Nous pouvons 
la reecrire sous la forme: 



A (x(3, 4, 1, 1) - Xo (3, 4, 1, 1)) G(l, 2) 



^(1,5) 



U 

5 



(7^(3,4,5,2) 
5 



(C.8) 
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ou 



Xo (3, 4, 1,1) = -^0(1,4)^0(3,1) 

Le deuxieme terme du cote droit de l'egalite est assez simple a evaluer de fagon 
locale. Dans un premier temps, on peut ecrire ce terme sous la forme suivante: 



_ 50(4, 3) 50(1,1 



-0(1* 2; 0) 



= X (3, 4, 1,1)0(1,2) 



i=0 



(V T (3)^(4)^(1)^}(1)^(1-)V}(2)) 



(^(3)^}(4)^ T (l)^(l)^(l+)^(2)) 



Le premier terme du cote droit de l'equation (C.8) est plus complique. De 
fagon locale, on peut trouver une expression simple en utilisant le fait que la 
self-energie peut s'ecrire: 



17(1, 2) = -UG-\1, 3K % (2, 2, 3, 2±)a a 



Comme presente au chapitre 2 lors du calcul de la self-energie, equation 
(2.25) et (2.28), nous regardons la contribution provenant des differents termes 
locaux pour determiner A de l'equation (C.8). II faut faire attention car plusieurs 
ordonnances sont possibles: 3 = 1*, 2 = 1 ± , 3 = 2 ± et 4 = 1 ± . II y a done 12 
termes differents: 6 pour l'ordonnance respective des indices 1, 2 et 3 multiplie 
par un facteur deux provenant de l'ordonnance des indices 1 et 4. Des lors, on 
obtient l'approximation: 



2 . 



jr . _ _l/n 2 -4(n T n i + |(n 2 -4n(n T n i )) (2 + 3n)(n T ^)-5n 2 /4 

a«w - -) ^ - + 2^ n 

+ T*—^ G ^ 6) {—CT- + — —) [n ~ —n ) 

ou I est la matrice identite. En fait, il serait sans doute plus juste de presenter 
ce resultat comme etant un ansatz plutot qu'une approximation. 
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Prolongement analytique 



Comme nous l'avons vu au dernier chapitre, il est possible de determiner le 
poids spectral a partir de la connaissance de la fonction de Green en temps 
imaginaire. On peut etablir un lien entre ces deux quantites en utilisant 
la representation de Lehmann [110]. La relation ainsi obtenue est appelee 
representation spectrale de la fonction de Green en temps imaginaire. Celle-ci est 
donnee par la relation (D.l). 



Son domaine de validite est restreint a rG (0,/3) (L'utilisation de "(" au lieu de 
"[" signifie que la borne n'est pas incluse dans le domaine). 

Done, nous voulons determiner le poids spectral, A(uj), pour un ensemble 
de valeurs de uj en utilisant la connaissance de la fonction de Green, G(t) pour 
un ensemble de valeurs de r. II est a noter que pour des donnees Monte Carlo, 
notre connaissance de la fonction de Green n'est que statistique et qu'a chaque 
valeur de G(r) nous devons associer une incertitude, c(r). 

Le probleme auquel nous faisons face est que 1' equation (D.l) n'est 
pas reversible. Plus precisement on qualifiera ce probleme comme etant "mal 
pose". Dans la section qui suit, nous detaillons les differentes raisons pour 
laquelle ce probleme est dit mal pose. Nous mentionnons comment chacun de 
ces problemes pourrait etre regie. Dans la seconde section nous decrivons la 
fonctionnelle d'entropie introduite afin de regulariser le probleme. Enfin, dans la 
derniere section, nous proposons differents algorithmes permettant de determiner 
le poids spectral. Une comparaison de deux de ces algorithmes a ete presente par 
Touchette et Poulin [108]. 
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D.l Problemes mathematiquement mal poses 

Baym et Mermin [9] ont montre que si Ton connait la fonction de Green 
en temps imaginaire sur toute la plage (0, (3) alors il n'existe qu'une seule 
fonction A(u) permettant de satisfaire la relation (D.l) et ayant les proprietes 
d'un poids spectral [110]. Pour le cas qui nous interesse, on peut dire que nous 
n'avons pas suffisamment d' information pour determiner sans equivoque le poids 
spectral. Ce manque d'information provient de differentes raisons: premierement, la 
fonction de Green n'a ete evaluee qu'en un certain nombre de points en temps 
imaginaire, r^i G [1, N]); deuxiemement, nos donnees sont caracterisees par une 
incertitude statistique. 

Tel que mentionne par Schafer et Sternin [99], il existe plusieurs fagons de 
regulariser le probleme. D'abord, la discretisation des frequences permet de reduire 
le nombre de solutions possibles. Cependant, ga n'aurait pas de sens de reduire le 
nombre de frequences a un point tel que la solution soit unique car la solution 
n'aurait probablement pas une allure physiquement acceptable. De plus, elle 
ne respecterait peut etre pas certaines proprietes des poids spectraux telle la 
positivite. La discretisation des frequences ramene le probleme a la solution de la 
relation: 

M 

G(t<) = ^K^u^A^), ie[l,N] (D.2) 

3=1 

ou on a introduit le noyau 

„ — TLO 

K <r>") = TTT& {D - 3) 

Sous forme matricielle, la relation (D.2) prend la forme simple G = KA. 
Ainsi la solution du probleme devient A = K~ 1 G. Cependant, comme nous avons 
mentionne, lors de la discretisation on doit garder un nombre de frequences plus 
importants que le nombre de donnees en temps imaginaire, M > N. Ainsi, la 
matrice K est une matrice rectangulaire et son inverse n'est pas unique. 

Ann de determiner une solution unique, on doit faire appel a d'autres 
methodes de regularisation comme la maximisation d'entropie. Cette methode 
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consiste a choisir parmi les solutions celle qui contient le moins d'information 
possible. Par information contenue dans une fonction, on parle du comportement 
de cette fonction. Plus une fonction est monotone, moins on peut en tirer 
d'information, le cas limite etant la fonction constante. Pour quantifier cette 
propriete on utilise la fonctionnelle d'entropie, S(A). Celle-ci mesure l'absence 
d'information contenue dans une distribution. 

Le choix de la methode de maximisation d'entropie repose sur le fait 
qu'etant donne que le poids spectral est evalue en un nombre de points plus 
eleve, celui-ci a la capacite de contenir plus d'information que ne peut en 
contenir la fonction de Green. Ce qui n'a pas de sens car on utilise la fonction de 
Green afin de determiner le poids spectral, ainsi l'information contenue dans le 
poids spectral devrait etre limitee par l'information contenue dans la fonction de 
Green et l'information que Ton a sur les proprietes des poids spectraux. Ainsi, on 
cherchera parmi les poids spectraux qui satisfont la relation (D.2) celui qui 
maximise la fonctionnelle d'entropie S(A). 

II est a noter que l'imposition de certaines conditions physiques pour 
le poids spectral peut reduire de fagon considerable le nombre de solutions 
possibles. Par exemple, on peut montrer que le fait que le poids spectral soit 
positif en tout point, impose que la fonction de Green ait un maximum unique 
dans l'intervalle (0,/3) et un comportement monotone de part et d'autre de ce 
maximum. La presence de bruit dans nos donnees MCQ peut faire en sorte que la 
fonction de Green ne respecte pas cette propriete (Ce probleme sera plus accentue 
pour des valeurs de At plus faibles). II faut alors ajouter parmi P ensemble des 
solutions admises, celles qui ne permettent pas de satisfaire de fagon exacte la 
relation (D.2). Une fagon d'y parvenir consiste a faire une decomposition en 
valeurs singulieres du noyau, eq. (D.3), et de tronquer les valeurs propres les plus 
faibles. On obtient alors un noyau qui permet une plus grande independance du 
poids spectral par rapport a de petites fluctuations de la fonction de Green. 

La decomposition en valeurs singulieres du noyau n'est pas l'unique 
methode permettant d'elargir le nombre de solutions. On peut egalement elargir 
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l'ensemble des solutions a toutes les fonctions possibles respectant certaines 
conditions physiques et a chercher parmi ces solutions, celle qui permet de 
minimiser 



On fait alors face au probleme de minimisation des moindres carres tout 
en maximisant la fonctionnelle d'entropie. En faisant appel aux statistiques 
bayesiennes [95], on peut ramener la determination du poids spectral a la 
recherche de la fonction permettant de maximiser la probability 

P(A\G, a) = e -ix 2 (A,G)+ aS (A) (£) 4) 

Le facteur a a ete introduit afin de donner plus ou moins d'importance au terme 
d'entropie par rapport au terme de moindres carrees. Le choix d'une valeur pour 
ce facteur est assez complexe et nous referons le lecteur a Particle de Jarrell et 
Gubernatis [54]. 



D.2 Fonctionnelle d'entropie 

Pour le choix de la fonctionnelle d'entropie plusieurs auteurs proposent la 
fonctionnelle d'entropie introduite par Jaynes (voir [51]) qui est mesuree par 
rapport a un modele par defaut, m(uj). 

S(A\m) = [ dw(A(u) - m{uo) - A{u)\n(A{u) / m{u))) 



Le choix de cette fonctionnelle vient de l'entropie conditionnelle d'une 
distribution continue. Etant donne que le poids spectral n'a pas une normalisation 
assuree, on ajoute a cette entropie conditionnelle le terme A — m qui fait en sorte 
que la normalisation de A devrait etre aussi rapprochee que possible de celle du 
modele par defaut. On peut verifier que notre entropie a un maximum unique 
donnee par A = m et qu'elle respecte la propriete de convexite de l'entropie. Pour 
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le demontrer, on suppose que le poids spectral est une quantite positive en tout 
point. II est a noter que cette entropie donne l'entropie de la distribution A 
mesuree par rapport a celle de m et sera done negative pour toute fonction A et 
nulle pour A = m. Le modele par defaut doit contenir le moins d'information 
possible, e'est pourquoi on choisit m(u)) = m. On peut verifier qu'en general, le 
terme en \ 2 dans (D.4) est dominant pour la maximisation de la probability 
(D.4) et qu'un changement du modele par defaut ou de l'entropie a peu d'effet 
sur le poids spectral obtenu. 

D.3 Algorithme pour la maximisation d'entropie 

Parmi les algorithmes introduits pour trouver le minimum de l'equation 
(D.4), on mentionne generallement les algorithmes de Bryan [17] et de Meshkov 
et Berkov [78]. Ces deux algorithmes sont deux methodes differentes pour la 
recherche du maximum de l'equation (D.4). II est a noter que ces deux methodes 
donnent des resultats quelques peu differents [108]. Pour l'etude du poids spectral 
dans cette these notre choix s'est porte sur l'algorithme de Bryan, car cet 
algorithme s'est avere plus rapide que celui de Meshkov. On a toutefois verifie 
que les resultats obtenus par chacun des algorithmes ne sont pas tres differents 
de fagon quantitative aussi bien que de fagon qualitative. Nous avons egalement 
developpe un nouvel algorithme donnant des resultats sensiblement identiques a la 
methode de Bryan. L'avantage de ce nouvel algorithme permet de separer les 
differents problemes relies a la technique de prolongement analytique et ainsi de 
faire ressortir 1'effet de la precision des donnees. Malgre que celui-ci ne fut pas 
utilise pour les resultats presentes au troisieme chapitre, nous tenons a presenter 
ce nouvel algorithme car nous pouvons ainsi detailler les differents elements des 
techniques de regularisation. 

S'inspirant des travaux de Geman et Geman [43], le nouvel algorithme vise 
a chercher parmi les solutions de l'equation (D.2), celle qui maximise l'entropie. 
La difficulte dans la recherche du maximum de la probability (D.4), vient de la 
presence de minimum locaux dans la partie provenant des moindres carres, x 2 ■ 



103 



Annexe D 



Cependant, le terme d'entropie est parfaitement convexe. C'est pourquoi notre 
algorithme consiste a considerer comme point de depart une valeur de a tres 
elevee pour ainsi n'accorder que peu d'importance au x 2 . Ce terme ne devrait 
induire qu'un petit deplacement dans la position du maximum. Trouvant la 
position de ce maximum, on change notre modele par defaut, on abaisse 
legerement la valeur de a et on cherche la nouvelle position du maximum 
de la probability. Ainsi, on peut introduire peu a peu dans le poids spectral 
l'information contenue dans la fonction de Green. Interpretant le parametre a 
comme une temperature que Ton abaisse, cette approche correspond a une 
methode de recuit simule. Une caracteristique importante pour l'utilisation de 
cette methode vient du fait que la fonctionnelle d'entropie a un comportement 
convexe en fonction du poids spectral [7]. 

Voici plus en detail le fonctionnement de la methode. Initialement, on 
considere une valeur de a tres elevee et on cherche la fonction A maximisant 

P (A\G, a) = exp(-i x 2 (A, G) + a S(A, m)) 

Le maximum de l'entropie est a A = m. Le terme a etant important, la 
solution au probleme de maximisation devrait etre tout pres du point A = m 
mais legerement deplacee du a la presence du terme de moindres carres, x 2 - 
Appelons cette solution A 1 . On peut dire que le x 2 ajoute dans le poids spectral 
une partie de l'information contenue dans la fonction de Green. On abaisse alors 
la valeur de a et on cherche alors la fonction A qui maximise 

Pi(A\G, a) = exp(-i x 2 (A, G) + ai S(A, A,)) 

On a ainsi encore une entropie convexe mais avec un minimum deplace a 
A = A-y. Le choix d'abaisser la valeur de a vient du fait que le nouveau modele 
par defaut, A x , contient plus d'information que n'en contenait m. Ann d'aller 
chercher encore une certaine quantite d'information, on se doit done d'abaisser la 
valeur de a. On itere ainsi jusqu'a convergence de la fonction A vers une solution 
A n — A n _ x . 
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Cette methode semble donner des resultats concluants. Cependant, il reste a 
verifier que la solution est unique et n'a pas de comportement singulier par 
rapport au choix du point de depart. II reste aussi a determiner revolution 
optimale pour les valeurs du parametre a. Typiquement, nous utilisons comme 
point de depart a = ^x 2 { m jG), puis a la premiere iteration a l = ^x 2 {A 1 ^G) et 
ainsi de suite. 

Le probleme de l'incertitude statistique est tenu en compte par le fait que 
nous utilisons une decomposition en valeurs singulieres du noyau et que nous 
tronquons les valeurs propres les plus petites. Pour decire cette procedure nous 
pouvons utiliser la notation matricielle oil nous avons N valeurs de r et supposer 
que notre discretisation des frequences conserve M > N points en frequence. 
Ainsi, A est un vecteur de dimension M, G est un vecteur de dimension N 
et le noyau, K, est une matrice de dimension N x M. On procede par la 
diagonalisation en valeurs singulieres du noyau: 

K = Q 1 T,Q 2 

La matrice E est une matrice rectangulaire n x M. Seules les composantes 
le long de la diagonale donnees par les points (1,1) a (N,N) sont non-nulles. Les 
dernieres composantes etant petites, elles representent une forte dependance du 
poids spectral sur de petites variations de la fonction de Green. C'est pourquoi il 
est approprie de tronquer la matrice S en remplagant ces derniers elements 
par des zeros. On note la matrice tronquee par E s et le nouveau noyau par 
K s = Q 1 T, S Q 2 - Cette modification du noyau augmente le nombre de poids 
spectraux permettant de satisfaire la representation spectrale (D.2). 
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Figure 1: Diagramme de phase des oxydes de 
cuivre avec phase supraconductrice (SC) et phase 
antiferromagnetique (AFM). La temperature T* 
indique la temperature a partir de laquelle on mesure 
la presence d'un pseudogap (PG) au niveau de Fermi. 
La ligne au-dessus de T* identifie une temperature 
a partir de laquelle certains oxydes demontrent des 
comportements differents de ce qui est prevu par la 
theorie des liquides de Fermi. 



Figure 2: Diagramme de phase en temperature et en 
remplissage du modele de Hubbard 2D (x=l-n). La 
ligne epaisse represente la temperature critique T BKT , 
la ligne mince la temperature critique champ moyen, 
T°, la ligne hachuree, la temperature a laquelle 
apparaitrait le pseudogap et la ligne pointillee, la 
temperature oil Ton observerait le crossover dans le 
comportement critique. Nous presentons a la figure 4 
la forme schematique du poids spectral aux quatre 
points a, b, c et d. 



Figure 3: Diagramme de phase schematique en 
temperature et en interaction du modele de Hubbard 
attractif 2D. La ligne pleine represente la temperature 
critique BKT et la ligne hachuree, la temperature 
critique champ moyen. 
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Figure 4: Poids spectral pour quatre points du 
diagramme de phase du modele de Hubbard attractif 
a deux dimensions. T* represente la temperature 
a laquelle apparaitrait le pseudogap et T BKT) la 
temperature du point critique BKT. 



Figure 5: Susceptibilite de paires en fonction du 
vecteur d'onde. Les carres representent les resultats 
MCQ tandis que les croix representent les resultats 
d'une approximation RPA ou la valeur de T est 
fixee par la valeur a q=(0, 0). Les vecteurs d'onde 
sont numerates selon les points du reseau. Ainsi 
(0,4) = (0,tt) 



Figure 6: Susceptibilite de paires en fonction de la 

frequence de Matsubara. Les carres representent les 

resultats MCQ tandis que les croix representent les 

resultats d'une approximation RPA ou la valeur de 

T est fixee par la valeur a q n = 0. Les chiffres sur 

l'abscisse correspondent a la valeur de n dans la 
frequence de Matsubara q n = 2-KnT. 
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Figure 7: La partie du bas presente le poids de 
quasiparticules tel que donne par la formule (3.2) et 
la partie du haut le facteur de structure de paires. 
Ces deux quantites sont presentees en fonction de 
la temperature et ont ete mesurees par simulation 
MCQ. Le poids de quasiparticule est donne au vecteur 
d'onde k= (0,7r) et le facteur de structure au 
vecteur q=(0, 0). Nous presentons les resultats pour 
differentes tailles de reseau avec interaction U = — 4 et 
pres du demi-remplissage, n = 0.95. 



Figure 8: Poids spectral obtenu par MME pour 
differents vecteurs d'onde d'un reseau 8 x 8 en fonction 
de la frequence. Conditions usuelles U = —4, n = 1 a 
une temperature T = 1/2. La partie de gauche presente 
les donnees MCQ obtenue avec une discretisation 
At = 1/10. La partie de droite provient du resultat de 
l'approche analytique (TPSC). On a inclu dans ce 
dernier un bruit gaussien du meme ordre de grandeur 
que celui des donnees MCQ. 



Figure 9: Poids spectral obtenu par MME pour 
differents vecteurs d'onde d'un reseau 8x8 a T = 1/4 
en fonction de la frequence. Les resultats presentes ont 
ete mesures a interaction U = — 4 pour un systeme 
demi-rempli n = 1. Comme la figure precedente, la 
partie de gauche montre les resultats provenant des 
donnees MCQ tandis que celle de droite provient des 
resultats de l'approche analytique. 
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Figure 10: Densite d'etats en fonction de la frequence 
pour trois temperatures differentes: T = 1/2, 1/3 et 
1/4. Nous avons utilise les conditions suivantes:?7 = — 4, 
n = 1., At = 0.1 pour un reseau 8x8. La partie 
de gauche provient du prolongement analytique des 
donnees MCQ et la partie de cedes de la methode 
analytique 



Figure 11: Effet de la precision des donnees MCQ 
sur le poids spectral obtenue par la methode de 
maximisation d'entropie. Nous considerons la fonction 
de Green mesuree par MCQ au vecteur (0,7r) avec les 
conditions habituelles: U = —4, n = 1, At = 0.1 sur un 
reseau de taille 8x8 a temperature T = 1/4. La partie 
de gauche montre la fonction de Green considered 
pour le prolongement. La courbe du bas montre les 
donnees obtenue et pour les deux autres nous avons 
ajoute un bruit gaussien (d'amplitude plus elevee pour 
celle du haut). La partie de droite de la figure nous 
montre le resultat du prolongement analytique. 



Figure 12: Fonction de Green en temps imaginaire 
pour des systemes de differentes tailles. On presente 
les donnees pour le cas U = —4, n = 1, k = (0,7r) 
avec pour le MCQ At = 1/10. La partie de gauche 
montre les donnees MCQ tandis que la partie de 
droite a ete obtenue suite l'utilisation des series de 
Fourier sur l'approche analytique evaluee en frequence 
de Matsubara. 
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Figure 13: Poids spectral pour des systemes de 
differentes tailles et au vecteur d'onde k= (0, 7r). La 
partie de gauche de la figure montre les resultats 
MCQ et la partie de droite les resultats de 
l'approche analytique (TPSC). Ces resultats viennent 
du cas U = —4, n = 1 avec pour les donnees MCQ 
At = 1/10. 



Figure 14: Contribution de la susceptibilite de paires 
uniforme a frequence nulle au facteur de structure de 
paires en fonction de la temperature. La valeur 1 
correspond au cas oil le facteur de structure serait 
totalement determine par cette frequence. Les points 
proviennent des mesures MCQ pour le cas U = —4, 
n = 0.95 et sont presentes pour differentes tailles de 
reseaux. 



Figure 15: Partie imaginaire de la susceptibilite de 
paires uniforme divisee par la frequence et en fonction 
de la frequence. Les resultats proviennent du modele 
de Hubbard avec U = —4, n = 1 sur un reseau carre 
8x8 avec une discretisation en temps imaginaire de 
At = 1/10. On presente les resultats aux temperatures 
T = 1/3, 1/4 et 1/5. 
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Figure 16: Fonction de correlation courant-courant 
dans la direction x pour differents vecteurs d'onde. 
Les carres representent q y = 0, iq n = avec differentes 
valeurs de q x , les losanges representent le resultat 
pour differentes valeurs de q y avec q x = et iq n =0 
et les cercles represente le cas limite iq n — > avec 
q x = q = 0. Les resultats presentes sont pour les 
conditions habituelles, U = —4, n = 1, At = 1/10 pour 
un reseau 8x8 a temperature T = 1/4. 



Figure 17: Spectre a une particule pour un reseau 
de taille 4x4 a basse temperature T = 1/8 et au 
vecteur d'onde (0, 7r). En comparaison, on montre le 
spectre obtenu a T = 1/5. Les valeur des parametres 
du modele sont U = —4, n = 1. avec une discretisation 
du temps imaginaire de At = 1/10. 



Figure 18: Longueur de correlation de paires 
comparee a la longueur d'onde de de Broglie en 
fonction de la temperature. Les donnees sont extraites 
des mesures MCQ pour un cas a interaction U = —4, 
et remplissage n = 0.95. Les losanges representent la 
longueur de correlation pour un reseau de taille 6x6, 
les triangles, celle d'un reseau 8x8, les croix celle 
extrapolee pour la limite thermodynamique et les 
cercles pleins, la longueur d'onde de de Broglie pour 
un reseau 8x8. 
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Figure 19: Comparaison entre le facteur de structure 
de paires (symboles gras) et le facteur de structure de 
charge (symboles vides) en fonction de la temperature. 
Les carres presentent les resultats pour le reseau 6x6 
et les triangles pour le reseau 8x8. La partie de 
gauche montre les resultats au remplissage n = 0.95 
et celle de droite pour n = 0.8. Pour les deux 
remplissages nous avons utilise une interaction U = — 4 
et une discretisation de At = 1/10. 



Figure 20: Poids spectral au remplissage n = 0.95 et 
au vecteur d'onde k = (0, n) et au remplissage n = 0.8 
pour le vecteur k= (0,37r/4). Les deux courbes ont 
ete obtenues par 1' application de la technique de 
maximisation d'entropie aux donnees MCQ. Nous 
presentons les resultats a temperature T = 1/5 avec 
interaction U = —4. 
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SOMMAIRE 



Dans cette these nous presentons une nouvelle methode non perturbative 
pour le calcul des proprietes d'un systeme de fermions. Notre methode generalise 
l'approximation auto-coherente a deux particules proposee par Vilk et Tremblay 
pour le modele de Hubbard repulsif. Notre methode peut s'appliquer a l'etude 
du comportement pre-critique lorsque la symetrie du parametre d'ordre est 
suffisamment elevee. Nous appliquons la methode au probleme du pseudogap dans 
le modele de Hubbard attractif. Nos resultats montrent un excellent accord avec 
les donnees Monte Carlo pour de petits systemes. Nous observons que le regime 
ou apparait le pseudogap dans le poids spectral a une particule est un regime 
classique renormalise caracterise par une frequence caracteristique des fluctuations 
supraconductrices inferieure a la temperature. Une autre caracteristique est la 
faible densite de superfluide de cette phase demontrant que nous ne sommes pas 
en presence de paires preformees. 

Les resultats obtenus semblent montrer que la haute symetrie du parametre 
d'ordre et la bidimensionalite du systeme etudie elargissent le domaine de 
temperature pour lequel le regime pseudogap est observe. Nous argumentons que 
ce resultat est transposable aux supraconducteurs a haute temperature critique 
ou le pseudogap apparait a des temperatures beaucoup plus grandes que la 
temperature critique. La forte symetrie dans ces systemes pourraient etre reliee a 
la theorie 50(5) de Zhang. 

En annexe, nous demontrons un resultat tout recent qui permettrait 
d'assurer l'auto-coherence entre les proprietes a une et a deux particules par 
l'ajout d'une dynamique au vertex irreductible. Cet ajout laisse entrevoir la 
possibilite d'etendre la methode au cas d'une forte interaction. 
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